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PHOSPHORYLATIVE AND NONPHOSPHORYLATIVE 
PATHWAYS OF ELECTRON TRANSFER IN 
RAT LIVER MITOCHONDRIA* 


By NATHAN ©. Kaptan, Morton N. Swartz,t Mary E. Frecu, snp MARGARET 
M. Crorri 
McCOLLUM-PRATT INSTITUTE, JOHNS HOPKINS UNIVERSITY, BALTIMORE, MARYLAND 
Communicated by B. H. Willier, June 15, 1956 


Separate pathways for the oxidation of reduced diphosphopyridine nucleotide 
(DPNH) and reduced triphosphopyridine nucleotide (TPNH) are known to exist in 
animal tissues. Although phosphorylation has been shown to be associated with 
oxidation of DPNH, the results with TPNH have been variously interpreted.! 
Hunter and Ford? have described rat liver mitochondrial preparations in which the 
oxidative and phosphorylative systems could be reversibly inactivated. Almost 
complete restoration of oxidative phosphorylation occurred upon the addition of 
diphosphopyridine nucleotide (DPN) and an appropriate substrate. We have 
found that such mitochondrial preparations do not contain either the oxidized or the 
reduced forms of DPN, but may contain traces of either TPNH or its oxidized form, 
TPN. It has been possible to study the single oxidative step from isocitrate to 
alpha ketoglutarate in these pyridine nucleotide—“depleted”’ mitochondria, by the 
addition of TPN. Using this approach, the question of phosphorylation associated 
with the oxidation of intramitochondrially generated TPNH has been investigated. 
In the course of this study, information has also been obtained which may be of 
significance in elucidating the function of TPNH and DPNH systems in regulating 
cellular respiration. 

MATERIAL AND METHODS 

Mitochondrial Preparation.—Rat liver mitochondria are ground gently with 6 
volumes of 0.25 M sucrose + 0.02 M tris(hydroxymethyl)aminomethane (pH 7.5) 
in a loose-fitting glass homogenizer. The homogenate is centrifuged three times at 
600 X g, for periods of 5 minutes each, to remove easily sedimentable materials. 
The homogenate supernatant is then centrifuged at 8,000 x g for 20 minutes and 
the mitochondrial pellet washed twice with 0.25 M sucrose + 0.02 M tris(hydroxy- 
methyl)aminomethane (pH 7.5). Final suspension in the same medium was made 
so that 1 ml. represents 0.6 gm. of liver. 

Depletion of Mitochondria._To 2.4 ml. of the above mitochondrial suspension is 
added 0.6 ml. of 0.25 M sucrose + 0.1 M potassium phosphate (pH 7.5). The 
mixtures are incubated by shaking in a tube for 15 minutes at 30° and then returned 
immediately to ice, prior to use in the Warburg experiments. 
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Warburg Flask Additions.The medium in the Warburg flasks contains the fol- 
lowing substances in amounts to give the indicated final concentration in the 3.0-ml. 
total mixture: 0.0017 .W ATP, 0.02 W NaF, 0.005 WM MgCl, 0.01 47 KCI, 0.005 / 
glyeylglycine (pH 7.4), 0.001 7 Na versenate (pH 7.5), potassium phosphate 
(pH 7.4), and 0.007 per cent cytochrome ¢ (Mann Research Laboratories). The 
reaction mixtures are kept cold, and the glucose, 0.02 M final concentration, and 
360 Cori units of yeast hexokinase are tipped in just as the flasks are placed in the 
30° bath. Other additions are indicated in the tables. After equilibration for 
10 minutes, the vents are closed. Oxygen uptake during the equilibration period is 
calculated from the 7-minute period after closing the vents, in a manner described 
by Hunter and Ford.2- Oxygen uptake is measured for 30-40 minutes, and upon 
removal of the mixtures from the manometer, 1.5 ml. of 18 per cent trichloroacetic 
acid is added as quickly as possible. Phosphate and alpha ketoglutarate are de- 
termined in the trichloroacetic acid extract. The alpha ketoglutarate is deter- 
mined by the method of Friedemann.* 


RESULTS 


Comparison of Phosphorylation Associated with the Oxidation of Mitochondrial- 
generated TPNH and DPNH.—Table | summarizes a typical experiment, comparing 


TABLE 1 


OXIDATIVE PHOSPHORYLATION IN UNDEPLETED AND 
DEPLETED MITOCHONDRIA 


DPN and TPN additions: 1.5 « moles per vessel 
Isocitrate and glutamate additions: 30 «4 moles per vessel 
Mitochondria addition: 1.7 ml. of final suspension per vessel 
Oxygen 
Uptake Phosphorus P:O 
Additions (u atoms) Uptake (u moles) Ratio 


Undepleted* mitochondria +glutamate 3.5 7.62 2.2 

Undepleted mitochondria + isocitrate 9.95 29.9 3 

Depletedt mitochondria +glutamate 0 0 

Depleted mitochondria +isocitrate 0 0 

Depleted mitochondria + DPN +glutamate 5.3 15.0 2.8 

Depleted mitochondria +TPN +isocitrate 3.0! 0.45 0.14 
* DPN content of undepleted mitochondria per Warburg flask is 20 wu 


+ DPN content of depleted mitochondria per Warburg flask is less an 1 wg., mitochondria having been pre- 
incubated with phosphate to remove pyridine nucleotides. 


the relative activity of TPN and DPN in promoting phosphorylation coupled to 
oxidation. As can be seen from the data, in the normal, ‘“undepleted”’ mitochon- 
dria, both isocitrate and glutamate oxidation result in phosphorylation. The 
“depleted’”’ mitochondria still possess the capacity to promote phosphate uptake 
when glutamate is oxidized. Isocitrate oxidation, however, does not result in 
phosphorylation. Alpha ketoglutarate accumulates almost ‘quantitatively as a 
result of the oxidation of isocitrate in the ‘depleted’ particles. In a number of 
experiments with the ‘‘depleted”’ mitochondria, the addition of TPN and isocitrate 
results in a P:O ratio of zero or close to zero. When extremely large amounts of 
TPN are used, a ratio as high as 0.4 has been obtained. This mitochondrial prep- 
aration can convert TPN and DPN at a slow rate, and the slight phosphorylation 
observed with the high levels of TPN may be the result of the formation of small 
amounts of DPN. 
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The addition of DPN, as well as TPN, and isocitrate to the ‘depleted’? mito- 
chondria resulted in phosphorylation (Table 2). This phosphorylation was as- 


TABLE 2 


RESTORATION OF PHOSPHORYLATION WITH ISOCITRATE 
OXIDATION BY DPN* 
a-Keto- 
Oxygen Phosphorus glutarate 
Uptake Uptake P:O Formed 
Additions w atoms) u moles) Ratio (« moles) 


Depleted mitochondria + DPN +glutamate 1.4 11.4 2.7 0.2 
Depleted mitochondria+TPN +isocitrate 3.2 0 0 2.4 
Depleted mitochondria +TPN + DPN +isocitrate 9.6 22.2 2.3 $3 


* Conditions the same as for Table 1 


sociated with an increase in oxygen uptake and in alpha ketoglutarate accumula- 
tion. The fact that phosphorylation is restored on the addition of DPN rules out 
the possibility that TPN is an inhibitor of oxidative phosphorylation. 

It is of interest to note that oxygen uptake is much greater with isocitrate than 
with glutamate in the untreated mitochondria (Table 1). In the “depleted”’ 
system the rate of glutamate oxidation is somewhat faster than isocitrate. How- 
ever, as indicated in Table 2, the addition of DPN to the TPN-isocitrate system 
results in a rate of oxygen uptake which is considerably faster than with glutamate 
and DPN. 

Pyridine Nucleotide Transhydrogenase Activity of Mitochondria.—The increase in 
respiration, resulting from the addition of DPN in the presence of isocitrate and 
TPN, suggested the possibility that pyridine nucleotide transhydrogenase might be 
involved in the stimulation. This enzyme, which has been found in animal tissues, 
catalyzes the following reversible reaction: 


TPNH + DPN = DPNH + TPN. 


Evidence for the role of the transhydrogenase in liver mitochondria is given in 
studies of cytochrome ¢ reduction, as illustrated in Figure 1. The rate with DPN 
and glutamate was found to be quite rapid when compared with the rate with TPN 
and isocitrate. TPN addition to the DPN-glutamate system actually slowed the 
reaction (Fig. 1). Addition of both DPN and TPN to the isocitrate reaction 
mixture produced a marked stimulation of cytochrome ¢ reduction. It is of in- 
terest that the rate of cytochrome ¢ reduction in the presence of isocitrate and DPN 
is somewhat faster than with isocitrate and TPN. This may be due to traces of 
TPN in the depleted mitochondria (see below). 

The data in Figure | indicate that there is a fast reduction of cytochrome ¢ by 
generated DPNH and a slow reduction by TPNH generated by isocitrate oxidation. 
The effect of the addition of DPN to the isocitrate reaction mixtures suggests a 
transfer of electrons from TPNH to DPN via the transhydrogenase with formation 
of DPNH, which could then be more efficiently oxidized by the cytochrome ec. 
The slower rate observed when TPN is added to the DPN-glutamate system can 
also be attributed to transhydrogenase activity, since there is more TPN than 
cytochrome ¢ present in the reaction mixture in Figure 1. 

Additional support for the above view has been obtained from studies on the 
reduction of the dye 2,6-dichlorophenolindophenol (Fig. 2). In contrast to cyto- 
chrome ¢ reduction, the rate of dye reduction with isocitrate and TPN is faster 
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Fic. 1.—Pyridine nucleotide transhydrogenase activity of rat liver mitochondria as indicated by 
cytochrome ec reduction. All reaction mixtures, 3.0 ml. total volumes, contained 2 mg. cyto- 
chrome c, 0.1 MW phosphate buffer (pH 7.5), 0.001 4 KCN, 0.003 1 MgCh, 0.05 ml. “depleted” 
mitochondria, and 500 wg. DPN or TPN, except when both nucleotides were added: 1,500 ug. 
TPN with glutamate and 1,500 ug. DPN with isocitrate were used. The reactions were started 
with 5 u moles of sodium glutamate or with 10 « moles of d/-isocitrate. 

















than with glutamate and DPN. The addition of TPN to the DPN-glutamate re- 
action and of DPN to the TPN-isocitrate system has no effect on dye reduction. 
These experiments indicate that the TPN-isocitrie dehydrogenase is a highly ac- 
tive system in mitochondria and that the limiting factor in the rate of isocitrate 
oxidation (in the absence of DPN) is the slow TPNH-—cytochrome ¢ system. More 
detailed experiments, directly showing the transhydrogenase activity, will be pub- 
lished elsewhere. 

DPN Isocitric Dehydrogenase.—The stimulation of the oxidation of isocitrate by 
DPN might be the result of a DPN-specific isocitric dehydrogenase rather than the 
transhydrogenase. Indeed, the addition of DPN in the absence of added TPN to 
the “depleted” mitochondria occasionally promoted the oxidation of isocitrate. 
It is our belief, however, that the presence of residual TPN and the transhydrogen- 
ase are responsible for this oxygen uptake. 

Plaut and Sung? report that rat liver mitochondria contain only traces, if any, 
of DPN isocitric dehydrogenase. In Figure 2 data are given which show that 
dye reduction with isocitrate and DPN is very slow. Furthermore, we have shown 
by direct test that the ‘depleted’ mitochondria do not reduce DPN when isocitrate 
is the substrate; under the same conditions TPN is rapidly reduced. 

These findings strongly indicate that DPN isocitric dehydrogenase is not involved 
in stimulation of oxygen uptake by DPN when isocitrate is the substrate. 
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Fig. 2.—Dye reduction by depleted rat liver mitochondria. Reaction mixtures were the 
same as those described for Figure 1, with the omission of the cytochrome ¢ and the substitu- 
tion of 49 yg. 2,6-dichlorophenolindophenol. 


DISCUSSION 
From the results presented in this paper, it is evident that TPNH and DPNH 
oxidation proceed through separate pathways in rat liver mitochondria. DDPNH 


oxidation through the cytochrome ¢ system is a rapid one and leads to phosphoryla- 
tion. On the other hand, the flow of electrons from TPNH to oxygen is a much 


more sluggish system, as indicated by the slow rate of cytochrome ¢ reduction. 
The TPNH pathway does not appear to be associated with the generation of 
energy-rich phosphate bonds. It cannot be stated with certainty that there is a 
complete absence of phosphate uptake in the oxidation of TPNH, since P:O ratios 
of 0.1 to 0.4 have been observed with some mitochondrial preparations. Whether 
or not these ratios are significant remains for further study. Since a number of 
preparations have yielded no phosphate uptake, it is possible that in the case of 
TPNH oxidation no phosphorylation occurs even between cytochrome ¢ and oxy- 
gen. The results we have obtained with TPNH may support the concept of 
Green et al.6 and Reif and Potter’? that separate electron-transferring systems 
may exist within the mitochondria. 

Figure 3 presents a scheme which illustrates the role of pyridine nucleotide trans- 
hydrogenase in linking the TPN and DPN systems. Phosphorylation may be 
obtained with TPNH oxidation only after electron transfer to DPN via transhydro- 
genase. Our data suggest that the phosphate uptake associated with TPNH 
oxidation may occur through such a mechanism (see Table 2). The transhydro- 
genase, by shifting electrons from TPNH to DPN and thence to cytochrome c¢, 
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will stimulate the rate of oxygen uptake, since the DPNH-cytochrome c¢ electron- 
transport system is much faster than the corresponding TPNH reaction. 

Recently Jacobson,’ of our laboratory, has shown that TPNH is the major pyri- 
dine nucleotide of carefully prepared rat liver mitochondria. Little DPNH is 
present in these mitochondria. The survival of TPNH during the isolation of 
mitochondria might be expected because of the relatively slow rate of oxidation 
of TPNH as compared with the oxidation of DPNH. 

Langdon’ has shown that TPNH can act as the reducing agent in the conversion 
of crotonyl CoA to butyryl CoA; DPNH cannot replace TPNH in this reaction, 
which leads to the synthesis of fatty acids (see Fig. 3). It is possible that the 
prime function of TPNH is to serve as a reducing agent. DPNH, on the other 
hand, is primarily oxidized and thereby gives rise to utilizable energy in the form 
of energy-rich phosphate bonds. The pyridine coenzymes may have different major 
functions. These different functions would then explain the required presence of 
two pyridine nucleotides rather than one in living systems. In this connection, 
it is of interest that Hoch and Lipmann" indicated the significance of the dual func- 

tion of reducing systems arising from 
Substrate — DPNH FAST Q the oxidation of substrates for both the 
a generation of energy-rich phosphate 
bonds and also for synthetic purposes 

as in fatty acid synthesis. 
Substrate — TPNH SLOW 0 The pyridine nucleotide transhydro- 
NO ~P genase has previously been postulated 
to serve as a regulatory mechanism for 
the conservation of oxidative energy.!! 
As illustrated in Figure 3, the shifting 
FATTY ACIDS of electrons from DPNH to the TPN 
Fic. 3.—-Pathways of electron transport. system would result in a lowering of 
phosphorylation; a shift in the direc- 
tion of the DPN system would increase the level of phosphorylation. Work is 


now in progress in this laboratory on possible factors which may regulate the 


activity of the transhydrogenase system. 


SUMMARY 


By removing pyridine nucleotides from rat liver mitochondria and then adding 
TPN and isocitrate, it has been possible to show that TPNH oxidation does not 
result in phosphorylation. 

The oxidation of DPNH (generated from the oxidation of glutamate) in the same 
preparation results in phosphorylation. Cytochrome ¢ reduction by mitochon- 
drially generated DPNH has been found to take place at a faster rate than with 
TPNH. The role of pyridine nucleotide transhydrogenase in rat liver mitochon- 
dria has been investigated, and a scheme that involves this enzyme as a regulator 
of respiration and of energy-releasing reactions is presented. 

* Contribution No. 146 of the MeCollum-Pratt Institute, aided by grants from the National Can- 
cer Institute of the National Institutes of Health (grant No. C-2374C), the American Cancer Soci- 
ety as recommended by the Committee on Growth of the National Research Council, and the Na 
tional Science Foundation. 
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t Postdoctoral fellow, National Cancer Institute, National Institutes of Health. Present ad- 
dress, Department of Medicine, Massachusetts General Hospital, Boston. 

1S. Ochoa, Ann. N.Y. Acad. Sci., 47, 835, 1947. 

* F. E. Hunter, Jr., and L. Ford, J. Biol. Chem., 216, 357, 1955. 

°'T. Friedemann, in 8. P. Colowick and N. O. Kaplan (eds.), Methods in Enzymology, Vol. 3 
(New York: Academic Press) (in press). 

''N. O. Kaplan, 8. P. Colowick, and E. F. Neufeld, J. Biol. Chem., 205, 1, 1953. 

5G. W. E. Plaut and 8.-C. Sung, /. Biol. Chem., 207, 305, 1954. 

®°D. E. Green, ef al., in W. D. McElroy and B. Glass (eds.), Inorganic Nitrogen Metabolism 
(Baltimore: Johns Hopkins University Press, 1956), p. 628. 

7A. KE. Reif and V. R. Potter, Arch. Biochem. and Biophys., 48, 1, 1954 

* K. B. Jacobson, dissertation submitted to the Johns Hopkins University (to be published ). 

*R. G. Langdon, J. Am. Chem. Soc., 77, 5190, 1955. 

 F. L. Hoch and F. Lipmann, these PRocEEDINGs, 40, 909, 1954. 

'N. O. Kaplan, Bacteriol. Revs., 19, 235, 1955. 

'2 8. P. Colowick, in W. D. McElroy and B. Glass (eds.), Phosphorus Metabolism, 1 (Baltimore: 
Johns Hopkins University Press, 1952), 297. 


CONDUCTANCE, VISCOSITY, AND VAPOR DENSITY FOR 
SEVERAL SALTS IN BROMINE AT 25° 
By Puitie L. Mercter* and Cuarzes A. Kraus 
METCALF RESEARCH LABORATORY, BROWN UNIVERSITY 
Communicated May 29, 1956 

Equivalent conductances reported by Moessen and Kraus! for several salts in 
bromine, and viscosities reported by Rice and Kraus? for the same salts in the same 
solvent, show that, at high concentrations, salts are largely or even completely dis- 


sociated into their ions. It seemed worth while to study such solutions in some- 


what greater detail. The chlorides and bromides form with one molecule of bro- 
mine rather stable complexes which melt at comparatively low temperatures. 
Compounds of the type R;NHX-Br. melt in the neighborhood of 35° and are 
readily handled at 25° in a supercooled condition. Here R represents an alkyl 
radical and XY a halogen atom. 

We have measured the conductance and the viscosity of solutions of Am;NHCI- 
Bro and Me;NHBr-Br. from comparatively low concentrations to the molten salts 
at 25°, and we have likewise determined the vapor densities over solutions of these 
two salts as well as of BusNBr-Bro and Me;NHCI-Bre. Mercier and Kraus* have 
determined the phase relations for these systems. 

Bromine was purified by the method described by Moessen and Kraus. The 
salts were prepared according to conventional methods. Bu,NBr was recrystal- 
lized from ethyl acetate, m.p. 117°-117.5°. Conductances of this salt in water 
closely checked those of Martel* with another preparation. Me;NHCl was re- 
crystallized from ethanol-methanol mixtures, Me;NHBr from acetone-methanol 
mixtures, and 7-AmNHCI from acetone-ether mixtures. 

Viscosities were measured by means of a modified Ostwald viscometer. The 
viscometer was provided with stopcocks, and the receiving bulb was provided with 
a l-cm. tube which was capped. Salt could be introduced through this tube. 
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Concentrated solutions, as well as the supercooled fused salts, were made up in a 
special preparation tube from which they were transferred to the viscometer. 
These solutions were diluted by adding liquid bromine. Dilute solutions were 
made up by adding weighed salt samples to a known quantity of bromine in the 
viscometer. The bromine content of the viscometer was determined by weighing 


on a sensitive balance. 
Four viscometers were constructed; they were calibrated against water and 
known sucrose solutions. Their efflux times are given in Table 1. The time 


TABLE | 
ErrLux TIMEs OF VISCOMETERS 
Efflux 


Viscometer Liquid Time 
Number Liquid Volume (Sec.) 


| Water 11.87 895.0 
2 40°7 sucrose 9.98 720. 
3 60°, sucrose 12.60 988 

i 60°, sucrose 10.34 140 


TABLE 2 


DENSITIES OF SOLUTIONS OF SALTS IN 
BROMINE AT 25° 
A. MesNHBr- Bre B. Am sNHCI- Bre 

Mole Fraction Density Mole Fraction Density 
of Salt (Gm/MI) of Salt (Gm/MI) 
0 3.102 0 3.102 
0.05347 3.020 04128 2.717 
07194 2.991 05050 643 
1026 2.941 1151 283 
.2018 2.778 1399 187 
2264 2.735 2909 777 
2621 2.680 3928 610 
3711 2.523 4452 555 
4336 2.446 5300 472 
6521 2.243 6390 398 
7882 2.149 7963 318 
0.9214 2.075 8568 295 
0.9804 250 

1.006 244 


Pe et pt pat pe ees NOD 


of efflux was determined as a function of the volume of solution in the viscometers, 
so that it was not necessary to adjust the volume of the solution to a definite value. 
A plot of volume versus efflux time made it possible to obtain the efflux times of the 
standard solution for any degree of filling. 

For Me;NHBr, only viscometers 1 and 2 were required; for Am;NHCI all four 
viscometers were employed. In going from one viscometer to the next, concen- 
trations were adjusted to overlap, thus precluding the possibility that the two 
yielded different results. 

The more concentrated solutions were made up by distilling bromine onto a 
known weight of salt in the viscometer. The melting points of highly concentrated 
solutions are above 25°, but, when melted, the solutions may readily be handled in 
a supercooled condition. After having determined the viscosity of the most con- 
centrated solution, more bromine was added and the viscosity of the solution deter- 
mined. This process was repeated until the limiting capacity of the viscometer 
was reached. 

For the more dilute solutions, known weights of salt were added to a known 
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weight of bromine until a desired concentration was reached. Viscosity values are 
shown in the third column of Table 3 for rounded mole fractions of salt which is 
assumed to be complexed with one molecule of bromine. 


TABLE 3 
CONDUCTANCES, VISCOSITIES, AND VAPOR DENSITIES OF SOLUTIONS OF SALTS IN BROMINE 
COoNDUCTANCE- Rev. Vapor- 

MoLe ViscosiTy VAPOR DENSITY 

FRACTION CONDUCTANCE Viscosity PrRopwuct DensiT? LOWERING 
Me;sNHBr 

0.0125 0 

O14 0.22 900 0.008 

O16 O:3 877 020 

835 042 

0215 0 é 769 076 
490 . 222 
327 307 
166 391 
O11 472 
474 753 
224 883 
124 935 
059 969 
050 974 
034 982 
023 988 
O17 991 


0.02 
0.04 
0.06 
0.08 
10 
14 
16 
18 
20 
30 
10 
50 
60 
70 
80 
90 
0 


i 
a te — 
oN Oo 


Pr 


0330 


a a) 


0418 
0639 
OS48 
106 
129 
153 
183 
215 
254 
Am sNHCI 


HED AAM- OD 
—-nNnNwnnnn--— 


40 

40 

25 

75 0093 035 

04 909 003 
26 ~ 
26 0305 130 : 

0 95 0481 190 735 094 
0 80 104 290 406 266 
0 331 305 
0 73 205 355 005 475 
0.25 05 386 105 622 675 
0 680 654 445 391 796 
0.3: .470 

0 160 998 459 184 .904 
0 458 

0.42 470 ! =} 

0. 0.245 00 490 090 953 
0 0.190 59 492 052 973 
0 0.140 32 465 O31 .984 
0 0.110 05 445 019 .990 
0 0.090 83 435 0.015 .992 
1.0 0.075 80 0.435 0.013 0.993 


0 
0 
0 
0 
0 
0 


1 
2 
3 
3 

4 

' 
4 
3 
2 


Oo — —_ 
> 


~ 
> 


Ot ie de Oo NO DO 


Densities of solutions were determined by means of a pycnometer. Values are 
shown in Table 2 for different mole fractions of salt which is assumed to contain 
one molecule of bromine. Vapor densities were determined by weighing a bulb of 
about 350-ce. capacity, filled with vapor which was in equilibrium with a known 
solution. These solutions were prepared in a cylindrical tube of about 100-ce. 
capacity which was provided with a stopcock and a small capped tube through 
which salt could be introduced. Bromine was distilled into this tube under vacuum 
conditions. After distillation, the tube was warmed and shaken until its contents 
became homogeneous. 
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The vapor-density bulb was then attached to the tube containing the solution 
and exhausted. Then the solution was cooled to dry-ice temperature, connection 
was made with the vapor-density bulb, and the whole apparatus was thoroughly 
exhausted in order to remove all traces of air. The solution was then brought to 
25° and maintained there with occasional shaking until equilibrium was estab- 
lished. The vapor-density tube was then detached and weighed. As was shown 
by measurements of vapor densities at different pressures of bromine, the vapor 
density, at the low pressures of the present experiments, is a linear function of 
pressure. Accordingly, the relative vapor-density lowering (do — d)/dy is equal 
to the vapor-pressure lowering (py — p)/po. The density of saturated bromine 
vapor at 25° was found to be 1.915 gm, liter. 

By weighing under hexane, the densities of Am;NHCI and Me;NHBr were found 
to be 0.933 and 1.570, respectively. For BusNBr and Me;NHCI we adopted the 
density values of Rice and Kraus, 1.225 and 1.099, respectively. 

The conductance cell consisted of two cylindrical tubes of about 50-mm. diam- 
eter which were connected at the bottom with a 10-cm. length of tube about 7 mm. 
in diameter. The cylindrical tubes were provided at the top with ground-glass 
stoppers which carried the electrodes at the bottom of the two large tubes. The 
cell was provided with suitable stopcocks and could be attached to a source of 
bromine and the vacuum system. The cell had a constant of approximately 36, 
which was independent of the position of the electrodes or the volume of solution in 
the cell. Means were provided for stirring the solution. Concentrated solutions 
were made up by adding bromine to a known weight of salt; dilute solutions, by 
adding salt to a known weight of bromine. The cell contents were determined 
by weighing on a sensitive balance. 

The results of our measurements of conductance, viscosity, and vapor-density 
lowering for Am;NHCI and Me;NHBr are shown in Table 3, where concentration 
is expressed as mole fraction of salt (containing one molecule of bromine) at rounded 
values as interpolated from large-scale plots. Values of equivalent conductances 
are given in the second column, of viscosity in the third column, and of the con- 
ductance-viscosity product, An, in the fourth column. Vapor densities are given 
in the fifth column and relative vapor-pressure lewerings in the sixth column. 
Vapor densities and vapor-density lowerings for BusNBr and Me;NHCI are given 
in Table 4. 

In Figure | values of A are shown plotted as a function of mole fraction, x, of salt 
assumed to be Me;NHBr-Bre and Am;NHCI-Bro, respectively. Points shown on 
these, as on all other plots, are experimental, not interpolated. 

The conductances of the bromide and the chloride pass through maxima of 37.0 
and 4.26 at x = 0.11 and 0.075, respectively. For the fused salts at x = 1, the 
values of A are 7.45 and 0.075, respectively. It is of interest to note that the con- 
ductances of Me;NHCI-Br. and BuyNBr-Bre pass through maxima of 15.75 and 13.9 
at x = 0.11 and 0.08, respectively.! These salts are all rather strong electrolytes 
in bromine solution. It will also be noted that while the conductance of Am;N HCI 
is markedly lower than that of Me;NHBr at the maximum, that of Bu,NBr is almost 
as large as that of Me;NHCI. The reason for this, as we will see, is due to the 
much lower conductance of the chloride with respect to that of the bromide ion. 

The viscosities of solutions of Me;NHBr and Am;NHCI are shown graphically 
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in Figure 2. Both plots show a slight inflection in the concentrated regions. The 
viscosities of the fused Me;NHBr-Bro and Am;NHCI-Br2 (2 = 1) are, respectively, 
27 and 620 times that of bromine (9.435 & 10~* poise). The high viscosity of the 
chloride is due to the presence of the large tri-?soamylammonium ion. At a con- 
centration of 1.0 molar, the relative viscosities of solutions of Me;NHCI1, Me;NHBr, 
Am;NHCI, and BuyNBr are, respectively, 1.6, 1.7, 3.4, and 3.9. 


TABLE 4 


Vapor DENSITIES OF SOLUTIONS OF SALTS 
IN BROMINE 
A MesNHCI B. BuiNBr 
Rel. Vapor- Rel. Vapor- 
Mole Vapor Density Mole Vapor Density 
Fraction Density Lowering Fraction Density Lowering 


0.04 é 0.04 882 0.017 
0.06 0.008 06 819 .050 
0.08 O8 727 098 
0.10 842 O38 10 599 165 
0.14 676 125 14 254 345 
0.16 515 209 16 048 453 
0.18 316 313 18 837 563 
0.20 160 394 0.20 674 0.648 
0.30 651 660 
0.40 393 795 
0.50 236 877 
0.60 149 922 
0.70 100 948 
0.80 071 963 
0.90 O57 970 
1.0 0.048 0.975 


eet et 


— 
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The conductance maximum and the rapid decrease of conductance at higher 
concentrations are due to the increasing viscosity of the solutions. At low concen- 
trations, where the viscosity of the solutions differs little from that of the pure 


solvent, we are able to approximate the degree of ionization of the electrolyte by 


means of the equivalent conductance of the solution, for the values of Ag may 
always be approximated by means of Walden’s rule. At higher concentrations, 
where the viscosity of the solution becomes much greater than that of the pure 
solvent, this can no longer be done. However, the degree of ionization may again 
be approximated by means of Walden’s rule. According to this rule, which is quite 
generally applicable, the conductance of the solution is proportional to the reciprocal 
of the viscosity. Thus the change in the product, An, yields an approximation of 
the change in the degree of ionization. If ionization is complete, the conductance- 
viscosity product should remain substantially constant or vice versa. If the 
product remains constant on changing the condition of the system, the electrolyte 
may be considered to be completely dissociated into its ions. If the product de- 
creases, the degree of ionization must be decreasing and vice versa. 

In Figure 3 are shown plots of the Ay products for Am;NHCI-Bre and Me;- 
NHBr-Bro, as a function of mole fraction. For these salts as supercooled liquids, 
An has the values 0.470 and 1.89, respectively. For Bu,NCl in ethylene chloride, 
Aon = 0.512. The value of An for Am;NHCI-Br, in the liquid state is much the 
same as that for similar salts in ethylene chloride at infinite dilution. This accords 
with the assumption that in the fused state the electrolyte is completely dissociated 
into its ions. 

The value of An for Me;NHBr-Br, is unexpectedly high; for Me,N Br in ethylene 
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chloride, Aon = 0.600. It seems that the bromide ion is a much better conductor 
than the chloride ion in these fused salts as well as in their solutions in bromine. 
We do not have data on the conductance and viscosity of Me;NHCI in the fused 
state, but from the data of Rice and Kraus? it appears that An approaches a value 
of about 1.0 in the limit. So, also, the conductance of BuyNBr at about 1 molar, 
the highest concentration measured, is practically as large as that of Me;NHCI, 
although the viscosity of the bromide solution is several times greater than that 


of the chloride. 
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Fic. 1.—Equivalent conductances of salts in bromine. 
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If Walden’s rule holds even roughly, it follows that in the fused salts, as well as 
in their solutions in bromine, the complexed bromide ion has about four times the 
conductance of the corresponding chloride ion. In the fused salt the bromide ion 
consists of the bromine atom with a single charge for the whole. The ion may well 
have a linear structure. In any case, any one of the bromine atoms may carry the 
charge at any given instant. The ion thus behaves much like the hydrogen, or, 
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better, the hydroxyl, ion in water. On the other hand, in the chloride ion, Cl-Bre, 
the charge is definitely localized in the chlorine atom; there is no exchange of 
charge with the adjoining bromine atoms. Accordingly, the chloride ion has a 
normal conductance. 

It has been shown on the basis of phase studies that bromide and chloride ions 
form a series of complexes with bromine molecules.* In some cases these com- 
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plexed ions contain as many as four bromine molecules. As bromine is added to 
the fused salt, it goes, for the most part, to form complexes, and but little remains 
as free bromine in solution. This is clearly shown by the very low vapor pressure 
of the concentrated solution, as illustrated in Figure 4. Up to a mole fraction of 
about 0.4, practically all the bromine is combined with the halide ion. Thus we 
are not dealing with a single fused salt; rather, we are dealing with a mixture of 
salts of varying composition. This accounts for the fact that the An product does 
not remain constant or decrease on initial addition of bromine to the fused salt. 
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If we know the value of the conductance-viscosity product of the salt when com- 
pletely ionized (we shall designate this as (An)o), We can approximate the degree of 
its association. The degree of ionization is given by the expression 


y = An/(An)o 


and the degree of association is | — y. The error in the evaluation of 1 — y by 
this method will be large when | — y is small. However, when | — y is large, a 
fair approximation may be made. 

We cannot determine just when ion association sets in. For the present calcu- 
lation, we shall assume that it sets in at x = 0.4 for Me;NHBr and at 2 = 0.35 for 
Am;NHCl. The corresponding values of (An)o are, respectively, 2.12 and 0.470. 

In Table 5 are shown values of | — y for the two salts computed on the basis of 
the method just outlined. Ion association sets in as soon as the concentration of 


TABLE 5 
DEGREE OF ASSOCIATION, | — y, OF SALTS 
IN BROMINE 
A. MesNHBr B. Am:sNHCI 
(An)o = 2.12 atz = 0.4 (An)) = 0.47 atx = 0.35 
zx 1-— ¥ £ 1-4 
} 0.00 0.35 0.00 
3 095 30 06 
2 30 25 14 
16 42 20 25 
.10 66 15 38 
06 8&2 .10 
04 90 O08 
02 0.96 0.04 


0 


free bromine becomes significant. Trimethylammonium bromide complexes with 
a total of three molecules of bromine or two molecules above that of the fused salt, 
R,NHX-Bro. At x = 0.33, there are two molecules of bromine above that of the 
fused salt. Since two of them may be complexed to the bromide ion, it seems 
probable that only about one molecule of bromine is present in solution in the free 
state. It is possible that complexes are formed containing more than three mole- 
cules of bromine. The degree of association at « = 0.3 is approximately 10 per 
cent. Seemingly, a very slight disturbance of the structure of the liquid electro- 
lyte causes a marked interaction between individual ions with the consequent for- 
mation of ion pairs. 

The high stability of the complex with one molecule of bromine is shown by the 
low vapor pressure of bromine over the fused salts. For a bromine content just 
below that corresponding to one molecule of bromine, the relative vapor densities 
lowerings for Me;N HCl and Me;N HBr were found to be 0.975 and 0.991, respectively. 
For Am;NHCI at « = 0.965, the relative vapor density lowering was found to be 
0.993. The last-named compound is the most stable of the three. 

The marked stability of the higher complexes is illustrated by the relative vapor- 
density lowerings, as shown in Figure 4. The vapor pressure falls off most sharply 
for BusNBr and least sharply for Me;sNHCI; Am;NHCI is next to the tetrabutyl 
compound; evidently the nature of the substituted ammonium ion has an appre- 
ciable effect on the stability of the complex with bromine. 

The broken lines in Figure 4 correspond to the vapor-pressure changes if Raoult’s 
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law held. They are drawn on the assumption that the salt component contains 
only one molecule of bromine. This, of course, is not true, for we know that 
stable compounds of higher order are formed. Nevertheless, a comparison of the 
observed vapor-pressure lowerings with those predicted by Raoult’s law yields 
information concerning the state of the electrolyte in dilute solutions. In all cases 
the vapor-pressure lowering in dilute solutions is much below that required by 
Raoult’s law. This indicates that the molecular weight in dilute solution is greater 
than the formula weight. In short, it indicates that in the more dilute solutions 
the dipoles interact to form quadripoles and perhaps more complex structures. 

In the present investigation, conductance measurements for Me;NHBr were car- 
ried down to a molar concentration of only 0.9. The data therefore do not enable 
us to follow the association process to lower concentrations. However, Moessen 
and Kraus! have measured the conductance of Me;NHCI to 0.036 N, and Rice and 
Kraus? have determined viscosities over the same concentration range. The 
highest concentration to which measurements were carried was 3.0 N, which cor- 
responds to x = 0.33. The value of An at this concentration was 0.976 and was 
leveling off rapidly. For this salt we shall assume a value of 1.1 for (An)o. This 
value may be somewhat high, but it is a sufficiently close approximation for our 


present purpose. 

In Table 6 are given values of A, 7, An, and y, for Me;NHCI at the lowest con- 
centration measured and the concentration of maximum conductance. This con- 
centration corresponds to x = 0.10. In the last line are given values of the ratio 
of the various quantities. At a concentration of 2.01 molar (x = 0.10), the salt is 


TABLE 6 


COMPARISON OF y VALUES FOR Me;NHCI at 
Two CoNcENTRATIONS 
Cc A n X 108 An 100 y 
0.0362 1.217 X 10°? 0.950 1.16 X 10° 0.0106 
2.011 15.75 2.50 0.394 36 
Ratio 55.5 1,216 2.62 3,300 3,300 


dissociated into its ions to the extent of 36 per cent. At 0.0362 N, y is only 0.0106 
per cent. The degree of dissociation has decreased in the ratio of 3,300/1, while 
the concentration has decreased in the ratio of 55.5/1. Obviously the degree of 
association of the ions increases enormously as the concentration of the solution 
decreases. 

There is an equilibrium between ions and ion pairs, and there is a relation be- 
tween the concentrations of the reactants and that of the product. But the rela- 
tion between these quantities is not that postulated by the law of mass action; it 
is just the reverse. It is, of course, the reverse of what occurs in dilute electrolyte 
solutions quite generally. There seems to be no escape from the conclusion that, 
in solutions of these electrolytes in bromine, the degree of ionization decreases with 
decreasing concentration, contrary to what we find in more dilute electrolyte solu- 
tions. 

Without doubt, if the conductance measurements in bromine were carried to 
lower concentrations, a point would be reached where the conductance would pass 
through a minimum and thereafter increase in accordance with the law of mass 
action. All typical electrolyte solutions in solvents of lower dielectric constant 
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exhibit such behavior. It is of interest to note that potassium iodide in iodine 
140° has a minimum at 0.013 N and a maximum at 1.14 N.! 

While data are lacking for more dilute solutions in bromine, we have extensive 
data for such solutions in benzene from | N to very low concentrations, where dis- 
sociation increases with decreasing concentration in accordance with the law of 


mass action. These data overlap those for the more dilute bromine solutions, par- 
ticularly in the case of Me;NHCI. What characterizes these solutions in benzene 
as in bromine is the enormous increase in ionization with increasing concentration. 

The conductance, that is, the ionization, of tri-tsoamylammonium picrate in 
benzene® closely parallels that of Me;NHCI in bromine over the range for which 


we have data for salts in the latter solvent. 
In Table 7 are given values of A for Me;NHCI in bromine and Am;NHPi in ben- 
zene at comparable concentrations. As may be seen from Table 7, the ionization 


TABLE 7 


COMPARISON OF Me;sNHCI 1n BroMINE 
and Am;NHPi IN BENZENE 
A C Ratio A Ratio 
A. Me;NHC1 in Bromine 
0.443 5.58 12: 157 
0.036 0.0122 
B. AmsNHPi in Benzene 
0.35 0.0135 
0.03 1.4% 1-4 


11 

of Am;NHPi in benzene increases 900-fold for an 11.7-fold increase in concentra- 
tion. This is approximately twice that of Me;NHCl in bromine over much the 
same concentration interval. It follows that the phenomena occurring in benzene 
are much the same as those in bromine. We may therefore employ data available 
for benzene solutions to assist in interpreting phenomena occurring in bromine 
solutions. We have extensive data relating to electrolytes in benzene. 

Tri-7soamylammonium picrate is a relatively weak electrolyte in benzene. Its 
dissociation constant® is 4 & 10~-*'. Its conductance passes through a minimum® 
at C = 1 X 10-* N with a conductance of 2.7 K 1077. Its conductance passes 
through a maximum’ at 1.08 N at a value of 0.0475. At 0.01 N the picrate is 
associated to quadripoles to the extent of 6 per cent and to 58 per cent at 0.1 N.* 
In the concentration range 0.01-0.1 N, the conductance, and therefore the ioniza- 
tion, increases 329-fold. In this concentration range, the molecular species are 
ion dipoles with some quadripoles, and the ions are simple ions and triple ions and 
possibly some ions of higher order due to interaction of quadripoles with simple 
ions. 

The conductance in the immediate neighborhood of the minimum may be ac- 
counted for on the basis of triple-ion formation, but the enormous increase in con- 
ductance between 0.01 and 0.1 N cannot be accounted for in this way. Moreover, 
complex ions can be formed only by interaction of simple ions with ion dipoles or 
quadripoles, and the large increase in conductance therefore requires a large in- 
crease in simple ions. It would seem that the mass-action equilibrium, which 
applies at lower concentrations, suffers a reversal in the neighborhood of the con- 
centration of the conductance minimum, as a result of which the dissociation of ion 
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dipoles increases with increasing concentration. It needs to be pointed out that 
this does not apply to the association of dipoles to quadripoles. 

It may prove instructive to consider the behavior of another type of electrolyte 
in benzene. Tetrazsoamylammonium thiocyanate is a relatively strong electro- 
lyte? in benzene (K = 1.93 X 10~").6 Its conductance passes through a mini- 
mum at 2 X 10-° NV, where A = | X 10~*. Beyond the minimum, the conduct- 
ance increases in a series of steps.’ It rises sharply, levels off, and then rises 
again. The plot of log A versus log C shows three rises and two plateaus. At 
1.14 N, the conductance is 0.338. Without doubt, if the concentration could be 
carried somewhat higher, the conductance would pass through a maximum, as do 


other similar salts.’ 

As shown by measurements of molecular polarization,’ the ion dipoles begin to 
form quadripoles at a concentration as low as 5 X 10° N. At 1 X& 107% N the 
association number of the salt, n, the number of formular weights per mole, is al- 
ready 4.5. From here on, 7” increases regularly to a value of 26 at 0.123 molal. 


Thereafter, n decreases sharply to a value of 10.5 at 0.43 molal. 

It is clear that triple ions cannot be present to an appreciable extent at concen- 
trations above 0.001 NV. Neither can ion dipoles be present to any considerable ex- 
tent; otherwise, the values found for the molecular weight of the solute could not 
be so high. The question is: What is the nature of the ions that conduct the 
electric current? It seems that there could be only two answers to this question. 
First, conductance is due to simple ions in equilibrium with charged aggregates. 
Second, conductance is due to charged fragments of the larger neutral aggregates. 

If the first of these hypotheses is true, the ion dipoles in equilibrium with the 
simple ions must be very highly dissociated. It is, of course, possible that only one 
kind of ion dissociates from the aggregates. In that case, we would have simple 
ions of one sign and complex ions of the other. This possibility is not excluded. 
It might serve to account for the stepwise conductance increase. 

If the second hypothesis is true, the solutions must contain a great variety of 
what we may call “fragment ions.’”’ These ions would have to be in equilibrium 
with a great variety of neutral aggregates. One difficulty with this hypothesis is 
that it does not serve to account for the large conductance increases found in the 
case of such salts as Am;NHP1. 

We shall not attempt to answer the questions that we have raised. Rather, our 
purpose is to call attention to a problem which, up to now, has been altogether 
neglected. 

There is one more case that we shall consider briefly. Seward has measured 
the conductance and the viscosity of BuyNPi in butanol at 93° from lower concen- 
trations to the fused salt. 

He found that there is no minimum in the conductance; it decreases continuously 
from the lowest concentration to the fused salt. The conductance-viscosity product, 
however, shows a pronounced and rather sharp minimum at 0.08 NV, and A has the 
same value for the fused salt as it has for the solution at infinite dilution. The evi- 
dence seems to be conclusive that, at concentrations above 0.08 NV, the degree of 
dissociation of the electrolyte increases with increasing concentration. The ion- 
pair equilibrium suffers a reversal at a relatively low concentration. This much is 
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clear; the phenomenon is not complicated by the formation of triple ions, as in other 
cases that we have examined. 


* The investigation here reported was supported by the Office of Naval Research under Con- 
tract N7onr-35809. This paper is based on a portion of a thesis presented by Philip L. Mercier 
in partial fulfilment of the requirements for the degree of Doctor of Philosophy in the Graduate 
School of Brown University, June, 1955. 
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STUDIES OF TYROSINASE PRODUCTION BY A HETEROCARYON 
OF NEUROSPORA* 
By N. H. Horowitz AND MARGUERITE FLING 


KERCKHOFF LABORATORIES OF BIOLOGY, CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA, 
CALIFORNIA 


Communicated by G. W. Beadle, May 31, 1956 


The tyrosinase of Neurospora crassa exists in two forms that differ in their thermo- 
stability.!_ This difference is inherited in a simple Mendelian way, the stable form 


of the enzyme being determined by a gene, 7°, and the labile form by its allele, 7”. 


The locus is on the right arm of the mating-type chromosome, 18 map units (un- 
corrected) distal to adenineless-35203 but proximal to albino-15300. As far as is 
known at present, the only other property in which the two forms of the enzyme 
differ is in the activation energy of thermal inactivation, which appears to be higher 
for the thermolabile form.2,_ Measurements of several functional properties, includ- 
ing Michaelis constants for two substrates and relative turnover numbers, have not 
shown any significant difference between the two enzymes.’ 

Since the T°-T” alleles are recognized by a qualitative rather than a quantitative 
effect on the enzyme, an opportunity was presented for investigating the question 
of whether the alleles act independently in heterocaryons or whether they interact 
so as to produce a hybrid form of the enzyme or so as to show dominance of one 
allele over the other. At the same time it was possible to carry out a test of the 
hypothesis that tyrosinase synthesis is mediated by self-duplicating cytoplasmic 
particles. The principle of this test will be explained in connection with the experi- 
mental results. 

Strains and Methods.—The heterocaryon was produced by fusion of hyphae from 
T* and 7” strains carrying, respectively, the mutant genes adenineless-35203 and 
histidineless-C140. The strains had previously been made nearly isogenic by six 
successive backcrosses of ad-7* progeny from the cross ad-T* X hist-T” to the hist- 
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T” parent. The nutritional markers served to stabilize the heterocaryon on minimal 
medium and to mark the nuclei in the second part of the experiment. Both homo- 
caryotic components of the heterocaryon were of mating type A. The constitution 
of the heterocaryon was verified by showing that both original homocaryons could 
be recovered from a single hyphal tip.* Conditions of growth and of extraction and 
purification of the enzyme were as described previously.!. Enzyme activities were 
determined colorimetrically with pt-DOPA as substrate.! Thermostability tests 
were carried out at 59° + 0.02° C. and pH 6. 

Tyrosinase Production by the Heterocaryon.— Measurements of the thermostability 
of tyrosinase produced by the heterocaryon and by the two parent homocaryons are 
shown in Figure 1. It is seen that thermal inactivation of the enzyme from the 





WwW 
a 
> 
— 
> 
— 
oO 
<« 


hist-TL 











1 1 i 1 1 1 1 1 1 1 1 1 
10 20 30 40 50 60 70 80 90 100 110 120 
MINUTES AT 59° 

Fic. 1.—Thermal inactivation of tyrosinase from two homocaryotic 
strains and the heterocaryon between them. The half-life of the 
enzyme from ad-7’S is 70 minutes in this experiment; the half-life of the 
enzyme from hist-7” is 3.4 minutes. In the case of the heterocaryon, 
the curve drawn through the experimental points is the theoretical for 
a mixture containing 52.3 per cent of tyrosinase of half-life 70 minutes 
and 47.7 per cent of tyrosinase of half-life 3.4 minutes. 


homocaryotic strains is a first-order process, the half-lives of the two forms at 59° 
being 70 and 3.4 minutes, respectively, in agreement with previous determinations. ! 
Inactivation of tyrosinase from the heterocaryon does not follow the first-order 
course expected in the case of dominance of one allele or of interaction to form an 
enzyme of intermediate stability. Rather, it closely follows the curve predicted on 
the basis of a mixture of the two original enzymes. The relative proportions of the 
two activities are deducible from the first-order tail of the curve, which represents 
the inactivation of the stable component of the mixture. By extrapolating this 
portion of the curve back to zero time, the initial amount of thermostable activity 
‘an be read off the y-axis; the amount of thermolabile activity and the expected 
course of the thermal inactivation are then easily derived. In this way it can be 
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shown that the curve of Figure | is that predicted from a mixture containing 52.3 
per cent thermostable and 47.7 per cent thermolabile tyrosinase activity. 

Resolution of the Heterocaryon.—To test the theory that production of the enzyme 
is mediated by self-duplicating cytoplasmic particles, the heterocaryon was re- 
solved into its genetic components by cutting off and subculturing hyphal tips; a 
fraction of such tips can be expected to be homocaryotie by chance. Since there are 
no effective crosswalls in the hyphae of Neurospora, the genetically unlike nuclei of a 
heterocaryon share a common cytoplasm. On the hypothesis of self-perpetuating 
cytoplasmic agents, homocaryotic cultures 
derived in the manner described might 
therefore be expected to resemble the 
parent heterocaryon in their tyrosinase 
production, rather than the original homo- 
caryotic strains. 

Of 1,220 hyphal tips taken, 638 survived 
on transfer to minimal medium supple- 
mented with adenine and histidine. Of 
these, 6 proved to be homocaryotic, as 
shown by the fact that they failed to 
grow on minimal medium, Transfer to 
medium supplemented with adenine or 
histidine showed that 3 of the isolates were 
adenineless (and therefore 7°) and 3 were 
histidineless (therefore 7”). These isolates 
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were next grown in large cultures, and 
their tyrosinase was extracted, purified, 
and tested. The results, plotted in Figure 
2, show that all the homocaryons produced 


> 


only one form of tyrosinase, and this in 
accordance with the genetic constitution 








a a a of the nucleus. The presence of | per cent 
a oicabaaee nd Pra oe of thermostable activity in the thermolabile 
Re oe Rea ey preparations would have been detected, or 
sinase from six homocaryons obtained by the presence of 10 per cent of labile activ- 
pera erlaen, oan ae Ty/Mat-T, ity in the stable preparations. “The half 
fit half-lives of 70 and 3.1 minutes, respec- lives determined for the six enzyme prepa- 
tively. rations were 3, 3.4, and 3.5 minutes for the 
enzyme from the histidineless cultures and 
55, 69, and 74 minutes for the enzyme from the adenineless cultures. The stable 
form is generally more variable than the labile form in thermostability tests, owing 
to the fact that it is influenced to a greater extent by impurities in the prepa- 
rations.' The effective impurity seems to be principally tyrosine, traces of which 
are apparently released by a proteolytic enzyme which contaminates some of the 
tyrosinase preparations. Like other tyrosinases, Neurospora tyrosinase is inacti- 
vated in reacting with its substrates, and the presence of tyrosine in the prepa- 
rations causes an apparent decrease in the thermostability of the stable form; the 
labile form is hardly affected, since it. is thermally inactivated too rapidly at 59° to 
engage in many catalytic cycles. 
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Discussion.—The absence of allelic interaction in the synthesis of tyrosinase by 
the heterocaryon parallels findings which have been made in the genetically similar 
case of sickle-cell anemia in man. It has been found that sickle-cell heterozygotes 
produce both the normal and the abnormal form of hemoglobin.‘ The same thing 
is true for the ;- and £.-lactoglobulins of cattle, which, on the basis of gene-fre- 
quency analysis of the data of Aschaffenburg and Drewry,° appear to be determined 
by a pair of allelic genes. It remains to be seen whether this will turn out to be 
a general rule for genes governing the structure of proteins, as would be expected 
on the template theory of gene action. 

The failure to detect cytoplasmic transmission of tyrosinase-forming capacity in 
hyphal tip transfers is consistent with the known Mendelian inheritance of the 
thermostability difference. In order to account for the rarity of known cases of 
cytoplasmic inheritance compared with chromosomal heredity, it has been suggested 
that plasmagenes may be subject. to decay in the germ line. The present experi- 
ments lend no support to this idea, in so far as it draws a distinction between the 
germ line and the soma. Decay of the hypothetical plasmagenes in our case must be 
presumed to take place in asexual as well as sexual generations. It may be pointed 
out that the present experiments are not subject to the interpretation that persist- 
ence of the hypothetical plasmagenes depends on the presence of particular chromo- 
somal genes (as in the case of kappa in Paramecium’), since both T-alleles support 
tyrosinase synthesis. 

In view of the fact that the RNA (ribonucleic acid) of tobacco mosaic virus has 
recently been shown to possess genetic properties,® the continued failure to detect 
genetic activity in the cytoplasm—which contains RNA—in anything approaching 
the variety of activities associated with the chromosomes poses an interesting prob- 
lem. One possible answer is that two structurally different forms of RNA exist, 
one capable of self-duplication and the other not.® 

Summary.—aA heterocaryon between two strains carrying, respectively, the alleles 
T* and 7”, governing the thermostability of tyrosinase, was found to produce both 
forms of the enzyme. These alleles thus appear to act independently in the synthesis 
of the enzyme. The heterocaryon was resolved by means of hyphal tip transfers, 
without going through the sexual stage. No evidence was found for cytoplasmic 
transmission of tyrosinase-forming ability. 


We wish to acknowledge the invaluable assistance of Mrs. Anne Kilim Nazarian. 
Dr. R. R. Jennings participated in the mapping of the T-locus. We are indebted 
to Dr. Dan Lewis, of the John Innes Horticultural Institution, for stimulating 
discussions. 


* Aided by grants from the Nutrition Foundation, Ine., and the Atomic Energy Commission 
(AT-04-3-41). 
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THE GLUCOSE CONTENT OF THE DEOXYRIBONUCLEIC ACTDS 
OF CERTAIN BACTERIOPHAGES* 


By Roperr L. SINSHEIMER 
PHYSICS DEPARTMENT, IOWA STATE COLLEGE 
Communicated by M. Delbriick, April 19, 1956 


The presence of glucose in the deoxyribonucleic acid (DNA) of the even-numbered 
bacteriophages of the T series was discovered by Jesaitis' and independently by 
Sinsheimer.2. This glucose was shown to be associated with the hydroxymethyl- 
cytosine (HMC) in these nucleic acids as a glycoside substituent to the hydroxy- 
methyl group (Sinsheimer;? Volkin®*). 

Further studies have indicated that the proportion of glucose is not the same in 
T2L and T4 bacteriophages and that the proportion of glucose seems to be connected 
with other recognizable traits of the bacteriophages (Streisinger and Weigle‘). 
However, the mode of inheritance of these properties is unusual. 

Methods.—All phage stocks were grown on Escherichia coli ‘‘B”’ according to the 
usual methods (Adams°) in the glycerol-casamino acid medium described by Fraser 
and Jerrel.6 The stocks were purified by treatment with pancreatic deoxyribonu- 
clease, followed by two cycles of differential centrifugation. DNA was obtained 
from the purified viruses by osmotic shock (Anderson et al.”), followed by high- 
speed centrifugation to remove ‘‘ghosts,”’ deproteinization according to the method 
of Sevag et al.,° and alcoholic precipitation. 

The molar proportion of glucose was determined by the ratio of the glucose 
content, as determined by the use of the anthrone reagent for glucose (Morris®), 
to the phosphorus content (Allen'®). A correction for the color produced by the 
deoxyribose in the anthrone test was made by subtraction of the color produced 
by a calf thymus DNA solution (Gary and Klausmeier'') of the same phosphorus 
content. Recovery experiments with added glucose gave quantitative results. 

Enzymatic Degradation of Bacteriophage DN A.— lf DNA from any of these bacterio- 
phage strains is degraded successively with pancreatic deoxyribonuclease and puri- 
fied venom phosphodiesterase, about 65 per cent of the phosphorus of the digest 
can be recovered as mononucleotides.?, The remainder is present as di-, tri-, and 
larger polynucleotides, all of which appear to contain HMC, as indicated by their 
ultraviolet absorption spectra. 

If the DNA is that of a wild-type T2L, five mononucleotides are obtained: 
thymidylic acid, deoxyadenylic acid, deoxyguanylic acid, and two nucleotides of 
HMC, one with and one without a glucose substituent. The molar proportions of 
these mononucleotides are presented in Table 1, together with the proportion of 
rach nucleotide in the T2L DNA which is obtained as a mononucleotide. — It is clear 


TABLE 1 
MONONUCLEOTIDES FROM Enzymatic Dicest or T2L DNA 
HMC- 
Thymidylic Adenylie Guanylie HMC Glucose 
Per cent of P of digest in 
mononucleotide 22.8 22:2 12.5 2:2 3.0 
Per cent of nucleotide in 
digest recovered = as 
mononucleotide 70 68 69 é 24 


502 
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that the proportion of the glucose-substituted HMC mononucleotide found among 
the 65 per cent mononucleotide fraction is much lower than that of the other 
nucleotides, and thus, conversely, the enzyme-resistant polynucleotides must be 
enriched in this component. It seems likely that the presence of the glucose sub- 
stituent, appropriately situated, confers resistance to phosphodiesterase action 
upon the polynucleotide. 

The presence of the HMC nucleotide lacking glucose in such digests is clear 
evidence for incomplete glucose substitution of the HMC. 

Analysis. The nucleic acids of T2L and T4 bacteriophage and of the progeny 
of various crosses of these viruses have been analyzed for glucose content and 
have been enzymatically degraded to ascertain the possible presence of hydroxy- 
methyleytosine lacking glucose substitution. The results are presented in Table 2. 

TABLE 2 
Presence of Presence of 
HMC Nucleotide, — Glucose- 


Moles Glucose Not Glucose- substituted 
Virus Strain 100 Moles P substituted HMC Nucleotide 


T2L rthtu* 12.6 + 
T4r‘*u 16 _ 
BX9 T2 r*h*ut 17 
BX9 T2 r*htu 16: 
V-G T2 rhu* 12 
V-G T2 rh*u* 16 


+ 


CO wl | 


T4 differs from T2, among other traits, by the presence of a locus u conferring to 
its carrier a high resistance to inactivation by ultraviolet light (Luria;!* Strei- 
singer). The strains B X 9 T2 are strains of T2 which, after one mating with 
T4, have been backcrossed to T2 nine times to make them as isogenic as possible 
with T2.* '% It can be seen that they have acquired the “glucose property” of 
T4 and that it is not the u locus which is responsible for the glucose content. The 
V-G T2 stocks resulted from a single burst of a cross of T2rh+ with B X 9 T2rh 
with delayed input of the second phage (Visconti and Garen'*). These experiments 
are described in more detail in the next paper. Under the conditions of delayed 
input, it can be seen that not all the progeny particles of the cross acquire the glu- 
cose property of T4. 

Discussion.—In the wild-type T2L bacteriophage some (23 per cent) of the 
hydroxymethyleytosine is not substituted with glucose, while in T4 all the hydroxy- 
methyleytosine is glucose-substituted. The inheritance of this ‘‘glucose property” 
is unusual in that it does not seem to be linked with any of the markers r, h, or u 
but does seem to be linked with the “bar’’ markers, which, in a cross T2 & T4 
with simultaneous infection, are transmitted to the vast majority of the progeny.‘ 

It is possible that the bipartite nature of the DNA of T2, suggested by Brown 
and Martin and by Levinthal,'® may be responsible for this unusual inheritance 
pattern. 


Cohen (Science, 123, 653-656 [1956]) has recently presented data suggesting a 
relationship between the glucose content of bacteriophage DNA and the r property. 
The data presented here do not support this relationship. 


* The research described in this paper was supported by grants from the American Cancer 
Society, acting through the Committee on Growth of the National Research Council, and from 
the United States Public Health Service. 
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PROPERTIES OF BACTERIOPHAGES T2 AND T4 WITH 
UNUSUAL INHERITANCE* 
By G. STREISINGERT AND J. WEIGLE 
‘ipieesbiat 5 cuenmenete eat tniiiiaainn 
Communicated by M. Delbriick, April 19, 1956 


The first experiments on infecting bacteria mixedly with two different phages 
(Delbriick and Luria;' Delbriick?) were done with pairs of unrelated phages. 
They led to the discovery of exclusion, i.e., the finding that any one bacterial cell 
will yield mature particles of only one or the other of the parental types, never both. 
The nature of the exclusion mechanism is not understood. 

Later it was found that exclusion is a function of the degree of relatedness of the 
two infecting phages. Phages which differ by only a few mutational steps are 
perfectly compatible (see Hershey*); not only do both phages multiply, but they 
also undergo extensive genetic recombination. Phages somewhat more distantly 
related show an intermediate behavior (Delbriick and Bailey‘), the details of which 
have never been worked out. 

The phages T2 and T4 belong to a very large group of related phages (Adams’). 
They are the phages whose genetic maps are better known (Hershey and Rotman ;* 
Doermann and Hill’) than those of any other phage. They are also the pair for 
which genetic recombination and partial exclusion‘ were first observed. In the 
present paper we will describe the partial exclusion of T2 by T4 and will show that 
the “excludability” of T2 is itself a genetic trait of T2 which is completely excluded 
from the progeny. The excluding powey of T4 is thus inherited in a novel way: 
it is imparted to the overwhelming majority of the progeny particles, irrespective 
of what other markers of either of the parents these particles inherit. Moreover, 
this excluding power of T4 is associated with two other characteristics. One of 
these is the characteristic that the particles carrying it, in contrast to those which 
do not carry it, plate with high efficiency on certain bacterial indicator strains. 
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The character is thus easy to score. The scoring of the other characteristic is 
laborious, but the character is of utmost theoretical interest: the particles with the 
excluding power of T4 have a different deoxyribonucleic acid (DNA) from those 
not carrying this power. In their DNA, 100 per cent of the hydroxymethyleytosine 
(HMC) is substituted with glucose, while, in the particles with the exclusion proper- 
ties of T2, only 77 per cent of the HMC is substituted with glucose (Sinsheimer’). 


MATERIALS AND METHODS 


The phages used in these experiments were’ T2 L and some of its mutant de- 
rivatives carrying the markers'’-'? h and re, and T4. The bacterial strains were 
Escherichia coli B and its mutant'? B/2, as well as derivatives!* '4 of E. coli K-12 
called 112-12 and C-600. The general phage techniques are those described by 
Adams." 


EXPERIMENTAL 


The bacteriophage properties having an unusual inheritance pattern are as 
follows: 

1. Efficiency of Plating (e.0.p.) on Strain 112-12.—T4 when piated on 112-12 
gives as many plaques as when plated on B, although the plaques on 112-12 are 
somewhat smaller. Plaques produced by r and r+ phages on 112-12 are indis- 
tinguishable. T2 plated on young cultures of 112-12 gives only about 6 per cent as 
many plaques as on B; the plaques which do appear are so small as to be almost 
invisible. The T2 plaques which do appear on 112-12 contain phage identical to 
the parental T2 with respect to the low e.o.p. on 112-12. Strain C-600 exhibits 
similar relations to T2 and T4. Plaques formed on B by T4, upon replication with 
velvet (Lederberg and Lederberg;'® Baylor’) on plates seeded with 112-12 and 
preincubated for an hour or on plates seeded with C-600 and preincubated for 
three hours, give plaques, while T2 plaques replicated under the same conditions 
give no plaques when replicated on either strain 112-12 or strain C-600. 

2. The Unequal-Yield Effect.—In a cross of T2 X T4 on B bacteria with equal 
multiplicities, the markers introduced by T4 occur among the progeny particles 
with a frequency three to ten times higher than those of T2 (Streisinger''). 

3. The Glucose Content of DNA.—In the DNA of T4, 100 per cent of the base 
5-hydroxymethyleytosine (HMC) is glucose-substituted, while in the DNA of T2 
only 77 per cent of the HMC is glucose-substituted. The studies on the definition 
of DNA are reported by Sinsheimer in an accompanying paper.* 

In across between T2 and T4 most, if not all, of the progeny acquires the efficiency 
of plating of T4 (Table 1). T2 (T2 bar) designates T2 phages which have acquired 
the e.o.p. of T4 through a cross with T4 or with another T2. It is shown below 
that the phages which have acquired through a cross with simultaneous infection 
the e.o.p. of T4 on strain 112-12 have also acquired the other two properties. 
Collectively the three properties will be referred to as the “bar properties.”’ 

In the experiments to be described on the transmission of the bar properties, we 
utilized as donor of these properties, instead of T4, a T2 strain which had been 
made more or less isogenic with T2. Such a strain had previously" been obtained, 
for the study of the u locus, by proceeding broadly as follows: First, T2 was made 
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isogenic with T2 for most loci, except the uw locus of T4 (controlling the level of 
ultraviolet sensitivity) through a series of eight backcrosses.'5 The « locus was 
replaced by the u* allele (from T2) in the last backcross. The progeny of the last 
(ninth) backcross will be symbolized as Bx 9 T2. 


TABLE 1 
SINGLE Bursts or A Cross or T2 X T4* 
No. Plaques No. Plaques 
Not Not 
No. T4 No. T2 Replicating Burst No. T4 No. T2 Replicating 
Plaques Plaques on C-600 No. Piaques Plaques on C-600 


l 30 0 6 24 7 0 
11 38 0 7 106 133 0 
36 28 0 8 18 32 0 
42 5 0 9 65 55 0 


a) 
5 26 1] 0 10 43 59 0 


* The cross was performed on bacteria of strain B grown to 108/ml, centrifuged, rg nded in buffer, and 
aerated for 90 minutes (to assure simultaneous infection; see S. Benzer, J. Bacteriol., 63, 59-72, 1952). The in- 
fected bacteria were distributed into individual tubes, so that the mean number of infec ted. B per tube was 0.26. 
The contents of the individual tubes were plated on B + B/2 after incubation for several hours at 37° C. On these 
plates T2 particles form turbid plaques and T4 particles form clear plaques 


In a cross of T2 XK Bx 9 T2 the inheritance of the bar properties could thus be 
studied within an almost identical genetic background. The properties of T2, 
Bx 9 T2, and T4 are contrasted in Table 2. It can be seen that the three bar 
properties introduced only once in Bx 9 T2, nine crosses earlier, had persisted and 
had remained together. The properties studied are not closely linked to the u, A, 
or r loci, since, no matter what marker is carried by the progeny of the last back- 
cross, the phages are identical with respect to the three properties (Table 2). 


TABLE 2 
ErricieNcy oF PLaTiInG, UNEQUAL-YIELD Errect, AND PER Cent orf GLUCOSE 
SuBstTiruT1IoN CoMPARED IN T2, Bx 9 T2, anp T4 Srrains* 
Ratio of Test 
Efficiency of Marker to T2 
Plating on Marker in Crosses Per Cent HMC 
Strain 112-12 with Approximately Substituted 
Strain (Per Cent) Equal Input with Glucose*® 
T2rthtut 12 1 77 
Bx 9 T2rth*u* 100 6T 100 
3.6T 
Bx 9 T2rth*u 100 6t 100 
Bx 9 T2rhu* 100 10.6T Not tested 
T4r*u 100 3.0t 100 
* The efficiency of plating indicates the ratio of the number of plaques on strain 112-12 to the number of plaques 


on strain B. The plaques of T2 on strain 112-12 are barely visible. 
The ratio of the number of markers of the strain studied to the number of markers of T2, among the progeny of a 


cross between the strain tested and T2, defines the unequal-yield effect. For instance, for Bx 9 T2r*th*tu? the 


cross (Bx 9 T2r*th*tu*) X (T2rhu) was performed with equal input multiplicity, and the ratio of h* to h in the 
yield measured. 


+ The multiplicity of T2 to T2 in these experiments was equal, but the total multiplicity varied. It was noted 


in later experiments that the total multiplicity in such crosses had an effect on the ratio of T2 markers to T2 markers 
in the yield. Some of the variability in this series of numbers may be due to this factor. 


In the crosses described thus far, the three bar properties seem to be imparted to 
every particle in the progeny. It seemed important to find out whether the bar 
properties were transmitted together under conditions in which not every progeny 
particle is converted into the bar form. These conditions can be realized by giving 
the non-bar parent a head start of several minutes. By delaying the input of the 
T2 parent until 8 minutes after the initiation of development of T2 (Visconti and 
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Garen'’), single bursts are obtained which contain a large number of non-bar T2 
particles and very few T2 particles which are like T4 with respect to e.o.p. on strain 
112-12. A typical burst of a cross T2rh+ X T2rh is described in Table 3. Two of 
its particles were analyzed, one h* replicating on 112-12 and one h not replicating on 
112-12. The first one which was thus bar with respect to e.o.p. proved to be bar 
for glucose (Sinsheimer*) and for the exclusion properties. The other which was 
non-bar for e.o.p. was non-bar for glucose, non-bar for excludability, but was bar 
for excluding ability. 


TABLE 3 
YIELD OF A SINGLE BactertuM INFECTED WITH T2rh* AND SECONDARILY WITH 
T2rh, Yretp1nG Boru T2 anp T2 Parricies* 

No. T2h* No. T2h 
No. h* No. h Replicating Replicating 
Particles Particles on 112-12 on 112-12 

136 3 1 0 
* The cross was performed on bacteria of strain B grown to a concentration of 2 X 10’/ml, centrifuged, resus- 


pended in buffer at a concentration of 5 X 10%, and aerated for 90 minutes. The bacteria were infected with an 
average multiplicity of two T2rh* particles per bacterium, and, after 4 minutes for adsorption, were diluted by a 


factor of 2 into concentrated nutrient broth and aerated at 37° C. After 8 minutes of development, T2rh particles 
were added at an average multiplicity of 400 per bacterium. Serum was added after 4 minutes to inactivate 
unadsorbed phage, and, after 2 minutes for serum action, the infected cells were distributed into individual tubes. 
After incubation for several hours at 37° C. to permit lysis to occur, the contents of the individual tubes were 


plated on B + B/2 to score the number of A particles. The number of T2 particles was scored by replication onto 
plates seeded with strain 112-12 and preincubated for 1 hour. There were 0.15 bursts per plate. 


The total exclusion of the non-bar properties and the partial exclusion of other 
markers from the progeny of a T2 X T2 cross with simultaneous infection may be 
related. In order to understand how these total and partial exclusions take place 
and to find out what connection exists between them, an experiment on the unequal 
vield was performed and is described in Table 4. 


TABLE 4 
Tue Unequat-Yie_p Errect in A Cross or T2rh X Bx 9 T2rtht* 
Fraction of 
Infective : 
oe As Ratio Centers Effective 
MULTIPLICITY OF of h/h* Yielding Multiplicity Ratio 1n YIELD OF 
T2rh Bx 9 T2rtht in Input Noh of h h/h* r/r* 
3.4 4.8 0.71 0.33 1.1 0.15 0.15 
* The cross was performed on bacteria of strain B grown to a concentration of 10’/ml, concentrated to 108/ml, 
and suspended in broth containing 10~* M potassium cyanide (to assure simultaneous infection; see S. Benzer and 
F. Jacob, Ann. Inst. Pasteur, 84, 186-204, 1953). After unadsorbed phage particles had been inactivated by anti- 
serum, infective centers were plated before lysis on B + B/2. The rest of the culture was diluted and incubated for 
2 hours at 37° C. The progeny thus obtained was plated on B + B/2. : i 
The effective multiplicity m of the T2 parent was calculated from the fraction of infective centers, z, which 
liberated no particles with an h marker: z = e™. 


The unequal yield of markers carried by the two parents of a cross of T2 X T2 
could be due to the complete elimination of a fraction of the T2 parental particles 
or to a relative decrease, by each particle, of its contribution of markers to the 
progeny (due to slower growth, for instance), or both effects might occur. 

The effective multiplicity m of the T2 parent in a cross of T2rh X T2rth+ can be 
calculated from the fraction x of infected bacteria which yield progeny not having 
the h marker of the T2 parent, according to the Poisson equation: x = e~”. It 
can be seen from Table 4 that a large fraction of the h markers of the T2 parent are 
completely excluded in mixed infections with T2, since the effective multiplicity of 
the T2 parent is less than one-quarter of the input multiplicity. The ratios of h to 
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h*+ and r to r* markers in the progeny of this cross are identical (0.15), and similar 
to, though slightly less than, the ratio (0.23) of the effective T2 multiplicity (1.1) to 
the T2 multiplicity (4.8), indicating that the unequal yield in this cross is due 
principally to exclusion. 

DISCUSSION 

In simultaneous infections of B cells with T2 and T4 the markers of T2 are 
partly excluded; in addition, the excluding power of T4 is transmitted to most, if 
not all, of the progeny. Thus the “excludability” of T2 is itself completely ex- 
cluded. Strains resulting from a cross of T4 X T2, and subsequently backcrossed 
to T2 nine times, so as to be virtually isogenic with T2, still retain the excluding 
power of T4. 

Two other properties seem associated with the excluding power of T4: a high 
efficiency of plating on certain strains of bacteria as compared to a low efficiency of 
T2, and a difference in the DNA of the phages. The HMC of T4 is 100 per cent 
glucosylated, while that of T2 is 77 per cent glucosylated. 

The backcrossed T4 derivatives which are more or less isogenic with T2 exhibit 
the first two properties of T4; the few strains tested contained 100 per cent gluco- 
sylated HMC as well. These properties are called “bar.” 

When infection of bar phage is delayed until several minutes after infection with 
non-bar phage, mixed bursts are obtained, containing bar and non-bar phage, in 
contrast to simultaneous infection. A particle which is bar with respect to one of 
the properties (e.0.p.) proved to be bar with respect to the other properties as well. 
A particle with the non-bar e.o.p. proved to be non-bar with respect to glucose, bar 
with respect to excluding ability, and non-bar with respect to excludability. These 
findings open the way to a study of the relations between the bar properties. 

The partial exclusion of markers carried by the non-bar parent and the total 
exclusion of the non-bar property itself, in crosses with simultaneous infection, can 
be understood with the aid of a more or less conventional genetic model, as well as 
with other models. For instance, if the non-bar properties were determined by loci 
which are eliminated if—and only if—there is present with them in the bacterium a 
bar locus introduced by another phage, only particles carrying the bar loci would be 
observed among the progeny. Markers carried by non-bar particles could, how- 
ever, be incorporated into vegetative particles with the bar loci by means of recom- 
bination and thus be ‘‘rescued”’ in a manner analogous to the rescue of markers of an 
ultraviolet inactivated particle in crosses with unirradiated phage. Loci close to 
the eliminated ones would, in this case, be expected to be rescued more rarely than 
distant loci. The study of the rescue of differently linked markers could yield 
information concerning this model. Such a model allows also an understanding of 
the separability of the bar properties. 

It is tempting to try to relate exclusion of the non-bar properties to the different 
glucose content of non-bar and bar DNA. Cohen? has shown that DNA not 
containing HMC, and thus not containing glucose, is broken down more completely 
by certain enzymes than is DNA containing glucosylated HMC, and Sinsheimer* 
has shown that glucose-containing HMC is less readily separated from polynucleo- 
tides than non-glucosylated HMC. Thus non-bar DNA, containing less glucosy- 
lated HMC than bar DNA, could more readily be subject to a breakdown mecha- 


nism. 
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The gross chemical and biological differences found between T2 and T2 phage 
may be due to structures other than classical genetic loci. A structural differ- 
entiation of T2 DNA has been suggested by Levinthal,?! who finds one large and 
several smaller pieces of T2 DNA, at least half of the larger piece being conserved 
throughout replication. A fractionation of the DNA of T2 into two fractions 
differing in their base ratio has been obtained by Brown and Martin”? after a rela- 
tively gentle removal of DNA from the protein coat. 

It is likely that a structurally bipartite T2 DNA is also functionally bipartite, 
consisting of material with the usual markers and of material with other functions. 
The difference in the glucose content of bar and non-bar DNA may be found to 
reside exclusively in the fraction of DNA not carrying markers. In that case the 
difference between bar and non-bar particles would be due not to loci but to a 
structurally separate piece of DNA. Construction of a detailed model involving 
such a separate piece at this time seems premature. 

The material described here promises to be useful for an understanding of the 
mechanism of partial exclusion and thus for an understanding of the establishment 
of the metabolic regime leading to the production of the DNA of a particular phage, 
in competition with that of another phage and of the bacterium. 

SUMMARY 

1. Among the differences observed between the phages T2 and T4 are the follow- 
ing: 

a) Markers of T2 are partly excluded in crosses of T2 & T4. 

b) T4 plates with a high efficiency, while T2 plates with a low efficiency on 
certain strains of bacteria. 

c) One hundred per cent of the 5-hydroxymethyleytosine (HMC) of T4 is 
glucose-substituted, while only 77 per cent of the HMC of T2 is glucose-substituted. 

2. These three properties, unlike the usual genetic markers, are transmitted to a 
vast majority, if not all, of the progeny of crosses of T2 K T4. These properties 
are referred to collectively as the “bar’’ properties, and T2 particles having these 
properties are symbolized as T2. 

3. The progeny of a cross of T2 X T4, made subsequently more or less isogenic 
with T2 by repeated backcrosses, retain all three of the bar properties. 

4. If progeny from single bacteria are infected with T2 and several minutes 
later are infected with T2, the following observations are made: 

a) The progeny contain particles with the T2 efficiency of plating as well as 
particles with the T2 efficiency of plating. 

b) Some of the particles seem to show a separation of the bar properties. 


The authors wish to express their indebtedness to Dr. M. Delbriick for his fruitful 
criticisms and for his unfailing aid in the redaction of this paper. 


* This investigation was supported by a grant from the National Foundation for Infantile 
Paralysis. 
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THE PRODUCTION OF TWO CHEMICALLY DIFFERENT TYPES OF 
CHROMOSOMAL BREAKS BY IONIZING RADIATIONS* 


By SHELDON WoLrr AND Henry FE. Luiproup 
BIOLOGY DIVISION, OAK RIDGE NATIONAL LABORATORY, OAK RIDGE, TENNESSEE 
Communicated by Karl Sax, May 31, 1956 


The production of chromosome aberrations by ionizing radiations has been 
studied extensively, mainly because it is a particularly fruitful source of information 
on the quantitative relations between radiation dose and cellular damage. Pri- 
marily, radiation produces chromosomal breaks in direct proportion to the dose ad- 
ministered.'! Ninety-five per cent of these breaks, as calculated by Lea? for Trades- 
cantia, restitute, i.e., rejoin in their original configuration. The others may stay 
open as one-hit aberrations or rejoin with still other broken ends to form two-hit 
aberrations. These aberrations can be formed only if the two breaks are open 
simultaneously and in close proximity. If the radiation is administered at low 
intensity, thus allowing some of the breaks to restitute before others are produced, 
then fewer two-hit aberrations result than if the same dose were administered at 
higher intensity.* 

Intensity studies of this nature in 7’radescantia led Catcheside, Lea, and Thoday* 
to propose that two types of breaks are produced by the radiation: one group that 
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rejoins relatively rapidly and another that does so more slowly. Studies on the 
combined effects of near infrared and X-rays led Kaufmann and Wilson’ to a similar 
proposal for Drosophila chromosomes. 

Previous intensity experiments on the seed of Vicia faba indicated that X-ray- 
induced chromosome breaks stay open for at least two hours when the seed has been 
soaked in water.6 However, breaks produced either in BAL’ or in vacuo® remain 
open for only one-half hour. Wolff and Luippold*’ have postulated that since 
this type of break stays open for long periods of time and since respiration and ATP 
are necessary for its rejoining, it is the break of a strong, possibly covalent, bond. 
However, Mazia,'’ Steffensen,'! and Levine! have reported that chromosomes are 
held together by the ionic bonds formed by the divalent cations of caleitum and 
magnesium, and Steffensen'! suggested that radiation breaks these bonds. Since 
breaks of ionic bonds cannot be expected to stay open for long periods of time, it is 
extremely unlikely that any of this type were observed in previous experiments 
with Vicia. The present 
experiments were therefore 
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Seeds of V. faba were . = = to to 
soaked in distilled water for oes 
24 hours and then peeled Fic. 1.—Two-hit aberrations produced by 600 r at differ- 
and arranged in a culture ent intensities. 
dish with their embryos 
facing upward. The dishes were exposed to 600 r of X-rays. For irradiation 
under anoxic conditions, the seeds were transferred to a 2 X 10~-* M solution of 
BAL one-half hour before exposure. In those final experiments that combined 
treatments with Versene and radiation, a different crop of seeds was used. Their 
response to radiation was the same as the crop used in the previous experiments if 
the seed coats were chipped and the seeds were soaked for only five hours. All ir- 
radiations were performed with a G.E. Maxitron tube operated at 250 kvp. with 
3 mm. of Al filtration added (hvl., 0.4 mm. of Cu). The seeds were germinated 
as previously described'* and fixed in C. E. Ford’s modification of Flemming fixa- 
tive (unpublished). Feulgen squashes were made, and at each point 300-500 


DURATION OF IRRADIATION (min) 


metaphases of the first mitotic root-tip division were scored for rings and dicentrics. 

The intensity curves in Figure | are characterized by plateaus that represent the 
intervals within which variation of the exposure time has no effect on aberration 
yield. These plateaus can exist only for times within which no appreciable re- 
joining occurs, for, otherwise, fewer open breaks would coexist at low intensities 
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than at high and the number of two-hit aberrations would be correspondingly 
altered. The length of the plateau wherein there is no diminution of two-hit 
aberrations is therefore a measure of the time the breaks stay open. This time, 
which is dose-dependent,’ is at least two hours in seeds irradiated in water and one- 
half hour in those irradiated in BAL. If the duration of treatment exceeds these 
times, rejoining occurs and the typical intensity effect results. However, at radia- 
tion intensities higher than 200 r per minute, the number of aberrations produced 
in air (water) and under anoxia (BAL) increases above the plateau level. This 
indicates the presence of another type of break that rejoins (either restitutes or 
forms a two-hit aberration) very rapidly. If the radiation is administered at a very 
high intensity, many short-lived breaks are present in the cell simultaneously and 
are capable of yielding many two-hit aberrations. In seeds irradiated with the 
same dose, but at lower intensity, these breaks usually restitute, and only those 
that do not rejoin for a relatively long period of time remain. The latter are the 
breaks, described in previous reports from this laboratory, that are responsible for 
the plateau on the intensity curves, and are probably of a covalent nature. Thus 
two sharply distinct classes of breaks can be recognized, one with a lifetime of only 
about one minute and another that persists from one-half to two hours. 

The very rapid rejoining of breaks observed only at the high intensities is consist- 
ent with the hypothesis that they are breaks of ionic bonds. For a further test of 
this hypothesis, seeds were soaked for varying times in a 0.001 M solution of ethyl- 
enediamine tetra-acetic acid (Versene). This concentration was chosen because 
it was effective in the studies of Mazia' on Drosophila and Melanoplus that led him 
to postulate that the chromosome consists of small units bound together by calcium 
and/or magnesium ions. Table I shows that low concentrations of Versene can 
break the chromosomes of Vicia, too. This can be interpreted as supporting evi- 
dence, in plant material, for Mazia’s and Steffensen’s models of the chromosome. 

TABLE | 
PERCENTAGE OF ABERRATIONS PRODUCED BY 
SOAKING SEED IN VERSENE 
(0.001 M) 


Time in No. 
Total Versene Aberrant Aberrations 
Cells (Hours) Cells (Per Cent) 


650 0 6 0.94 + 0.38 

750 3 27 3.60 + 0.69 

250 6 36 14.40 + 2.4 

We then reasoned that, since a three-hour soaking in Versene produces some 

chromosome breakage, this period must be long enough to insure penetration of the 
chelating agent into the cell. If subsequent irradiation at low intensities then pro- 
duced “‘ionic”’ breaks in the chromosomes at the locations of calcium and/or mag- 
nesium ions, chelation would occur and the breaks, instead of restituting very rapidly, 
would stay open to become available for the production of aberrations. The results 
of a typical experiment of this type are reported in Table 2, which shows that pre- 
irradiation treatment with Versene in this instance gave a twofold increase in total 
aberration yield. It should be mentioned that the results of this type of experi- 
ment are unusually variable, both between slides and between experiments. Thus, 
although the increased aberration yield is qualitatively certain, the quantitative re- 
sults may be poorly reproducible. 
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Steffensen’s earlier observations!! that the amount of chromosome breakage pro- 
duced by combined calcium deficiency and radiation is more than the additive 
amount of both is consistent with these results. The same rationale used to in- 
terpret the Versene-radiation effect can also be used to interpret the caleium de- 
ficiency—radiation effect. 


TABLE 2 


PERCENTAGE ABERRATIONS PRODUCED BY COMBINED VERSENE 
(0.001 M7) AND RapiatTion TREAT- 
MENTS * 
Percentage Percentage 
Slide Total Normal Percentage Dicentrics 
No. Treatment Cells Cells Fragments and Rings 
50 96.0 2.0 2.0 
Versene, 3 hours , #0 100.0 0.0 0 
3 : . : 50 98.0 0 0 
Total) 50 98.0 33 67 
50 66.0 0 >.0 
50 0 0 0 
50 0 0 
600 r 50 r 0 0 
44 9.1 3.6 
6 56 r.6 85 
Total 300 74.6 3.6 3.3 
l 50 4- 0 32.0 
2 50 46 0 é 0 
3 Versene 3 hours 00 : 0 34.0 
: slus 600 r 50 ° 
5 prs 50 ’ 0 0 
6 50 0 0 
Total | 300 : 0 


¥ 


* Seeds soaked for a total of 5 hours. 


It is possible, therefore, to interpret our results on the basis of two chemically 
different types of chromosome break produced by radiation. One type is a break 
of covalent bonds that require energy to close. The evidence accumulated here 
about the other type, although not conclusive, is consistent with the hypothesis 
that it is the break of an ionic bond. 

It is interesting to note that the shapes of the Vicia intensity curves, after incor- 
poration of the high-intensity points, are also reflected in the curves obtained when 
Drosophila sperm are irradiated. Haas et al.'* have reported an increased number 
of translocations from irradiation at very high intensities. These translocations 
are analogous to the aberrations formed from “‘ionic’’ breaks in Vicia. These 
breaks disappear rapidly from the system by restitution, which is the reason they 
were never observed in previous experiments at relatively low intensities. There 
is no further intensity effect, because the other breaks, presumably of covalent bonds, 
stay open until the time of fertilization. These, then, are the equivalent of the 
breaks that give rise to the plateau of the Vicia curves. Since sperm respire very 
little, the breaks cannot rejoin until they are supplied with the requisite energy pro- 
duced by respiration of the fertilized egg. 

Thus in Drosophila sperm it may be postulated that the same two kinds of chro- 
mosome break exist and that the length of the plateau is determined not by the re- 
covery of a “rejoining system” as in Vicia* but by the availability of the energy 
produced in the egg. 
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Summary.—Experiments on the seed of V. faba indicate that radiation-induced 
chromosome breaks may have two chemically different natures. One seems to be 
an ionic break, which closes very rapidly, as would be expected if only electrical 
factors were necessary for its rejoining. The other is covalent, and it stays open 
for long periods of time and needs a source of energy for the biosynthesis of the bonds 
formed in its closing. Such apparently unrelated experimental results as those ob- 
tained by the use of calcium deficiency,'" '? chelating agents,” !? metabolic in- 
hibitors,* ° and high-intensity radiation’ can be united by this concept of two 
types of chromosomal damage. 

* Work performed under United States Atomic Energy Commission Contract No. W-7405- 
eng-26. 
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FREQUENCY OF SOMATIC MUTATION TO SELF COLOR IN MAIZE 
PLANTS HOMOZYGOUS AND HETEROZYGOUS 
FOR VARIEGATED PERICARP* 


By DonaLp R. Woopt anp R. ALEXANDER BRINK 


COLORADO AGRICULTURAL AND MECHANICAL COLLEGE AND UNIVERSITY OF WISCONSIN 


Communicated June 28, 1956 


Introduction.—Maize plants heterozygous for variegated pericarp and cob (P?"’) 
and a stable allele (e.g., P”’, colorless pericarp and red cob) bear ears which are 
more heavily striped than those from homozygous variegateds. The basis of the 
difference is the number of mutations of P’’ to P” (self colored pericarp and cob) 
in pericarp tissue. Emerson’ ? who first established these facts, estimated that the 
single P*’ allele in P’’P”” plants mutated to P’”’ about 2.8 times as frequently as 
each of the P’’ alleles in variegated homozygotes. Tests which we and our colleagues 
have made with P*’ alleles of various geographic origins and on several genetic 
backgrounds have revealed that the numerical relationship differs in different back- 
grounds and, rarely, may even be reversed. Frequently, however, the difference 
is comparable to that which Emerson observed. 

It will be shown in the present study that the difference in number of mutations 





Vou, 42, 1956 GENETICS: WOOD AND BRINK 515 


to self color in PP” and PP"? plants can be accounted for satisfactorily in terms 
of dosage action of Modulator (7p), a transposable element associated with the 
variegated allele. 

The results of previous investigations* are consistent with the view that the 
variegated allele is a dual structure, comprising the P’” gene and Modulator. It 
may be represented, therefore, by the extended formula P’Mp. Mp, in this posi- 
tion, inhibits the pigment-producing action of P’. The mutation of variegated 
to self color consists in removal of ./p from the P locus, ordinarily by transposition 
to another site in the genome. It has been demonstrated also that one or more 
such transposed Modulators in a variegated pericarp plant markedly reduces, in an 
orderly way, the frequency of P’’ to P’’ mutations.* 4 

This latter fact suggests the reason why mutations to self color are less frequent 
per P’’ allele in P’’/P"’ than-in P’’) P’’ plants. Each of the Modulators in the 
former genotype may be expected to reduce the frequency of P’’ to P’” mutations 
in the homologous P’° allele, whereas the single P’’ allele in heterozygotes is not 
subject to this modifying action. 

If the mutation rate of a P°" allele is equally influenced by an additional Mp in 
any position, including the locus of the homologous allele, then twice as many 
mutations to self color regularly will occur in homozygous variegateds as in hetero- 
zygous variegateds carrying a transposed Mp in hemizygous condition. This 
result would be expected because there are two P’ alleles which can undergo muta- 
tions to P’”’ in the former genotype as compared with only one in the latter, and one 
Modulator is present in each case, not counting the one at the mutating allele. 
The experimental data obtained in the present study verify this prediction. 

Materials and Methods.—The maize strains used were variegated counterparts 
of a commercial yellow dent inbred line known as W23 (P"’). The initial outeross 
to variegated pericarp had been followed by five backcrosses of the P°’/ P”’ segre- 
gates to W23, resulting in a high degree of uniformity in genetic background. Three 
genotypes then were selected from this inbred stock, differing in the expression of 


variegated pericarp as follows: (1) P?’.Mp/P*’, designated “medium variegated” ; 
(2) P’Mp/P’’ + tr Mp, —, designated “light variegated”; (3) P’Mp/P”Mp, 
designated “homozygous variegated.’’ Genotypes 1 and 2 are alike in possessing 
a single variegated allele, but genotype 2 carries a hemizygous transposed Modulator 
also. Genotypes 2 and 3% differ in having one and two variegated alleles, respec- 
tively, and each has two doses of Mp. Genotypes | and 3 are heterozygous and 
homozygous for the variegated allele, neither possessing a transposed Modulator. 


The ears borne by these three classes of plants are illustrated in Figure 1. It is 
seen that the three phenotypes are readily distinguishable, and, of course, they 
differ in the character of their offspring in testcrosses with P””. 

The light variegated phenotype is represented in Table 2 by two families. The 
transposed Modulators carried in these families resulted from independent mutations 
of P’’. The significance of the fact that the frequency of mutation of P*’® to P”’ in 
the two cases is closely comparable will be discussed later. 

Medium variegated, light variegated, and homozygous variegated ears from the 
converted W253 inbred stock were scored for mutations to self color. Size of the 
mutant area was recorded, and the results were summarized in terms of number of 
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red stripes or patches per 1,000 kernels. Size of stripe was determined with the aid 
of models. The latter were constructed in accordance with the findings of Ander- 
son and Brink® that size is a function both of time when the stripe was initiated and 
position on the kernel. The models were obtained by superimposing a circle, with 
its center on the point of silk scar and a diameter parallel to the germinal face, upon 
scale drawings of a kernel as viewed from the top. The circle was then divided into 
eight equal sectors. The length of the included are of sectors falling on the ab- 
germinal face of the kernel was greater than those falling on the germinal face be- 
cause the kernel develops asymmetrically in relation to the point of silk attachment. 
Stripes covering one segment were called one-eighth-kernel mutations and were the 
smallest scored. Stripes traversing two sectors were scored as mutations covering 
one-fourth of the kernel. Additional size classes of mutants recorded were those 
covering one-half kernel, one kernel, and two or more kernels. 
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Fic. 1.-Three genotypes of variegated peri- 
carp in maize. Left to right: medium varie- 
gated (P/P’"), homozygous variegated 
(P’/P), and light variegated (P°/P’" + 
tr Mp/-—). 


The mutations were distributed into five size classes, as the stripes did not always 
fall on the sectors described. These classes were as follows: one-eighth to less 
than one-fourth, one-fourth to less than one-half, one-half to less than one, one 
kernel, and two or more kernels. The classes are referred to by their minimal val- 
ues rather than by their mid-points because the minimal values are in geometric 
progression with the one-kernel class, which has no range. Mutant stripes covering 
less than one-eighth of a kernel, although frequent in medium variegateds particu- 
larly, are difficult to score and so were not recorded. 

Experimental Results. The distributions of offspring in the families from which 


representative ears were scored for somatic mutations to self color are shown 
in Table 1. The phenotypic grouping of the variegateds is somewhat arbitrary, 
since the “light”’ ears in one class of families, although similar, did not correspond 
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precisely with those in another. 
relevant to the present experiment. 
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“Very lights” are mutants of a special type* not 
These data demonstrate the homogeneous 


character of the stocks used, the regular distribution of the ?°’ allele, and the trans- 


mission of transposed Modulator independently of the P locus. 


GENOTYPE 
VARIEGATED 
PARENT 


Pr Per 


Pr /Per + tr Mp/— 
Prev / pr 


TABLE 1 


ACCORDING TO CLASS OF VARIEGATION 


Very 


PEDIGREE 
. Light 


6A-40-8 @ 


6-291-1 & W238 
Bk 293-1 @ 2 


7) 


] 


COLORED OFFSPRING 
Self 
Light Medium Colored 
20 
19 62 20 
23 20 2 
24 0 2 


DISTRIBUTION OF OFFSPRING OF DIFFERENT GENOTYPES 


COLORLESS 
PERICARP, 
Rep Con 


54 
43 
0 


TOTAL 


25* 
156t 
88 
28 


* Pov/ Pee offspring. 
¢ Per/Pwr and Pur/Pwr offspring. 


The number of mutations per 1,000 kernels for medium variegateds, light varie- 
gateds, and homozygous variegateds and their distribution according to stripe 


size are shown in Table 2. Differences between families within a given variegated 


TABLE 2 
DISTRIBUTION OF NUMBER OF MutatTIoNns oF P°® to P™ 
PER 1,000 KERNELS 


Est1- 
MATED No. No. 
> 


Toran Ears 2or 
KERNELS Scored More ] 
4,857 10 0.62 8.65 1 
4,351 92 9.42 1 
I 
l 


MvTATIons PER 1,000 KERNELS 
0.5 0. 25- 
0.125 
39.12 
32.64 
36.06 
22.06 
16.54 
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( 
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) 24.5 
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3 20. 3: 
12 
16 


pre 


Pwr 
9,208 76 Ol ( 
9,295 86 10 
9,736 at G.af 7¢:8 
19,031 58 6.73 10.7 
5,466 28 48 6.95 
10,344 10 2.51 t 
15,810 0.13 85 5 


86.88 


Pre Pre 


53.76 
26.45 


genotype appear to be due to fluctuation around a common mean and so may be 
combined in the analysis. The values plotted in Figure 2 are based upon the total 
mutations in each variegated genotype. 

The frequency of somatic mutation to self color in medium variegateds (P"’/ P”’) 
proved to be about 1.6 times as high as that in homozygous variegateds in these 
stocks. Although the observed difference is of the same order as that originally 
reported by Emerson? for these two genotypes, the magnitude, as such, of the dif- 
ference is not important. As pointed out in the introductory paragraph, the amount 
of the difference may be expected to vary, and does in fact vary, from one strain 
to another, depending both upon the particular P’’ allele concerned and the geno- 
typic background. On the other hand, with a given P’’ allele on a given genetic 
background, there should be a constant relation (2:1) between the mutation 
rates in homozygous variegateds and light variegateds if Modulator dosage is the 
controlling factor. 

As the data in Table 2 and Figure 2 show, mutations of P’’ to P’’, in fact, are 
almost exactly twice as frequent in homozygous variegateds (P°"/P*") as in light 
variegateds (P""/P”’ + tr Mp/—). The chi-square test shows a close fit on this 
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basis for each of the kernel-stripe classes which can be compared. Further con- 
firmation of the relationship is the fact that the ratio of the regression coefficients 
of the mutation frequencies for the homozygous variegated and light variegateds 
was 1.87, a value which is not significantly different from 2.0, according to the (test. 
Discussion.—A premise implicit in the design of this experiment was that trans- 
position of Modulator from the P locus to a position elsewhere in the genome is 
not necessarily accompanied by a change in state (mutation) of Modulator. The 
validity of this assumption is borne out by the fact that the mutation rates of P’" 
to P” in the two light variegated families scored, 6-291-1 and 6A-40-8, carrying 
Modulators independently transposed from the P locus to unknown but almost 
certainly different sites, were essentially the same. Inspection of the ears based 
on two additional light variegated muta- 
40, tions within the W238 inbred stock 

pYYpwr ; 
showed that they also were alike and 
indistinguishable in phenotype from the 
particular light variegateds used in the 
present experiment. This evidence sup- 
ports the view that the modifying effect 
of Modulator on P’’ to P’” mutations is, 
in fact, independent of the location of 

Mp in the genome. 

It would be expected, therefore, that 
Mpasa component of a P"’ allele would 
have the same action on mutation of a 
second P*° allele in the nucleus as a 
transposed Modulator. The data pre- 
sented in Table 2 fully substantiate this 


n 
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z 
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° 
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° 
= 
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expectation. The total number of P°° 
to P’” mutations of all sizes exceeding 








one-eighth kernel on the P’’/P”° ears is 

53.76 per 1,000 kernels. The corre- 

iting sti ip : ye as sponding value for the light variegateds 
r1G. 2.——Number of mutations of PF" to Pr > , r * Of AF ; ‘ 
per 1,000 kernels, distributed according to size (P(r Mp —) is 26.45, or almost 


of the mutant area. The abscissa is on a geo- exactly one-half. Both genotypes have 
metric scale, which gives an approximately 
linear relation. 


NUMBER OF KERNELS IN MUTANT AREA 


one Modulator in the genome, in addi- 
tion to the one directly associated with 
the mutating allele. In other words, a Modulator at one allele acts as a trans- 
posed Mp with respect to the other. The validity of this relationship is borne out 
by the further fact that, within sampling limits, the distribution of the mutations 
between the different stripe-size classes in the two arrays also conforms to expecta- 
tion on this basis. 

It is interesting to recall at this point the properties of the system which Rhoades® 
described controlling mutations of a; to A; in the presence of Dotted (Dé) in maize. 
Rhoades found that the number of such mutations was influenced by dosage of 
both a, and Dt. The frequency of mutation rose linearly with the dosage of a, but 
at an accelerating rate as the number of Dé factors was varied from one to 
three. 


, 
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Thus both Dt and Mp are nonlinear in their respective actions, but the dosage 
effects are of opposite direction. 

Summary.—The effect of Modulator on stability of the variegated pericarp 
allele (P°° or P’-Mp) in somatic tissue was studied in three otherwise near-isogenic 
genotypes designated “medium variegated” (P’Mp/P"’), “light variegated”’ 
(P"Mp/P”’ + tr Mp/—), and “homozygous variegated” (P’Mp/P”Mp). As 
Emerson had earlier observed in unrelated stocks, the single variegated allele in 
plants heterozygous for a stable P gene was found to mutate to self color approxi- 
mately three times as frequently as either of the alleles in variegated homozygotes. 
It was postulated that the lower rate in the homozygote was the result of the action 
of the Modulator component in each variegated allele on mutation to self color of 
the other variegated allele. The test of this hypothesis which was applied was a 
comparison of the frequency of mutation in homozygous variegateds, carrying two 
P*’ alleles embodying two Modulator units, and light variegateds, carrying one 
variegated allele and one transposed Modulator. If the mutation rate of a P* 
allele is equally influenced by an additional 7p in any position, including the locus 
of the homologous allele, then twice as many mutations to self color would be ex- 
pected in the homozygotes as in the light variegateds. The numerical results ob- 
tained on scoring the two classes of ears in question for mutations to the self color 
agreed closely with expectation on this hypothesis. 

* Paper No. 626 from the Department of Genetics, College of Agriculture, University of Wis- 
consin. This study was aided by grants from the research committee of the graduate school, 
of funds supplied by the Wisconsin Alumni Research Foundation. 

+ This paper is based on a portion of a thesis presented by Donald R. Wood in partial fulfilment 
of the requirements for the degree of Doctor of Philosophy in the Graduate School of the University 
of Wisconsin. <A sabbatical leave from Colorado Agricultural and Mechanical College for the 
year 1953-1954 is gratefully acknowledged. 
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ON MAUTNER’S EIGENFUNCTION EXPANSION 
By WiiuraM G. BADE AND Jacos T. SCHWARTZ 
UNIVERSITY OF CALIFORNIA AND YALE UNIVERSITY 


Communicated by N. Jacobson, May 29, 1956 


The present note had its origin in a paper of Mautner! on the general eigenfune- 


tion expansion problem for self-adjoint operators on a Lebesgue space and recent 
papers of Garding? and Browder* which apply Mautner’s results to elliptic differ- 
ential operators. 

Let (S, 32, v) be a o-finite measure space and 7' be a (possibly unbounded) self- 
adjoint operator in Lo(S, >>, v). A version of the spectral theorem asserts that 
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there exist vectors }fat such that Lo(S, 52, v) is the orthogonal direct sum >> Ha, 
where Ha is the manifold of all vectors F(7’)fa, where F € Lo(ua), wa = (1+) fa, fa). 
Moreover, the correspondence F'(7')fa — F(+) establishes an isometric isomorphism 
U of H onto the direct sum >>L.(ue) which diagonalizes 7’ in the sense that 
(UF (T)f)a(X) = F(A) (Uf )a(d) for fe D(F(T)). Following Mautner, we say that 7’ 
has an eigenfunction expansion if for any such diagonalizing decomposition of Lo(S, 
a v) there exist kernels Wa(s, \) such that the corresponding operator U’ is deter- 


mined by the equations 
(Uf)a(A) = S5 f(s) Wels, d)v(ds) (1) 


the integral existing in the mean-square sense‘ in /o(ue), and the inversion formula 


f(s) = SS2. (Uf)a(d) Wals, X)pa(dr) (2) 


holds, the integrals existing in the mean-square sense‘ in Lo(S, >, v) and the 
series converging in the topology of Le(S, >>, v). 

1. The Main Theorem.—Our main theorem establishes a sufficient condition 
that an operator 7 should possess an eigenfunction expansion. In later sections 
we discuss the relation of this result to those of Mautner and Garding. More com- 
plete details will appear in a forthcoming book.’ It should be noted that we as- 
sume that (S, >>, v) is o-finite but do not assume that L.(S, >>, v) is separable. 

1.1 Derrnirion: A bounded linear operator A in Lo(S, >, v) has| property 
(a) if there exists a sequence |S, of sets of finite measure covering S such that 

y-ess sup | (Af) (s)| < @, el ee 
Sn 
for each f € Lo(S, 2 v). The sequence Saf of sets will be called a covering sequence 
for the operator A. 

If A has property (a), the map Af — (Af!|S,,) of (Af) (-) into its restriction to 
S, is a closed and hence bounded linear map of Le(S, a v) into L.(S,, a v), 
forn = 1,2,.... In fact, this is an equivalent formulation of property (a). 

1.2 Tueorem. Let T be a self-adjoint operator in L.(S, >o, v) with the resolution 
of the identity E (+). If for each bounded Borel set e of the real line the projection E(e) 
has property (a), then T has an eigenfunction expansion. 

1.3. Lemma. Under the hypothesis of Theorem 1.2, there exists a common cover- 


ing sequence }S,| for all the projections E(e), i.e., a sequence |S, such that 


v-ess sup | (E(e)f) (s)| < ©, os (ee 
Sn 
for every bounded Borel set e and f € LoAS, ad v). 

Clearly, if | C,,} is a suitable covering sequence for E(e) and e € e, it is a suitable 
covering sequence for E(e,). Now let @ be a finite positive measure on (S, >>) 
such that @(e9) = O implies v(e)) = 0. Let e, = [—n, n], and select sets A, ¢ Ssuch 
that 


v-ess sup | (E(e,)f) (s)| < @, fe LAS, d, v), 
An 


while the 6-measure of the complement A,’ of A, satisfies 0(A,’) < 2-“+», If 





Vou, 42, 1956 MATHEMATICS: BADE AND SCHWARTZ 


Si = f? A,, then 0(8S,;’) < '/2 and v-ess sup | (L(e)f) (s)| < © for every bounded 
n=1 Si 

Borel set e and f € L2(S, : v). The construction of So, Ss, ... proceeds by indue- 

tion. 

We now construct the kernels Wa(s, A). 

1.4. Lemma. Let |S,{ be a covering sequence for S such that formula (3) holds for 
every bounded Borel set e and f ¢€ LAS, >>, v). Let g € LAS, pa v), and let u(-) = 
(E(+)g, g). Then there exists a function W defined on the Cartesian product of S and 
the real line P, measurable with respect to v X wu, such that, if e is any bounded Borel 
subset of R, 

(i) v-ess sup Jf. | W(s, d)|2 u(dd) < @, 
seSn 
(ii) (KH(e)F(T)g) (s) = S Wis, \)F()u(dr), Fe Le(u). 

Let e, = [—n,n]. It follows easily that the map A,: F(-) > F(T)g is a closed 
and hence bounded map of Lo(e,, «) into L.(S,, 52, v). Consequently A,* is a 
bounded map of L..*(S,, >., v) into Le(e,, wu). Now L,(8,, 42, v) is isometrically 
imbedded in L.*(S,, , we v), and the restriction B, of A,* is a bounded map of 
L\(S,, 3., v) into Le(e,, w). Since Lo(e,, u) is separable, it follows from a theorem 
of Dunford and Pettis® that there exists a (v X u)-essentially unique » X u-meas- 
urable complex-valued function W“” defined on S, X e, such that 


(iii) (Bf) A) = fs, f(s) W™ (s, \)o(ds), f € LS, >, »), rN €@n, 


and such that 
v-ess sup ff, | W™ (s, d)| %u(dd) = ||B,) 2 < &. 
Sn 
The function W is defined by the formula 
W(s, A) = W™(s, dA), (s, A) eS, X ey. 


Formula (ii) ‘s easily established, using Fubini’s theorem. Q.E.D. 

Now let fa be a family of vectors determining a diagonalizing decomposition of 
LS, >, v), and associate with each fa the corresponding kernel Wa determined 
by Lemma 1.4. It follows by relatively standard arguments (cf., e.g., Browder, 
op. cit.) that for each f ¢ Lo(S, }°, v) the integrals J’s, f(s) Wa(s, \)v(ds) and 
S a” (Uf)a(X) Was, )ua(dd) exist for almost all \ and s, respectively. Moreover, 


(Uf)a(A) = lim Ss. f(s)Wals, A) v(ds), f € LoS, > 3 v), (4) 


n—>o 
and 
( af) (s) = lim La (Uf)a(A) Wa(s, A) ual(dr), f € LAS, ti v), (5) 
n—> co 
where Pa is the perpendicular projection of L.(S, >>, v) onto Ha, the limits being 
in the topology of Lo(ue) and Lo(S, >>, v), respectively. Formula (2) follows from 


the fact that f = OP f. 


1.5. Coroutiary. A self-adjoint operator T in Lo(S, >>, v) has an eigenfunction 
expansion if there exists a sequence |S,,| of sets of finite v-measure covering S such that 


every fin f¥ D(T™) is v-essentially bounded on each of the sets Sy. 
m=1 





522 MATHEMATICS: BADE AND SCHWARTZ Proc, N. ALS. 


This useful corollary follows from the observation that (e)L(S, 30, v) ¢ D(T™) 
for every m when e is bounded. 

Mautner! showed that if, for some A, the resolvent (A; 7’) of a self-adjoint oper- 
ator 7 in L2(S8, >>, v) was an integral operator of Carleman type (i.e., 

(RA; T)f) (s) = Ss K(s, Of(Ov(dt), 
where K is product-measurable and K(s, +) € Le(S, >, v) for almost all s), then 7 
has an eigenfunction expansion. Since the range of R(A; 7) is D(7T), the domain 
of 7, it follows immediately from Mautner’s hypothesis that for each g = R(\; T)f 
e D(T) we have 
g(s)| < |f|0(s), 


where 6(s) = [fs K(s, t)|*v(dt)]*. On taking S, = 1s! 6(s) < nt, Mautner’s 
result follows from Corollary 1.5. Garding? generalized Mautner’s theorem for 
separable L2(S, >, v) to the case that some function F(7’) of T is a Carleman oper- 
ator, where F(\) does not vanish on o(7’).7. In the case that inf | F(A)! > 0, 


[—n, n] 
n 1,2,..., it is easily seen that Garding’s hypothesis implies that each of the 
projections E({—n, n]) has property (a). (Otherwise we obtain hypothesis (@) for 
all the projections E(e), where inf | F(A)! > 0, the set [A F(A) = 0} playing the 
=e 
role of the point A = © in Theorem 1.2. <A straightforward modification of 
Theorem 1.2 will then yield the more general result of Garding. ) 

2. Elliptic Operators.—Let S be an open subset of real Euclidean n-space, and 
let + be a formally symmetric elliptic differential operator with sufficiently dif- 
ferentiable coefficients. Then 7 determines a symmetric operator 7o(7) defined 
on the C® functions whose supports are interior to D. If 7 is any self-adjoint ex- 
tension of 7'o(7), Garding and Browder prove that an appropriate power of the 
resolvent of 7 is a Carleman operator and hence that 7’ has an eigenfunction ex- 
pansion. The functions Wa(+, A) belong to C” and (r — A)W,(+, A) = 0. Their 
construction of the Carleman kernel is based on a result of Fritz John which estab- 
lishes the existence of a local fundamental solution for 7. We shall now indicate 
how the theorems of Garding and Browder may be proved simply via our Corol- 
lary 1.5, without the introduction of fundamental solutions. We shall say a func- 
tion f e Lo(S, do, v) possesses the partial derivative 0f/Oz' in the mean-square sense if 


Be Pa iy 6 VE owe gy Be) PE oe Rh cate) = of 
h—0 : Ou' 
exists in the topology of Lo(S, >°, v). A function f belongs to the class H, if it 
has mean-square partial derivatives of all orders <p in each interior subdomain 
of S. If p > [n/2] + 1, then, by a result of Sobolev,’ each f in H, is continuous. 
It is known’ that if m is the order of 7, then D(7(r)*), the domain of the adjoint 
of 7o(7), consists of all functions f in H,, such that rf e Lo(S, >>, v), and To(r)*f = 
tf. Moreover, if (7'(7)*f, h) = (g, A) for all h € D(To(r)) and if g « H,, then f « 
Hysm. Thus, if T is any self-adjoint extension of 7 (7), the functions in D(7’) 
(¢ D((To(r)*)’) are continuous if mr > [n/2] + 1. It follows from Corollary 1.5 
that 7 has an eigenfunction expansion, for we may take for |S, any sequence of 
compact sets whose union is S. 
In a recent paper L. Hérmander” has shown that a certain class of nonelliptic 
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partial differential operators 7 with constant coefficients, called by him operators 


of local type, also have the property that each f « f D((T(r))") is continuous. 
n=1 

It follows then, by Corollary 1.5, that any self-adjoint extension of T(r) has an 

eigenfunction expansion; this was shown by Hérmander, using Garding’s method 

of proof. 

3. Carleman Operators. —Our next theorem shows that when Lo(S, >>, v) is sep- 
arable, the hypotheses for Theorem 1.2 ere actually equivalent to the requirement 
of Garding that an appropriate function of 7 should be a Carleman operator. 

3.1. THrorem. Let T be a self-adjoint operator in a separable space L(S, >, v), 
and let each projection E(e), bounded, have property (a). Then there exists a bounded 
function F defined on the real line such that inf |F(\)| > 0, n = 1, 2,..., and 


[—n, n] 
F(T) ts an integral operator of Carleman type. 

This result follows from a simple characterization of Carleman operators given 
in the next lemma. Since it may have independent interest, we state it for a pos- 
sibly nonseparable space. 

3.2 Lemma. Let (S, >o, v) be a o-finite measure space. A bounded operator in 
L(S, 2, v) is of Carleman type if and only if it has a separable range and has property 
(a). 

Proof: Suppose 


(Af) (s) = Ss K(s, Of(t)v(dt), fe LAS, = v), 


where K is product-measurable and 6(s) = fs | K(s, t)|?v(dt) < © ae. We 
have already noted such an operator has property (a). To see that A has a sep- 
arable range, construct operators A, with kernels K,, such that K,,(s, t) | < K(s, t)| , 
lim K,,(s, t) = K(s, l) v X v-a.e., and SS K,,(s, t)| 2v(ds)v(dt) << ©. Then each A, 


7 
is a compact operator and hence has a separable range. For EF e >, define (P xf) (s) 
= f(s)Xx(s). If sup 0(s) < ©, we note |Px(A,f — Af)| > 0, f ¢ LS, LY, »). 
E 

If E, = }s\0(s) < n}, then |Af — Py, Af| > 0, f « Lo(S, &, v). From this it 
follows easily that A has a separable range. 

A proof of the sufficiency may be patterned on the proof of Theorem 1.2. Since 
A has a separable range, we may find a sub o-field > of }> such that AL(S, }, v) 
¢ L(S, >>, ) ¢ LS, ¥, v) and Li(S, &, v) and L.(S, 5, v) are separable. The 
map C,: f — (Af|S,) of LA(S, >, v) into L.(S,, >, v) is closed and 
bounded. The restriction B, of C,,* to Li(S,, te v) is a bounded map into L.(S, 
>, v) having a separable range. Hence, by the theorem of Dunford and Pettis 
cited above, there exists a v X v-essentially unique measurable function K™ de- 
fined on S, * Ssuch that 


(B,f) () = Ss, f(s)K (s, Ov(ds), 


f e LA(S,, eae v), teS, 
and such that 


v-ess sup fs | K"(s, t)| 2v(dt) = ||B,||? < @. 


seéeSn 
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The proof may be completed as before, to show that 
(Af) (s) = fo K(s, Of(Ov(dt), = fe LS, &, v), 
where K(s, t) = K’” (s, t), (s, t) eS, XS. QED. 
To prove Theorem 3.1, let each projection /(e), e bounded, have property (@), and 
select a sequence {S,|, with the properties of Lemma 1.3, covering S. There exist 
numbers 7, such that 


v-ess sup | (2 ([—n, n])f) (s)| < VM, )f), s €S,. 


Define F(A) = > Mi 2> "xen, nl (A), Xn, n) denoting the characteristic func- 
n=1 


tion of the interval [—n,n]. Then 


> (E({—n, n])f) (s) 


Ad a oe n=1 
(F(T)f) ( OM, 


f € LAS, > Vv), 


the series converging uniformly on each of the sets S,. It follows that F(7) has 
property (a), and thus /(7) is a Carleman operator by Lemma 3.2. 

4. Linear Independence of the Kernels Wa.—If Lo(S, 32, v) is separable, it is 
known that there exists a diagonalizing sequence }f,{ for 7’, a finite Borel measure 
uw on the rea) line R, and a sequence R = e; D3 e D ... of Borel sets such that 
(E(e)fn, tn) = ple N e,),n = 1,2,.... Weshall call the corresponding decomposi- 
tion an ordered decomposition. Let W;, 7 = 1, 2, ..., denote the corresponding 
kernels whose existence is proved in Lemma 1.4. Then W,(s, \) may be supposed 
to vanish for \ ¢e;. Since W; is v X u-measurable, the function W,(+, A) is v- 
measurable for u-almost all A € e;. 

4.1. THrorem. The v-measurable functions W,(+, Xd) W,(+, A) are linearly 
independent for u-almost all d in en. 

The proof proceeds by induction. Suppose that the functions W,(-, A), 7 = 1, 
..., p — 1, are linearly independent for g-almost all A on e,-1, while the set 
Wi(-,A),¢ = 1,..., p, is not linearly independent for u-almost all A on e,. By a 
fairly complicated argument, it is shown there exists a Borel subset ¢ ¢ e, of posi- 
tive u-measure and Borel-measurable functions a; such that 

pl 
W,(*, A) = ‘> a (A) Wis, A), rN € 49, 
t=1 
the principal tool being a vector form of the well-known theorem of Lusin. By 
integrating this equation over appropriate sets (cf. Lemma 1.4 [ii]), we may contra- 
dict the orthogonality of the manifolds H,,..., H,. Complete details will be pub- 
lished elsewhere. 

We may define the multiplicity of a projection E(e) in the resolution of the iden- 
tity for 7’ as the number of sets e; such that u(e A e;) + O in an ordered decom- 
position which diagonalizes 7. The multiplicity of the identity will also be called 
the multiplicity of T. Let T be a self-adjoint extension of the symmetric operator 
7To(7) corresponding to an ordinary differential operator 7 of order n on an interval 
I. Theorem 4.1 may be used to show that the multiplicity of 7 is <n. The 
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kernels W,(-, A) belong to C”(/) and satisfy the boundary conditions at a fixed end 
point of J or, more generally, at an end point at which there is a purely discrete 
spectrum (in a suitable sense). As a consequence we obtain additional information 
as to the magnitude of the multiplicity of 7’: e.g., let + be a real second-order oper- 
ator on an interval. If the boundary conditions determining T include a boundary 
condition at a limit-circle end point, then 7 has multiplicity one." 
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ON THE DISTRIBUTION OF CERTAIN FUNCTIONALS 
OF MARKOFF CHAINS AND PROCESSES 
By D. A. DaruinG ANp A. J. F. Stegert* 
UNIVERSITY OF CHICAGO AND NORTHWESTERN UNIVERSITY 
Communicated by G. E. Uhlenbeck, June 7, 1956 
1. Funetionals of Markoff Chains.—Let x,, n = 0,1, ... be a Markoff chain 
taking on values in an abstract space 2 and having stationary transition probabilities 
P,(a, BE) = Priani,e Bla, = x}. Let Pola, E) = x(x), the characteristic function 
of the set FL, and suppose x» = x, a constant. Let V(y), y « Q, be real valued and 
measurable with respect to ¥, the Borel field of sets over which the P,’s are proba- 
bility measures. In this section we study the distribution of the functional 
U, = > V(z,). (1.1) 
j=0 


This distribution will be determined, in principle, if we know the function 


R(z, E) = > Efe*’"\a, = 2, 2, € E}P,(x, E)z", 0 < z <1, & real, (1.2) 


n=0 
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and Theorem | below gives R as the unique solution to Fredholm type equations. 
Set 


P,(z,E) = >> P,(z, Ez", 0 < 2< 1. 
n=0 
Then we have 
THeorEM 1. The function R of equation (1.2) satisfies uniquely either of the two 
integral equations 


R(x, E) PAa, BE) + SU — e7*®*™) Pa; dy) R(y, BE), (1.4) 
R(x, E) = PAz, E) + f A -e7 VO) R(x, dy) Py, E). (1.5) 


The proof of this theorem is accomplished by integrating the trivial identity 
n n ; a 
eS 2, P — 7 ai > {}i— eh (si) 6%, | (x,) 
P 0 


over the set x, € &, which has probability P,(z, £). Then, on multiplying by 2” 
and summing on n = 0, 1, ..., we obtain equation (1.4). Equation (1.5) follows 
similarly by starting with the identity 


e* E vias) = 1+ > (1 — eV yg, E Veep, 
j=0 

It is simple to show that if Ry satisfies equation (1.4) and Ry, equation (1.5), 
then necessarily Ry, = Rg,, and since the function R of equation (1.2) satisfies 
both of them, the uniqueness assertion of the theorem follows. Equations (1.4) 
and (1.5) might be termed the “forward” and “backward” integral equations for R. 

We notice that for fixed x, &, z, the function R of (1.2) is, as a funetion of F, a 
bounded, complex-valued, completely additive measure over §. 

We obtain somewhat more directly the distribution of U, by putting # = Q in 
equation (1.2). Thus, defining 


Q(x) a; : E} e* "| x, = z}2", 


n=0 


we obtain from equation (1.4) the following equation satisfied by Q: 


I ? EV(y) 
Q(x) = : + f (l-e VO) P (yz, dy) Q(y). (1.6) 


Ifz, =r+t+ytyet ... + Yn, 2% = x, where the y, are independent, identically 
distributed random variables with common characteristic function g(t) = E}exp 
(ity,)}, it will follow by taking the Fourier transform of equation (1.5) that 


S 2. &R(z, dy){e~#Y¥ — zg(t)} = e. (1.7) 


2. Funcetionals of Markoff Processes.—If x(t), t > 0, is a Markoff process with 
x(0) = xand P(2, E, t) = Prix(t + s) « E| x(s) = z} as its stationary transition 
probability, the results of Section | carry over pretty much the same, except that 
one uses Laplace transforms instead of generating functions. Suppose that for 
each fixed FE, P(x, F, t) is measurable in (x, 1) (with respect to ¥ and the Lebesgue 
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sets in 0 < ¢ < o) and that V(y) is real, measurable, and such that 
So E\| Vix(t))| | x(t) € BE} < ©, 0<t< ©, P(x, E,t)>0. 
Putting 


U(t) = Bc; V(a(r)) dr, 
R(x, E) = Jf,” ec “El e®"| 20) = a, a(t) € E} 
v(t, B) = R im e~ "P(x, E, t) dt, 
for s > 0, € real, we obtain the following theorem, analogous to Theorem 1: 
THroreM 2. The function R of equation (2.2) satisfies uniquely either of the two 
integral equations 
R(z, EF) = p,(z, BE) + 1 , i V(y) p(x, dy) R(y, EB), (2.4) 
R(x, E) = p(x, E) + i& SL Vy) R(a, dy) ply, E). (2.5) 


The proof is entirely analogous to that of Theorem 1. Equation (2.5) was given 
by Kac! for a symmetric stable process under different conditions and with a 
different method of proof. Blanc-Lapierre and Fortet? gave an analogue of 
equation (2.5) for a process with nonstationary increments in a somewhat different 
setting. However, it does not seem to have been noticed that R satisfies the 
pair of equations (2.4) and (2.5), and it is the interconnection between them (and 
the corresponding differential equations of the next section) which makes much of 
the theory simple and cohesive. 

Analogous to equation (1.6), we obtain for the function 


Q(z) = Ra, 2) = JS,” e~ “Ef e*®"’|2(0) = ax} dt 


the following integral equation 


| 
Q(x) = — + & S V(y) p(x, dy) Q(y). 


Theorem 2 shows that to each function V there is associated a unique R. Let 
now FR, and R, be associated with V;, V2. Then it is easily shown that 


R(x, E) — R(x, E) = i& S (Vely) — Vily)) Ril, dy) Rely, E), 


a formula which would be useful presumably in a perturbation calculation. 
Let a(t) be an additive process with E} exp(déxr(t))} = e*®. Then, analogous 
to equation (1.7), we obtain, by taking the Fourier transform of equation (2.4), 


S &'’ R(x, dy)(s — &V(y) — ¢(e)) = 


(ef. Kae,! who considered the case of a symmetric stable process, in which case 
g(t) = —|&/",0< a < 2.) 

3. Differential Equations Satisfied by R.—Let 2 = E", the n-dimensional Euclid- 
ean space, and suppose that P(z, E, t) has a density p(a, y, t), © = (UM, %, ..., 
tn), Y = (Yrs Yor -. + Yn) Over E”. Then in many cases p,(z, y) = So” e~“p(a, 
y, t) dt will be, as a function of x, the negative of the Green’s solution to Lp — sp = 
0, where 
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oO 


a; +> b , 

=2 * Ox Ox; ‘Ox; 

which vanishes for }>x,2 > © and is singular at x = y, a similar remark holding 
true for the adjoint equation L*p — sp = 0 as a function of y. 

Now if p has this property, it is easy to show that FR, considered as a measure 
over £", has a density r(v, y), and a purely formal application of the operator 
L — s to equation (2.4) and L* — s to equation (2.5) suggests that ris the negative 
of the corresponding Green’s solution to 

(L + 72V — s)r = 0, aig yo = ¥ieyY 4) 
(L* + 72&V — s)r = 0, y geog V = V (y).. (3:2) 


It seems difficult to find an existence proof for these differential equations (cf. 
Rosenblatt® for the Wiener process case), but it is possible to show that if p is the 
unique solution to (L — s)p = 0, (L* — s)p = 0, then the solution to equations 
(3.1) and (3.2), if it exists, is unique and is indeed the density corresponding to 
equation (2.2). 

Under similar circumstances the function g corresponding to equation (2.6) 
satisfies (L + 7&V — s)qg = —1, q bounded. 

4. An IIlustration.—-Let x(t) be the Uhlenbeck process; i.e., x(¢) is Gaussian, 
mean 0, stationary, Ltt with L(x @) a(t + s)) = exp (— | s| ). Put Viz) = 

osgnx, so that f°’ V(a(r)) dr + t/2 is the amount of time x(r7) is positive, 
0<7<t. This is i cae ad of the so-called “are sin law’ for the Wiener 
process. We have, namely, for this latter process Pr} Jy! '/2 sgn a(r) dr + t/2 < 
a} = (2/m) sin- 1 V7 a/t (see, e.g., Kac4 

In this case, setting x = 0, ¢ nit re .2) becomes 


d’r d 
(yr) + (/oté sgn y — s)r = 0, y ~ 0. 


dy? dy 


After some tedious calculations the Green’s solution can be obtained in terms of 
Weber functions, and we ultimately obtain 


co z : bg } ' s 
f r dy = q(0) = i) e —st EB} elf S35 2 sgn x (r)dr x(0) = 0} dt m rf e + 
0 - . 


/y it MEE tate Sate Sh Be 2+ 8/2 + it/4)P(8/2 — i€/4) 
8? + 1/4 &? T('/2 + 8/2 — 1&/4) (8/2 + 21€/4) + T('/2 + 8/2 + 21/4) (8/2 — 2&/4) 


In order to get the distribution explicitly, it is necessary to invert this last 
expression as a Laplace transform on s and a Fourier iar m on &, which 08 
very difficult. However, further calculations show that (1/+/t) fo! '/2 sgn a(r) dr 
is asymptotically normally distributed with mean 0 and variance '/» log 2,t—> @. 

We are preparing for publication some further examples of applied interest. 

* This work was initiated while the authors were at the RAND Corporation. The first-named 
author received subsequent support at the Statistical Research Center, University of Chicago, 
under sponsorship of the Statistics Branch, Office of Naval Research. The second-named author 


was a Guggenheim Fellow at the Institute for Advanced Study. 
1M. Kae, “On Some Connections between Probability Theory and Differential and Integral 
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COMBINATORIAL PROBLEMS IN THE THEORY OF GRAPHS. III 
By G. W. Forp anp G. E. UHLENBECK 
INSTITUTE FOR ADVANCED STUDY, PRINCETON, NEW JERSEY, AND DEPARTMENT OF PHYSICS, 
UNIVERSITY OF MICHIGAN 
Communicated May 9, 1956 


1. Introduction.—In two previous notes! we have presented solutions to a 
number of counting problems for linear graphs, Note I being devoted to graphs 
with labeled points and Note II to graphs with unlabeled points. The subject of 
this note will be the asymptotic behavior of these solutions for graphs in which 
the number of points is large. For definitions and notation we refer to Notes I 
and II. 

Although the solutions in Notes I and II were given in the form of functional 
relations between the various counting series, the counting series themselves 
were not considered as necessarily representing actual functions. Rather, the 
functional relations were looked upon as a convenient way in which to express 
the recursion relations between the counting-series coefficients. Certain of the 
star tree counting series, however, do represent functions analytic in some region 
about the origin, and in such instances we can use a knowledge of the analytic 
behavior of the functions to determine the asymptotic behavior of the counting- 
series coefficients. The methods used in this note are based on this idea and are 
a generalization of those used by Otter? for the case of Cayley trees. 

2. The Asymptotic Numbers for Star Trees with Labeled Points—Let T(p) be 
the number of rooted mixed star trees with p labeled points which are built out of 
stars from a given collection of types of star. Introducing the counting series 

a ” x? 
T(x) z T'(p) ot (1) 


we note that 7'(x) is equal to T(x; yi, ye, ...), defined by (I, 5), in which each 


Yq is set equal to unity, L.e., 
(2) 


Hence, from the implicit equation (I, 7) satisfied by T(r; y1, ye, ...), we see that 
T(x) is a solution of the implicit equation 


F(T, n= y —r7r exp [S’(T) | _— 0, 


where S’ = dS/dT, and 
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(4) 


in which q, 1s the number of points and s, the symmetry number of a star of type 
a, and the sum is over all types of star in the collection. If we for the moment 
assume that the sum over a is finite, i.e., the mixed star trees contain only a finite 
number of types of star, then /(7', x) is an entire function of T and of x. Since 
F(T, x) is zero for T = 0, 2 = 0, while OF /OT does not vanish for 7 = 0, x = 0, 
it follows from the implicit function theorem that relation (3) may be solved for 
T(x) which is analytic in some region about x = 0. The only singularities of 
T(x) in the finite plane will be at those values of x for which 
= = 1-—T7S8S"(T) = 0, 
oT 

where we have used relation (3) to eliminate x. In general, there will be many 
such values of 7’, but from the fact that S(z) is a polynomial with positive coefficients 
we can conclude that there is exactly one real positive 7’, which we denote by 7%, 
which satisfies equation (5), and that there are no 7’ satisfying equation (5) whose 
absolute value is less than 7». Since the coefficients in the power series (1) are 
all positive, it follows that, if xo is the radius of convergence of (1), then 7'(xo) = 7’. 
From expression (4) for S(x) we see that 0/07 will vanish for 7’ = 7) exp (70) 
only if 6 = 2x(m/d), m = 0,1, .«. d — 1, where d is the greatest common divisor 
of the set of integers [¢, — 1]. On the other hand, from the implicit equation (3) 
we see that 

T(xe'*/D) = gftalm/A T(z). (6) 
from which it follows that the only singularities of T(x) on the circle of convergence 
(a| = ao oceur at x = xo exp [i24(m/d)], m = 0,1, ...,d — 1. 

From equation (6) we see that the behavior of 7'(x) in the neighborhood of 
xX = 2X exp [i2x(m/d)| may be expressed in terms of the behavior in the neighbor- 
hood of « = x9; hence we confine our attention to the singularity at « = x. It 
can readily be shown that 


oF 
( a} i (S’’(To) + ToS’'’'(To)] ¥ 9, (7) 
OT?) 7 = "7, 


and hence the singularity of 7'(x) at x = 29 must be a branch point of order two. 
That is, in the neighborhood of x = 2x9, T(x) is a regular function of (a) — «)'” 
and may be expanded: 

T(x) = Ty — b(a — 2)" +.... (8) 
An expression for b may be obtained by differentiating relation (3) with respect 
to x and then taking the limit as x — x». The result is 


daqb?[S’'(T'o) + T)S'''(T5) | = T». (9) 


If we let ¢(p) be the number of free mixed star trees with p labeled points which 
can be constructed with stars from the given star collection and define the counting 
series ¢(x) similarly to equation (1), then we see from (I, 11) that 
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Pe d 
T'(2) x L(x), 
dav 


or, from (I, 14), 
tix) = T(x) — [T(a)S’(T(a)) — S(T(a))]. 
From these expressions we readily see that, near x = 29, 


To 2b 3/s 
tr) = Ua) — (to 2) em 2 


ro Ao 


Using the Cauchy integral formula to express 7'(p) and ¢(p) in terms of T(x) 
and t(x), respectively, and choosing for the contour a circle of radius r > x» indented 
around each of the cuts running from the branch points x» exp [¢24(m/d)] to in- 
finity, one can show that for large p the main contribution comes from the neighbor- 
hood of the branch points. Using equation (6), one sees that the contribution for 
the branch point 29 exp [727(m/d)] is the same as for the branch point x9 except 
for a factor exp [¢2mm(p — 1)/d]. The integral around 2p is readily evaluated by 
the method of steepest descent. Adding the contributions of all branch points 
then yields the following asymptotic expressions for 7'(p) and t(p): 

p—! 


not an integer, 
( 


p—! 


( 


an integer, 


l ry? 
_ Tp) ~ | 
p | 


' not an integer, 
_ ip) ~ b we ' (14) 
Pp: d x p + ” Pp / 2 ) - t . 
0 an integer. 
2V & d a 


As a simple application of these results, consider the case of pure Husimi trees 
in which the stars are m-gons and hence* 


| 
S(z) = i. (15) 
2m 


Using the general results above, we see that 


9 1/(m — 1) 9 1/(m — 1) 
T, = y ts 2 
» aati ’ to = ; 
(- _ ) (~, — =) 


9 m/(m — 1) 4 
bh: we (- - ) el/(m - sk (16) 


and, of course, d = m — 1. In this case an exact expression for T'(p) also may be 
obtained from (I, 8). 
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not an integer, 


p" p a 


l ; 
: = an integer. 
2"n! m— | 


One can readily verify that for large p expression (17) is asymptotically equal to 
expression (13) when the values (16) are substituted. 

A less trivial application is the case of arbitrary mixed Husimi trees. Here the 
collection of types of star is the collection of all polygons, including the line, so that 


(18) 


Although S(z) for this case is not a polynomial in z, it can be shown‘ that the 
discussion above still applies, provided we take 7) to be the smallest real positive 
value of 7 satisfying relation (5). After some numerical work, we obtain the 
results 


Ty) = 0.45631, X = 0.23874, b = 0.87170. 


In Table 1 the exact and asymptotic expressions for {(p) are compared for p = 1, 
59 MO: 
TABLE 1 
ComPARISON oF Exact AND AsyMpTOTIC EXPRESSIONS FOR ((p) FOR M1Ixep Husim1 TREES 
p 
1 2 3 4 5 6 

Exact U(p) l | 4 31 362 5,676 
Asymptotic ((p) I 1 3 28 332 5,297 

3. The Asymptotic Distribution of the Points Among the Slars in the Star 
Trees.—It is of interest to obtain some measure of the relative importance of a 
particular type of star in the collection of types of stars used to construct the 
star trees. One measure which can easily be calculated is the average fraction 
of the points in the star trees with p labeled points which lie on stars of a given 
topological type. To accomplish this, we introduce the counting series 


re) 
t,(z) = (q. — 1) E Ux; Yr, Yr, sl , (19) 
OY, y, = 1 


Yp = 


i.e., each y, is set equal to unity after the differentiation. From the definition! of 
(x; Yi, Yo, .--) it is clear that 


en, ee 
S’ de — 1)n, t(p; m, no, ... (20) 


t,(x) ; 
p=1 P> [np] 


is the counting series for the total number of points which lie on stars of type o 
in all the star trees with p labeled points.’ Using expression (I, 14) for ¢(x; 41, ye, 
.) and the implicit equation (I, 7) we readily see that 


| 
T(x). (21) 
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Introducing the expansion (8) for 7(x), we obtain the following expansion of ¢,(x 
in the neighborhood of x = 2»: 
Glen. Up, -- I/o 
t.(x) = t,(x%0) — = To" ~*b(ay — 2)" + .... 22) 


x 
go 


We define f,(p), the average fraction of the points which lie on stars of type o 
in the star trees with p labeled points, by® 


coeff x? : t,(2) 


: 23) 
(p — 1) coeff x”: t(x) ( 


f(p) = 
Asymptotically, using the results of section 2, we see that 


2 ali (24) 


GAdse — 1) 
fo(p) ~ = 


So 


If, as will nearly always be the case, 7) < 1, then we see from relation (24) that 
the distribution of the points among the stars in the star trees is weighted toward 
the stars with fewer points. As an example, the distribution for the case of arbi- 
trary mixed Husimi trees is given in Table 2. 
TABLE 2 
Asymprotic DistRIBUTION OF THE POINTS AMONG THE STARS IN Mrxep Husii TREES 
2 3 4 "a H * 7 i 8 

tm (per cent) 45.6 20.8 14.2 8.7 5.0 2.7 1.4 

4. The Asymptotic Numbers for Star Trees with Unlabeled Points —The dis- 
cussion in this section will be similar to that in Section 2 for the case of star trees 
with labeled points. Accordingly, we shall only sketch the arguments and, since 
the general formulas are rather complicated, restrict our attention to the case of 
mixed Husimi trees. The generalization to arbitrary mixed star trees is straight- 
forward. 

Let 0(p) be the number of rooted mixed Husimi trees with p labeled points. 
Introducing the counting series 


A(z) = Dd) O(p)z?, 
p=1 


we see that O(x)/x is equal to O(ys, ys, ...), defined by (II, 11) for the case of 
Husimi trees, in which each y,, is replaced by «”~', m being the number of points 
in an m-gon. That is, 

O(z) = O(a, 2% 4.) o4....-). 26) 
From the implicit equation (II, 12) satisfied by O(y, ys, ...), we see that O(x) 
is a solution of the implicit equation 


2 1 2 
’(0, xr) =0 — x exp ; > i bae tse o(e'))t = 0, 
k=1h 


m= 2 


where Z(D,,) is the cycle index of the dihedral group of order 2m,’? Z'(D,,) = 
(0/Of:;)Z(D,,), and Z(D,,, f(x)) denotes the function of x obtained by replacing the 
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variable f, by f(a*). In relation (27), 0(x") is considered as a function of x for all 
n > 1, so that (0, x) is an analytic function of © and x for |@| < 1 and |2| < 
& < 1, where & is the radius of convergence of the power series (25).* In order 
to show that £ is greater than zero, we note that 


coeff x”: Z’(D,,, O(a*)) < coeff x”: ef™ — D(z), (28) 


and, therefore, the power series (25) is bounded term by term by the power-series 
expansion of the function @(7) which is a solution of the implicit equation 

ns l 

$(0,z) =e-7 I = 0. (29) 

ang - oO 
Since the infinite product converges for |@) < 1, and since @ = 0 while 06/00 = | 
when 0 = Oand x = 0, it follows from the implicit function theorem that the radius 
of convergence of 6(x), which is less than or equal to  , is greater than zero. Hence 
& is greater than zero. 
The function 0(7) must have a singularity on the circle of convergence at « = £p, 

from which it follows that® 


ob | 
(*). , = 1-0) Y Z(D,, O(&)) = 0. (30) 


0 = Off) m=2 


Equations (30) and (27) together determine & and 0(&). Next, one can show 
that 0°?/00? ¥ 0 for x = & and 0 = O(&), and hence the singularity of O(x) at 
x = & is a branch point of order 2. That is, near + = &, 


O(x) = O(&) — Ble — x) + .... (31) 
An expression for 8 may be obtained by differentiating relation (27) with respect 


to x and taking the limit as z — &. The result is 


33%, >> [Z'’(Dn, O(ko) )+ O(&)Z’"’(D,,, O(Eo)) | = org + > D> gtO(é") x 
k=2 2 


m=2 2m=:; 


. > 
> &*(E)0) of Z' (Dns ace) | (32) 


> 


Z'"(Du, OE) +2 > 
k=1m = 
where 0 = d0/dx. 
If @(p) is the number of free mixed Husimi trees with p unlabeled points and 
@(x) is the corresponding counting series, then 6(2) may be expressed in terms of 
O(xr) by means of (II, 19): 


oo 


A(x) = O(x) — DO [O(x)Z’(D,,, O(2)) — Z(Dn, O(zx))]. (33) 


m 2 


Expanding both sides of equation (33) about + = &, we see that, because of equation 
(30), the coefficient of ( — a2) ’* vanishes, and we obtain 


A(x) ama 4(&) —_ (Eo — ZI + 6( —_ x) . + see (34) 


where 
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28 


Bey ; 


[ t > DY Hoee)Z’’ (Da, Os") 4 
k 2m 
ye 
k 


The asymptotic expressions for #(o) and O(p) are obtained exactly as in Section 
2. The results are 


aa &)"*O(Eo**) of, Z'(Dn, ac) | (35) 


I m= 2 


B P . 
O(p) ~ £, : sn (36) 
2V 4 
36 


6(p) =~ y 
1V 


a as (37) 


The numerical values are as follows: 
0.22215, O(&) = 0.45881, 6B = 4.395, 6 = 11.46. (38) 


In Table 3 the exact and asymptotic expressions for O(p) and @(p) are compared 
forp = 1, ..., 8. 
TABLE 3 
Comparison or Exact aNpD Asymptotic EXPRESSIONS FOR O(p) AND 0(p) FOR Mixep HusiM1 
TREES 

Pp 
3 5 7 8 
3 26 297 1,066 
3 30 365 1,347 
2 
I 


Exact O(p) l 
Asympt. 0(p) l 
Exact @(p) l 

l 


9 2: 63 18S 


Asympt. @(p) 5 : 45 146 


1G. W. Ford and G. E. Uhlenbeck, “Combinatorial Problems in the Theory of Graphs. I,” 
these PROCEEDINGS, 42, 122, 1956; G. W. Ford, R. Z. Norman, and G. E. Ublenbeck, ‘“‘Combi- 
natorial Problems in the Theory of Graphs. II,’’ these ProcEEDINGs, 42, 203, 1956. We shall 
refer to these as “Note I’’ and “Note II,’’ respectively, and equations in them will be referred to 
as “(I, 17),’’ ete. 

2 R. Otter, “The Number of Trees,’’ Ann. Math., 49, 583, 1948. 

’ The case m = 2, Cayley trees, is an exception. For Cayley trees S(z) = 27/2, Ty = 1, re = 
e~!, and 6? = 2e. 

4 The conditions on S(z) for which the discussion of Sec. 2 applies will be considered in the next 
note in this series. 

6 There is a certain arbitrariness as to whether the number of points associated with a star of 
type o in a star tree is dg Or dg — 1. The choice here is the more convenient. 

6 Note that, due to relation (I, 6), > fo(p) = 

o 


7 To be specific: 


l 9)/2 pee 
' | (irr ~ "" +,™"), mary 
I / 
Z(Dn) = » e(s) fs + 


2m om 


l = 
[scm = an m odd, 


where g(s) is the number of integers less than s and relatively prime to s; g(1) = 1. 

* Clearly & < 1, since the counting series 7(z) for graphs of labeled points is bounded by 0(2), 
so that f& < 29 <1. That @(p) > 7T(p)/p! follows from the fact that to every graph with un- 
labeled points there correspond s graphs with labeled points, where s is the number of symmetries 
of the graph. 

® The primes on Z always refer to differentiations after f,. 





COMPATIBLE SEMINORMS IN A VECTOR LATTICE* 
By CasPpER GOFFMAN 
UNIVERSITY OF OKLAHOMA 


Communicated by S. Bochner, May 16, 1956 


Let XY be a vector space. A locally convex topology 7 in X is given by a system of 
seminorms pg, ae A, for which X, is a Hausdorff space. We suppose that there is an 
order relation w in X for which X, is a vector lattice. We suppose, further, that the 
order dual, X,’, which consists of all bounded linear forms on X,, is a total subspace 
of the algebraic dual X’ of X, i.e., for every x ¥ 0 there is a bounded f such that 
f(x) # 0. In particular, this implies that X,, is Archimedean. A seminorm p is 
compatible with wif |x| 2 |y) implies p(x) = p(y). A locally convex topology r is 
compatible with w if its seminorms may be chosen so as to be compatible with w. 
For general information see N. Bourbaki, Eléments de mathématique.! 

Under the above assumptions, we observe that for every w there is a finest 7 
compatible with w. In particular, we show that there is at most one Banach 
topology compatible with w, and, if it exists, it is the finest locally convex topology 
compatible with w. Thus, for various function spaces, the norm in use is seen to 
yield the only possible Banach space compatible with the natural order relation of 


the space. 

We first note that if p is a seminorm compatible with w, then the continuity semi- 
norm set of p is contained in X,,’. For, if f(x) is a linear form on X such that, for 
some M, | f(x)| < Mp(2), for alla e X, then, for every y > 0, lz| <y implies p(z) S 


p(y), so that 
f(z)| < Mp(z) < Mp(y), 


which means that f(x) is bounded. We thus have the fact that if + is compatible 
with w, then X,’ ¢ X,’. 

We next note that there is a finest convex topology 7 compatible with w, namely, 
the one given by all the compatible seminorms. We point out that for this topology 
X,’=X,’. For, if f(x) « X, is positive, the function p(x) = f(| |) is aseminorm 
compatible with w for which f(z) is in the continuity seminorm set of p(x). But 
every bounded linear form on X, is the difference of two positive linear forms. 
Hence, for the finest topology + compatible with w, we have X,’ > X,,’,sothat X,’ 
= X,’. Fora total subspace Y ¢ X’, we call the finest locally convex space X, 
for which X,’ = Y a Mackey space, and summarize the above remarks in 

THEOREM |. Jf 7 is compatible with w, then X,’ ¢ X,’. If 7 is the finest convex 
topology compatible with w, then X,’ = X,,’. If r is compatible with w, X , is a Mackey 
space,and X ,’ = X.,’, then r is the finest convex topology compatible with w. 

In applying this theorem to Banach spaces, we first prove 

Lemma 1. Jf X,, is a vector lattice and r: |\x'| is a compatible norm for which X , is 
a Banach space, then, if \x,| 1s a Cauchy sequence of increasing positive elements, x = 
lim x, implies x = Xp, for every n. 


Proof: Suppose that (x, — «)+ > Oforsomen. Then (x, — x)+ = (x, — x)* for 
every p>n. But 
536 
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so that 
2 |li(r4, — x)t 


Since (2, — «+ > O, it follows that {2,} does not converge to x. 

THEorREM 2. Jf X,, is a vector lattice and r: | x is a compatible norm for which X, 
is a Banach space, then it is the finest convex topology compatible with X 
vector lattice may be normed as a Banach lattice in at most one way. 

Proof: Since X,is a Mackey space and X,’ ¢ X,,’ 


w? 


o»  Thusa 
we need only show that every 
positive f(z) on X,, belongs to X,’.. Suppose, then, that f(a) is not continuous but is 
— ves . ‘ » : | “s 
positive. Then, for every n, there is an x for which) x2) S land |f(x)| 2 n. Since, 


| 


by compatibility, | |2)) = 2) and since 
f(a+ +a-) = If(x z-)| = |f(a)}, 


it follows that there is an x > Ofor which |2) < Land | f(x)! =n. Now let {2z,} be 


a sequence of positive elements such that, for every n, 


l and S (Xn) 


l l l 
Yn = =X + Tot... t+ 


9) 92° on 


Then {y,} is an increasing sequence of positive elements. It isa Cauchy sequence 
since 


l l l 
Ynip — Yn — Un41 + see + Ln+p < 


Qn+l Qn+p On? 


for every n, p. Hence, if y = lim y,, then, by the lemma, y 2 y,, so that 


n> © 


. 3. RS, 
Stay > fF) > Ie) a 


n° 
- 


for every n, which is impossible. 

For all the spaces of real functions and sequences of analysis the order relation is 
the one for which the positive elements are those for which f(2) 2 0 everywhere, or 
almost everywhere, as the case may be, or x, 2 Oforalln. As examples of Theorem 
2, we mention only the two most important cases. 

Corouuary |. For the vector lattice of continuous functions on a compact set S the 
only compatible Banach space topology is gwen by 


= max | f(x) 


reS 


Corouuary 2. For the vector lattice L, of pth-power integrable functions, p 
on a measure space S, the only compatible Banach space topology is gwen by 


fl = (S| f@)| dul? 
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There are important vector lattices whose associated finest compatible convex 
topology is not a Banach space topology, and we hope to be able to discuss the gen- 
eral situation in the future. 


* Supported by National Science Foundation Grant NSF G-2267. 


iv, Paris, 1953; Livre VI: Intégration, chap. ii, Paris, 1952. 


A FORMULA FOR SEMISIMPLE LIE GROUPS 
By Harisu-CHANDRA 
INSTITUTE FOR ADVANCED STUDY 
Communicated by Marston Morse, June 26, 1956 


In order to save space we shall follow strictly the notation of an earlier paper.' 
Define fo, Po as usual,? and suppose that bo = bo A fo + bo A po. Then we say that 
hy is fundamental if bp 9 fy is maximal Abelian in fo. Any two fundamental Cartan 
subalgebras of qy are conjugate under G. 

THEeoreEM |. Let 6, be a connected component of hy’. Then there exists a real number 
c such that 

lim F (HH; O(m)) = ef(0) (H « bi) 
H->0 
for every f €C(qo). Moreover, c = 0 tf bo is not fundamental. 

THEOREM 2. Suppose that ho is fundamental and b, h, are all the distinct 
connected components of by’. Let c; be the real number of Theorem \ corresponding to 
b(l<j<r). Thenyn+t+at+...+c¢ 49. 

The following result plays an essential role in the proofs of the above theorems. 
Let dH denote the Euclidean measure on by and dk the Haar measure on the com- 
pact analytic subgroup Ky of Gy corresponding the fo. 

LEMMA 1. Assume that hy © fo, and for any g € C( ho) put 

g(X) = S exp (iB(X, kH))x(H)? & g(sH) dk dH (X € qo) 
KoX bo seW 
Then the integral 


F-(H) = (HH) S ce G(a*H)de* 


is convergent for H € bo’, and there exists a constant ¢ # 0 such that 


> «(s) Fs(sH’) =c Sy, » e(s) exp ((B(H’, sH))r(H)g(H) dH 
seW seW 
for all H’ € bo’ and g € C(bo). 

Let G be any connected Lie group whose Lie algebra is qo, and let A be the Cartan 
subgroup of G corresponding to bp. 

LemMa 2. Let & be the centralizer in G of an element ay ¢ A. Denote the natural 
mapping of Gon G = G/E by «> & (a@ eG). Then we can find a neighborhood B of 
ay in A with the following property. Given any compact set w in G, there exists a com- 


pact set Qin G satisfying the condition that if xar~' €w for some ae Band x € G, then 


Te. 
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By combining this with Theorem 2 of Paper 1, we can obtain certain results on 
G as follows. \ being an indeterminate, let D(x) (x « G) denote the coefficient of 
N in det {(A + 1)J — Ad(xr)}, where / is the rank and J the identity mapping of 
q. Let A’ be the set of those points a « A where D(a) # 0. We regard A’ as an 
open submanifold of the Lie group A. Note that GA = G/Ao = G*, and put 
a** = vax! (a € A, x € G), where x* denotes the coset +A in G*. Let C.°(@) be 
the set of all indefinitely differentiable functions on G which vanish outside a com- 
pact set. Then the integral 


oa) = |D(a)|"* fg f(a’) dx* 


converges for fe C."(G@) and a € A’, and @¢, 1s of class C” on A’. Let 8 be the uni- 
versal enveloping algebra of gq and the subalgebra of 8 generated by (1,6). Then 
elements of 8 and § can be regarded as differential operators on G and A, respec- 
tively. Let 3 denote the center of 8 and y the isomorphism of 3 into described 
in Paper 2 (p. 394). Then one can show that ¢,, = y(z)¢;(feC.°(@),ze3). Let us 
call a subset of G bounded if its closure is compact. Consider the space @o(A’) of all 
functions g on A’ of class C” satisfying the following two conditions: (1) g is zero 
(on A’) outside some bounded set; (2) 7,(g) = sup |g(a; v)| < © for every v eS. 
ae A’ 
We topologize C."(G@) in the usual way! and @)(A’) by means of the seminorms 
rT, (ve ®). 

THEOREM 3. gy € Co( A’) for f € CL" (G) and f > ¢, ts a continuous mapping of 
C.°(G) into @o(A’). Moreover, for any compact set w in G we can find a bounded 
subset B of A’ with the following property. If the carrier of f is contained in w, then 
ob, = O outside B. 

Corotuary. Let da denote the Haar measure on A. Then if g is a measurable 
function on A which is bounded on every compact set, the mapping T,: f—> J’ 4 go da 
(feC.(G)) is a distribution on G. 

The distributions of the form 7',, where g is a character of A, are closely related 
to the characters of G (see Paper 2). 

S() is isomorphic to under a mapping which preserves every element of 5, 
and so for any q e€ S(b) we get a differential operator on A, which we shall denote 
by O(q). Then Theorems 1 and 2 have the following analogues on G. 

TuroreM 4. Let A, be a connected component of A’ whose closure contains 1. 
Then there exists a constant c such that 

lim (a; O(m)) = cf(1) (a € Aj) 
ail 
for all f eC. (G). Moreover, ¢ = 0 if bo is not fundamental. 

THeorEM 5. Now suppose that bo is fundamental and A,,..., A; are all the dis- 
tinct components of A’ whose closures contain 1. Let ec; be the constant of Theorem 4 
corresponding to A; (1 <j <r). Then not every c; can be zero. 

Theorem 5 gives a simple formula for f(1) in terms of ¢,. In case G is compact, 
the proof of this formula is quite trivial. Moreover, I had proved it earlier® for (1) a 
complex semisimple group and (2) the 2 X 2 real unimodular group. More recently 
Gelfand and Graev® had extended it to the n X n real unimodular group. 

! These PROCEEDINGS, 42, 252-253, 1956. This will be referred to as “Paper 1.”’ 

2 See Trans. Am. Math. Soc., 75, 187, 1953. 
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3 See Bull. Am. Math. Soc., 61, 389-396, 1955. This will be referred to as ‘Paper 2.”’ 
4 See L. Schwartz, Théorie des distributions, Vol. 1 (Paris: Hermann & Cie, 1950). 

5 See Trans. Am. Math. Soc., 76, 514, 1954, and these PROCEEDINGS, 38, 339, 1952. 

6 Doklady Akad. Nauk S.S.S.R., N.S., 92, 461-464, 1953. 


ANALYTIC GROUPS OVER COMPLETE FIELDS* 
By Jun-tcut [Gusa 
DEPARTMENT OF MATHEMATICS, JOHNS HOPKINS UNIVERSITY 
Communicated by Oscar Zariski, June 14, 1956 


1. IJntroduction.—lf we accept the standpoint that real and complex fields are 
completions of algebraic number fields by their Archimedean valuations, analytie 
groups, for instance, over completions by non-Archimedean valuations appear with 
equal right as real and complex Lie groups. We can see immediately that most 
of the basic principles of the theory of Lie groups! can be carried over to analytic 
groups over an arbitrary complete field. We shall illustrate this by proving the 
existence of canonical co-ordinate systems on commutative analytic groups. This 
is a simple remark, but it is of significance because it contains Mattuck’s generaliza- 
tion of Lutz’s theorem? and the local isomorphism theorem of p-adic additive and 
multiplicative groups’ as special cases. 

2. Commutative Analytic Groups.—By a field we shall understand a field of char- 
acteristic 0. A complete field is a field which is complete with respect to the 
metric topology induced by a real-valued valuation of the field.4. We observe that 
fundamental existence theorems, for instance, the existence theorem for a completely 
integrable system, can be proved by the method of ‘calcul des limites’”® over 
an arbitrary complete field. It was not by chance that the notion of completeness 
was discovered by Cauchy in proving various existence theorems by this “calcul 
des limites.” 

Now, if k is a complete field, analytic manifolds over k, and hence analytic groups 
over k, can be defined as in the classical case. Also, if @ is an analytic group of 
dimension n over k, there are n independent forms of Maurer-Cartan, say w1,.. . , 
wn, on Y. We know that @ is commutative if and only if w; are closed for i = 
i,...,m.7 Weshall assume that © is commutative. 

Pick a system of co-ordinates, say (7, ..., %,), on @ at the neutral element e 
such that 2,(e) = O for? = 1,...,n. Then we can express w; in terms of the 


n 

co-ordinates 21, ..., 2, as w(x, dx) = }> A;,(x) dx; fori = 1,...,n. Here A;,,(x) 
j=1 

are convergent power series in the x, ... , 2», and we have det (A;,,(0)) # 0. Con- 
sider the system of differential equations 


Ou; 
= A;,;(x) 
Ou; 
for 7,7 = 1,.-.,n. Since w, are closed, this system is completely integrable. 
Therefore, u; can be expressed as convergent power series without constant terms 
in the 1, ...,2,. Since det(A,,(0)) # 0, we can take (wm, ..., u,) as a system of 
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co-ordinates on ® at «. Moreover, in terms of the co-ordinates uw, ..., u,, we have 
w,(u, du) = du, fort = 1,...,n. 

Now pick an element o9 of @ which is very close to «. The functions u;(oo-0) = 
v,(o) are defined and analytic in some neighborhood of ¢; hence there exist conver- 
gent power series f;(w4, ..., u,) such that v,(¢) = fi(m(o),..., un(o)) for? = 1,..., 
n. Since w; are invariant by translations, we get w;(v, dv) = w,(u, du), ie., Of:/ 
Ou; = 6; for7,7 = 1,...,n. This shows that f;(u) = a; + u; with some constants 
a, in k, and we have u;(o9) = vi(e) = f:(0) = a;, and hence 


U;(ao*o) = U;( Go) -+- u;(o) 


fora = 1,...,n. In other words, the mapping ¢ > (14(¢), ... , Un(o)) gives an 
analytic local isomorphism of @ to the n-dimensional vector group k” over k. 
If the valuation of k is Archimedean, we can extend this local isomorphism uniquely 
to a homomorphism of k”" into , and the proof depends only on the Archimedean 


character of the valuation. 

3. Mattuck’s Theorem.—Let G be an algebraic group variety of dimension n 
defined over a complete field &. If we denote by @ the group of rational points of 
G over k, by the theory of simple points,* @ becomes an analytic group of dimension 
noverk. Hence, if G is commutative, & is locally isomorphic with k". In particu- 
lar, if we take as G the multiplicative group of the universal domain containing k, 
we see that the multiplicative group of & is locally isomorphic to the additive 
group of k. On the other hand, if the valuation of & is non-Archimedean, k” con- 
tains small subgroups isomorphic to the n-dimensional vector group over the ring 
of integers of k. Furthermore, if k& is locally compact and nondiscrete, and if G 
is complete as a variety,’ G is compact; hence, by an elementary theorem on top- 
ological groups, any open subgroup is of finite index. Here the compactness of G 
follows from the facts that any complete variety is “dominated” by a projective 
variety” and that the set of rational points over / of a projective variety is com- 
pact. 

* This work was partially supported by the National Science Foundation. 

1C. Chevalley, Theory of Lie Groups, Vol. 1 (Princeton, N.J.: Princeton University Press, 
1946). 

2A. Mattuck, “Abelian Varieties over P-adic Ground Fields,’ Ann. Math., 62, 92-119, 1955. 

3 Cf. C. Chevalley, Class-Field Theory (Nagoya University, Japan, 1954), p. 59. 

* See, e.g., B. L. van der Waerden, Algebra, Vol. 1 (Berlin; Springer, 1955). chap. x. 

5 See, e.g., E. Goursat, A Course in Mathematical Analysis: Differential Equations (English 
trans. by Hedrick and Dunkel) (Boston, Mass.: Ginn & Co., 1917), pp. 51-53. 

6 Cf. C. Chevalley, Theory of Lie Groups, pp. 68-69, 100. 

7 Tbid., pp. 152-155. 

§ Q. Zariski, “The Concept of a Simple Point of an Abstract Algebraic Variety,’ Trans. Am. 
Math. Soc., 62, 1-52, 1947. 

9A. Weil, Variétés abéliennes et courbes algébriques (‘Actualités sci. et ind.,’’ No. 1064 [Paris; 
Hermann & Cie, 1948]), p. 25. 

” This remark is due to Chow. 





ABSTRACT HOMOTOPY. IV 
By DanieL M. Kan 


WEIZMANN INSTITUTE OF SCIENCE, REHOVOTH, ISRAEL, AND COLUMBIA UNIVERSITY, NEW YORK, 
NEW YORK 


Communicated by P. A, Smith, April 6, 1956 


1. IJntroduction.—The usual definition of the homotopy groups of a simplicial 
complex involves only its underlying topological space and disregards the simplicial 
structure. Our main result is the possibility of defining the homotopy groups of a 
simplicial complex using only its simplicial structure. The results will be stated in 
terms of ¢.s.s. complexes.' They clearly also apply to simplicial complexes, because 
every simplicial complex may be converted into a ¢.s.s. complex by a (partial) 
ordening of its vertices. 

The possibility of defining the homotopy groups of a simplicial complex in terms 
only of its simplicial structure was implicit in the results of a previous note.’ 
However, the definition implied by the present results is simpler and resembles the 
definition of the homology groups. 

2. C.s.s. Groups.*—A c.8.s. group G is a ¢.s.s. complex such that, for every integer 
n> 0: (a) G,, the set of its n-simplices, is a group; (b) the functions e': G,—~>G,_, 
and n': G,_;— G, are homomorphisms (written as operations on the right) for 0 < 7 
<n. 

2.1. EHvery c.s.s. group satisfies the extension axiom.” 

Let G be ae.s.s. group. Define for each integer n > 0a group G, ¢ G, by 

n-1 


~ 
G, = n_ kernel e' 
i—O0 


Then o €G, implies ce"e G, 1. Let the (right) homomorphism e": G, > G,_1 be 
defined by e" = e"G,. Then image @" is a normal subgroup of kernel é"~', i.e., G = 
iG,, "| is a (not necessarily abelian) chain complex. Its homology groups are 
H,,(G) = kernal #"/image #'*?. 

Let G be ac.s.s. group, and let 1 ¢G@, be the identity. Then m(G, 1) may be turned 
into a group, and, in view of 2.1, the homotopy groups of G may be defined,‘ and we 


have 

2.2. 2,(G,1) = H,(G) for each integer n > 0. 

3. Paths and Loops.—Let K be a c.s.s. complex, K, the set of its n-simplices. 
Let ¢ « Ky be a fixed 0-simplex, the base point of K, and let ¢” = gn... 9. If 
re K,,4;, then we denote by re,,; its (n + 1)th vertex re"... &. Let oe K,. An 
n-path of K ending at o is a sequence (a1, . . ., ¢2,) of (n + 1)-simplices of K such that 


gie"t! = 9", F2j-1€n41 = F2j€n4iy for 8 Ee 
o2,e"t'! = a, o2je"*! = oojie"*!, for bey 7. 

It is called a reduced n-path if in addition 
02541 F ony, 02) F ooj41, for all 7 


A (reduced) n-path ending at ¢” is called a (reduced) n-loop. 
A tree of K is a connected subcomplex 7 ¢ K which contains no reduced loops. 


542 
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A maximal tree T ¢ K isa tree with T) = Ko. Roughly speaking, a (maximal) tree 
isa connected |-dimensional subcomplex containing no closed loops (but. all vertices). 
4. The Homotopy Groups.—Let K be a connected ¢.s.s. complex with base point 
¢. We now define a ¢.s.s. group G as follows. Choose a maximal tree T ¢ K. 
Let G, be the free (non-Abelian) group with a generator @ for each o ¢ K,4, and a 
relation ¢ = 1 whenever + ¢ 7,4: 0r 7 = pn” for some p ¢ K,, and define the (right) 

homomorphisms e': G, — G,_; and n‘: G,1— G, by 
oe’, gn’ = on’, for0 <i<n; 

= ge"(ge"t!)—!, én” = an”. 

It is easily verified that G so defined is a ¢.s.s. group. 

As in Note III, let jk: K — S|K)| be the natural monomorphism of K into the 
simplicial singular complex of its geometrical realization. Let g, = (jK)¢. Then 
our main result is 

4.1. For each integer n > 0 there exists an isomorphism gx: ,(S\K\, ¢,) A Hy 
(G). This isomorphism is natural on the category of c.s.s. complexes with base point. 

Thus the homotopy groups of |K) may be expressed in terms of the free groups G,, 
and the homorphisms e': G, ~ G,_:. If K is finite in every dimension (in particu- 
lar, if K isa simplicial complex), then G,, is finitely generated for all n. 

5. The Homology Groups.—Define a c.s.s. group A as follows. Let A, = G, 
[G,,, G, |, and let the homomorphisms e' and 7‘ be those induced by the correspond- 
ing homomorphisms of G. Thus A is ““G made abelian.”” Denote by k: G — A the 
natural epimorphism and by k: G— A the induced chain map. Let he: 2,(S/K|, 
¢,) > H,(S'K|) be the Hurewicz homomorphism. Then we have, for each integer 
n> 0: 


5.1 There exists a natural isomorphism ax: H,(S\K\) @& H, (A). 
5.2. Commutativity holds in the diagram 


z,(SiK\, ¢.) —“— 8.10) 


hx re 
+ 


| 
BiSiK}|) ———» BaAAd 


6. The Maximal Tree T.—We define a c.s.s. complex B as follows. An n-simplex 
of Bisa pair (y, ¢), where yeG, andaeK,. Its faces and degeneracies are given by 
(y, ale’ = (ye'”, ce’), (y, «)n' = (yn' on’), for0 <i<n; 
(y, ae" = (ye'G, ce”), (y, o)n” = (yn", on”). 
Then the ¢.s.s. map p: B— K defined by p(y, ¢) = coisa fiber map. We identify G 
with the fiber of p over ¢ by setting y = (y, ¢”) for y €G,. 
For each reduced n-path b = (1, . . ., ox) of K ending at a, define t(b) « B, by 


t(b) = (y, o), where y = (G2j-1'Ge;). 
j=l 


6.1. The function t defines a one-to-one correspondence between the reduced n-paths 
(resp. reduced n-loops) of K and the n-simplices of B (resp. G). 
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If we identify the simplices of B with the reduced paths of K by this corre- 
spondence ¢, i.e., if we consider B as a complex with the reduced n-paths of K as its 
n-simplices and G as a ¢.s.s. group with the reduced n-loops of K as its n-simplices, 
then 

6.2. The fiber map p: B— K with fibre G over ¢ ts independent of the choice of the 
maximal tree T. Changing the maximal tree merely means changing the bases of the 
free groups G,,. 

7. The /somorphism gx.—Application of the functor’ Hx” to the fiber map p: 
B— K yields a fiber map Ex* p: Ex” B— Ex” K with fiber Ex” G over py = e’y¢. 
Clearly Ex” G is also a c.s.s. group. It may be shown that Ex” B is contractible, 
i.e., r,(Ex° B) = Oforall n. The exactness of the homotopy sequence of a fiber 
map then implies that Ox: r,(Ex” K,~) & m,1(Ex” G, 1) for all n. 

Let ¥, = (jEx” K)p. Then the isomorphism gs: ,(S|K|, ¢.) & H,-1(G) is de- 
fined as the composition of the isomorphisms in the following diagram. 


©_(Ex” G, 1) —*— «,(Ez” K, ¥) oe mn(S|Ex” K|, ps) 


~ | @ | 
le G* | sle K* 


y 1 
H,3(G) = ,-1(G, 1) 7, (S/K|, ¢s) 


8. The Groups G,.—Because G, ¢ G, and G, is free, it follows that G,, is free. 
‘ . ws » | " . 
A (in general infinite) set of generators of G, can be found as follows. For each inte- 
ger? with 0 <7 < n, define a subgroup G,,' ¢ G, by 


G,' = G, kernel €& 9 ... 9 kernel €'-'. 


It is easily verified that the function e'n' is a coset function of G,,'*+! in G,', i.e., that, 
for o, 7 €G,', or! €G,'+! implies ce'n' = re'n' and conversely. Asa set of genera- 
tors of G,° = G, is given, now iterated application of the Schreier theorem® yields 
a set of generators for G, = G,". 

Thus the problem of computing the homotopy groups of a finite simplicial complex 
is brought back to the group theoretical problem of computing the homology groups 
of a (non-Abelian) chain complex with a countably generated free group in every 
dimension. 


1 Cf. 8. Eilenberg and J. A. Zilber, Ann. Math., 51, 499-513, 1950. 

2D. M. Kan, these PROCEEDINGS, 42, 419, 1956. This will be referred to as ‘‘Note III.’’ 

3 The results of this section are due to J. C. Moore. 

4 Analogous to the definition of the homotopy groups of a cubical complex which satisfies the 
extension axiom; see D. M. Kan, these PRocEEDINGs, 41, 1092-1096, 1955. 

5 In this section we freely use the notation and results of Note III. 

6 Cf. A.G. Kurosch, Gruppentheorie ( Berlin, 1953). 





FOUNDATIONS OF QUASIGROUPS 
By SHERMAN K. STEIN 
Communicated by P. A. Smith, May 31, 1956 
A quasigroup is a multiplicative system xy = 2 satisfying three conditions: 
(i) x and y determine z, (ii) x and z determine y, (iii) y and z determine x. The 
most thoroughly studied quasigroups satisfy the condition of associativity 
ab-c = abe 


and are called groups. 

As observed by E. L. Post and others, conditions i, ii, and iii are symmetric; 
the operations of ‘solving’ and “evaluating”? play similar roles. Therefore 
quasigroups come in clusters of six. One could think of the six multiplica- 


‘ “ee 


tions ry = z, ye = 2, 7 = y, 2 = Y, ye = 2X, zy = x as conjugates of the 
original quasigroup xy = z and representing different points of view rather than 
different systems. However, if the original quasigroup Q satisfies an identity v, 
then the conjugates of Q may satisfy quite different identities which could be 
termed conjugates of v (see Theorem 1). If v is equal to its conjugates it is termed 
invariant. 

The operation of conjugation provides a method and point of view for the study 
of quasigroups and identities which illuminates much of the work on the axiomati- 
zation of groups by A. Tarski, M. Ward, and H. Furstenberg, for example, and sug- 
gests that the study of conditions other than associativity may be very fruitful. 
The following theorems illustrate this point of view. 

THEOREM 1. The identity ab-ac = be is a conjugate of associativity. 

The proof of this theorem is included because it is short and illustrates well the 
relation between conjugate quasigroups. The associative identity a-be = ab-c 
can be written: be = x, ax = y, ab = z imply zc = y. In the conjugate quasi- 
group obtained by interchanging the second and third letters, this is equivalent 
to: bx =c,ay = x, az = bimply zy = ¢, or simply az-ay = zy, which is equivalent 
to ab-ac = be. 

. Thus the theory of groups is equivalent to the study of systems satisfying the 
identity ab-ac = be. This fact suggests that systems satisfying other identities 
(even those usually in opposition to associativity) may merit study. 

THEOREM 2. A quasigroup satisfies the identity ab-be = ac if and only if it is a 


group which is a power of the group of order 2 or, equivalently, if all its conjugates are 


groups. 

In general, however, of all its conjugates, only the transpose of a group is again a 
group. 

THEOREM 3. There exists no finite quasigroup of order 4k + 2 satisfying the iden- 
tity avbe = ab-ac. For all other orders there do exist such quasigroups. 

The identity a-be = ab-ac or “left distributivity”’ has been studied by C. Burstin 
and W. Mayer, M. Knaster, and T. Mituhisa. 

THEOREM 4. A quasigroup of order 4k + 2 containing an element which is equal 
to its square cannot have a transitive group of automorphisms. 

Many quasigroups are homogeneous in the sense that their group of automor- 
phisms is transitive. The only homogeneous group is of order 1. This illustrates 
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a second symmetry which quasigroups can possess but which groups (with more 
than one element) lack. 

Definition: A quasigroup is medial if it satisfies the identity ab-cd = ac+bd. 
(Medial quasigroups, under different names, have been considered by I. M. H. 
Etherington, D. C. Murdoch, R. H. Bruck, M. Knaster, K. Toyoda, and T. 
Mituhisa.) 

TueoreM 5. Lach of the following conditions on Q is equivalent to the condition 
that Q be medial. 

1. For any algebra A and homomorphisms f, g: A — Q the product fg: A > Q is 
a homomorphism. 

2. Q satisfies an invariant identity of the form ab-cd = (t(a)+t(b)) (t(e)+t(d)), 
where t is a nonidentical permutation of a, b, ¢, d. 

3. Q satisfies (7) xa = by, za = bw imply xw = zy, and (ii) the simultaneous equa- 
tions ad+x = y+de,ab+x = y+be havea solution x, y. 

4. Qsatisfies: ab = a’b’,cd = c’'d' imply ac-bd = a'c’+b'd’. 

5. Qsatisfies: pq = rsimplies xp-yq = xr-ys. 

THEOREM 6, A smooth convex planar curve K is a conic tf and only if the follow- 
ing quasigroup ts medial. Let T be a fixed tangent to K with point of contact A. 
If a, b are in K — A, let ab be the point of contact of the tangent to K passing 
through the intersection of T with the line through a and b (define aa = a). 

Proofs of these theorems, along with a general development of the theory and 
implications of conjugation, will appear elsewhere. 


SOME MOLECULAR COLLISION INTEGRALS FOR POINT 
ATTRACTION AND REPULSION POTENTIALS* 
By M. A. Exrasonq, D. E. StoGryn, ANnp J. O. HirsCHFELDER 
UNIVERSITY OF WISCONSIN NAVAL RESEARCH LABORATORY 


Communicated June 1, 1956 


If the pairs of molecules in a dilute gas interact with an intermolecular energy 
of the form @ = d/r’, the collision integrals which determine the transport proper- 
ties can be written! 

Qius (7)"(5 ” an P +4 ) (2 ‘6)) {((6) 1 
YO) = | — d? ‘ 2 — ( A‘ (6). 
Qu LT (s (1) 
Here yu is the reduced mass of the colliding molecules. The following constants 
are evaluated: 
Repulsive interaction potential, d positive: 
A‘)(2) = 0.397601, A®(2) = 0.527843, 
A“(3) = 0.3115, A®(3) = 0.3533; 
Attractive interaction potential, d negative: 
A‘ (2) = 0.806907, A (2) = 0.710970, 
A®(3) = 0.6750, A®(3) = 0.4641, 
A(6) = 0.434, A®(6) = 0.328, 
A®(6) = 0.502, A“(6) = 0.420. 
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Although the A‘(6) have been previously tabulated and used for repulsive 
potentials with 6 > 4, the notion of the A‘ (6) for attractive potentials is new. It 
is somewhat surprising that these collision integrals should be finite. For large 
impact parameters the trajectory is deflected toward the center of force. As the 
impact parameter is decreased, the molecules tend to spiral about each other more 
and more, until, for a critical value of the impact parameter, true orbiting occurs. 
For an impact parameter slightly smaller than the critical value, the molecules 
spiral a few times about each other as they approach each other with a relative 
velocity which becomes progressively greater. Since our hypothetical molecules 
have no rigid cores, the two molecules pass through each other (r = 0) with infinite 
velocity and then separate, the outward leg of the trajectory being geometrically 
similar to the inward path. The calculations of the collision integrals for these 
attractive potentials had to be carried out numerically. 

The second virial coefficients for either 6 = 2 or 6 = 3 would be infinite,? but the 
transport coefficients are finite. The constants for 6 = 2 were readily evaluated 
in closed form. For the repulsive potential with 6 = 3 an expansion developed by 
Mott-Smith* was used. For the attractive potential with 6 = 3 the A‘(3) were 
evaluated numerically making use of elliptic integrals (of the first kind) for the 
angles of deflection. For é = 6 the A‘?(6) were taken directly from the Z“ ® inte- 
grals which Mason‘ tabulated for the Buckingham (6-exp) potential in the limit of 
low temperature where the exponential repulsive part of the potential does not 
affect the transport integrals. 

The collision integrals are defined in the following way.’ The collision of two 
molecules is considered as the equivalent one-body problem. A particle with re- 
duced mass u = myms/(m, + ms) approaches a fixed center of force with initial 
speed g and impact parameter b. The potential energy of interaction is given by 
& = d/r’, where d is a constant and r the distance of the incoming particle from 
the scattering center. Figure 1 shows the trajectory. The particle is deflected 
through an angle of deflection x, and the distance of closest approach is rp. 


h b | lI 6 
a. Yn = , Yo = ’ 
r , y bd |} 


Vm 


Now, if we define 


the angle of deflection may be written as a function of yo: 


ym(yo) ] 1 6) —'/: 
x(Yo) = © — 206, = © — 2f [ i (4) | dy. (3) 
0 5 \Yo 


The upper limit on the integral is the first positive root of the polynomial in the 
denominator of the integrand. 

Then the collision cross-sections for the various transport properties, in the case 
of repulsive potentials, are given by® 


éd 1" 
Qe = 2-| : | AM(6), (4) 


"/ omg? 
where the molecular collision integral, A‘? (6) is given by 


AM(6) = fo’ (1 — cos! x) yodyo. 
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Fic. 1.—A typical trajectory for the equivalent Fig, 2. 


A typical trajectory for the equivalent 
one-body problem with a repulsive potential. 


one-body problem with an attractive potential, 
in the case where 8 < Bri, and the particle passes 
through the center of force. 


When the potential is attractive, d = —a and the potential has the form 


Also, let us define 


(7) 


It is easily seen that with 8, which is a real, positive quantity, equation (3) becomes, 
for an attractive potential, 


ym(B) 1 67 -'/2 
x(B) = r — 26, = 7 — 2 f [ —-y? + (“) | dy. (8) 
0 6 B 


If 8 is greater than .rit, the polynomial in the denominator of the integrand has 
a first positive root, y,(8). However, if 8 is less than Berit, Ym(B) = ©, corre- 


sponding to the fact that the incoming particle actually passes through the center 
of force, r, = 0, and then goes out again. In this case the angle of deflection is 


given by 


a 1 y 694 -'/2 
(8) —260, = -2 f E —-yt+ ( ) | ly 
| " ou 


rather than by equation (8). Such a trajectory is shown in Figure 2. 
For the attractive potentials the various collision cross-sections are given by® 


ba \? 
Q” = 2s ( ) A(6), (10) 
2ug~ 


where 


AM(6) = fo’ (1 — cos! x)B dB. (11) 
I. CENTERS OF REPULSION 
For repulsive potentials, when 6 is 2 or greater, the A‘?(4) are finite. The case 


of 6 = 1, however, is a potential of sufficiently long range that the integral for A‘ (1) 
isinfinite. In this case, the angle of deflection is 
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x(yo) = 2sin~! [(1 + 4y2) 7"), (12) 


and the integral in equation (5) diverges. 

a) 6 = 2; Repulsive Potential.—The repulsive, inverse-square potential is of 
interest because the integrations may be carried out analytically.’ First of all, 
the angle of the deflection from equation (3) is 


ym l cid 
x(yo) = 3 — 2f [ — ( +; ) | dy. (13) 
0 2Y0” 


Let x = 2yo?; then y, = (1 + 1/x)7'” 


x(x) 2 (1 + =] “ Ty (1 + ‘Y | ( a) : 
x(x) 7 E _ ( + -) “| 


With this change of variables, equation (5) becomes 
AM(2) = '/4 fo (1 — cos! x) de. (15) 


Now let « = 2?/(1 — 2*); then x = w(1 — 2). Substituting in equation (15) and 


l sin wz 
ANG) © >"/y + . = d (16) 
4 Jo 1 


integrating by parts, 


_ r (' sin 2xz 
A®(2) = — f ~ dz. (17) 
0 


t 


ae 
But it is easy to show that 


1 sin n7rz iy as 
dz = '/.(—1)" [Sa(2nx) — 2Si(nz)], 
0 


l g7 
Z 


where® 
: r sin z 
S7(x) -f dz. (19) 
0 < 


A“ (2) = 0.397601, (20) 


Thus we obtain 


A® (2) = 0.527845. (21) 


b) 6 = 3; Repulsive Potential.—In the case of the inverse-cube repulsive potential, 
we were able to make use of a method which has been developed by Mott-Smith* 
to evaluate the A‘(6).. Through the application of Biirmann’s theorem® he has 
found that equation (5) can be written as 


AM(§) = 27? ‘fr (1 — cos! x) F(x) sin x dx, (22) 


where the /'(x) sin x is given by 
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. p—1-: 
F(x) sin x = cos * >, dD, (sin x) 


9 
- 2 =0 


The D, are constants, depending only on 6. 


9 


oo p! \dz’ , 2 ~/\sin '/oH(z) 


with the H(z) given by 


tc (ts s s\ I'('/.6t + '/.)T(—s + '/2 
Me) = - > ( ) #3 «-»'(') ( (—s+ 2) 
AY 


s=1 t=0 t T('/.6t — s + 1) 


For 6 = 3 we obtain the results given in Table 1. 


TABLE 1 
Tue Morr-Smira CorrricrENts Dp FoR 6 


Dp r Dy, 

0. 25438 3 —0. 00282 
0. 12487 4 0.00317 
0.01019 


Substitution of equation (23) into equation (22) and integration give the general 
expression 


nat l 
AM(6) = —2-74 DU p, & (-r2e (7) : 


p=0 n=1 


8p + 26n — 2/° 


n 


The two cases which are of interest are then 


: 
A(D(S — 94/3 
A(3) = 2 >>; D, F = | 


- l I 
A@(3) = 973 ) _ : 
“) Ly, Sige 


Thus we obtain 
A‘)(3) 0.3115, (29) 
A ()(3) 0.3533. (30) 
Il. CENTERS OF ATTRACTION 
When the potential is attractive, there is a singularity in the expression for the 


angle of deflection. The Mott-Smith formulation is no longer applicable, so we 
proceed as follows. 


There exists a critical value of b = b,,i, for any relative kinetic energy */2ug? such 
that for b less than b,,;, there are no positive roots to the equation 
fy) =1-y + ce = 0, (31) 


where 
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and for b = Derit, © = Cerit. This critical value may be obtained by noting that 
equation (31) has a double root at the critical value. Therefore, both f(y) and f’(y) 
vanish for the same positive root, Yerit. 


fy) = y[-2y? + cy] = 0. (33) 


One may combine equations (31) and (33), to solve for this root 


6 i 
Yerit = ey 


and then go back and obtain the critical value of c, 


se —5/2 9) 16/2 —-1 
Corit = 26 16 = 2] . 


6 1/2) — (1/6) 
=, = h/8 of 
Berit = 2 (, v3 ,) ‘ (36) 


Thus 


a) 6 = 2; Attractive Potential.—For the attractive case, with 6 = 2, the angle of 
deflection is, from equations (8) and (9), 


{1 — 1/287] ~1/2 1 —Il/, 


and 


Thus 


x(B) 


x(B) 


B<'/s (40) 


In spite of the infinite value of the angle of deflection for all values of 8? less 
than '/», one may still obtain the collision integrals. Equation (11) is split into two 
parts: 


AM(2) = [fF is (1 — cos! x)6B dB + fo 2 (1 — cos! x)B dB 
= G + A, (41) 


The second of these, the H, are obtained in closed form. Taking the expression 
for x from equation (39), and letting « = [1 — (28?) —'], it is found that 


I T l 

HY = ‘os? —~ ad . 
2 |, ie 2 (= = ) 
l : ) 

H® = { sin? mx d ( ). 
4 0 x — | 
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These expressions may be integrated by parts to give 


HW = T { sin mv asi 
ie l— x 


H® = — +i sin wee ; 
4 1 l — 


These may be rewritten as 


© Gi x *1 gi zt 
H® = oly wee dx “4 = = dare, 


4 \ 0 l— a 0 ] —_ 
- 7 ” sin 2rr “I sin 2rxr | | 
H® = +- dx + | dx; . 7) 
4 \ f l aah xz 0 ] _ x? f \ a 
The first integral in each case is found in Bierens de Haan," and the second is 
given by equation (18). Thus 
La | ( 


HO = ; y Leta) sin r — Si(r) cos r] + 1/2 Si(2x) — Sim) ¢, (48) 


H® =x T 4 = (Cie) sin 24 — Si(2m) cos 24] + '/2Si(4r) — Si(2m)f. (49) 


) 2 


Therefore, 


H® = , Si(2n) = 0.556907, 


H® = _ Side) 0.585970. 


The values of the integrals G” may be approximated in the following way. 
GY = 1/5 fo? (1 — cos! x)dB?. 


If the upper limit on the integral for the angle of deflection is taken to be an arbi- 
trarily large, but finite number z, 


AS A ee I i I 
aie ee 5A ) 1 35 7 ) lige f ies : i): 


With this expression for the angle of deflection, the integrand of equation (52) is a 
very rapidly oscillating function of 8°, oscillating between 0 and 2 for G and be- 
tween 0 and 1 forG®. As z becomes very large, the peaks become very close, and 
the average value of the integrand becomes very close to 1 for the G and to '/» 
for the G. Thus, in the limit that z becomes infinite, corresponding to the actual 
kinetic theory problem, 


GM = 


19 
7°?) 
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Thus we obtain 
A“(2) = G® + H™® = 0.806907, 
A (2) G® + H® = 0.710970. 


b) 6 = 3; Attractive Potential.—The calculations of the collision integrals for the 
inverse-cube attractive potential are somewhat more difficult.!' Figure 3 shows 


® 10 1) 12:13 14 15 6 17 18 19 20 29 22 23 ‘ or 


. tt 
geen?) ° 2734567 8€ 9 ON 219 4 15 te 


Fig. 3.—Angle of deflection y versus 8? for at- Fig. 4.—(B%rii — 87) versus 26, 
tractive potential with 6 = 3. when B < Berit. 
the angle of deflection x as a function of 6°. As in the previous case, there is a 
singular point. For 6 = 3, Berit = (*/4)'® = 0.95318 and y,,(8) varies from unity 
for B = © tO ¥m(Berit) = 3". Because of the singular point it is convenient to 
divide the integration of A‘(3) into two parts: 
AM(3) = C, (3) + €,(3), (58) 

where 

C(3) = | 2 fo. ae {1 —_ cos! (20) } dB?, (59) 

C(3) = Vo S314 [1 — (—1)! cos! (262)] dB. (60) 
Here 6, and @) are given in equations (9) and (8), respectively. The variation of 
Breit — B? as a function of 26, is shown in Figure 4. 

In evaluating C,‘0(3), it is convenient to change the variable of integration’? 
from 8? to q®, where q is minus the (only) real root of the polynomial in the de- 
nominator of equation (9) : 

B? = ¢g*(3(1 — q))-”. (61) 


Since the real root is negative in the range 0 < 8? < (°/4)', gis positive. Making 
the change of variables, 


q ‘/e x q? q =f 
q= ; f a ie dy, (62) 
lL—¢ 0 Lf. 1 


which is an incomplete elliptic integral of the first kind.'!* Thus 


Th ‘/s : 
i, = . a gF( ’ hk), 
l-¢¥ 
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where 
ws 
y?(3 — q’) 
a (3 — 24°) ((L — @) B — @)] 
cos! [q? — 2 + {(1 — 4g) B- @)}"). 


Since 


= ’ 


d(q*) ; (1 — q’?)” 


d(B?) gts (3 — q’) 


the first integral, Cy{0(3), is now given by 


3 — q’) 


l */ ( 
° 1 — cos! (26;) =, d(q?). 7 
2-3 J, | settles la — q’) A") (67) 


- 


C,(3) = 


Figures 5 and 6 show the integrands for / = 1 and for] = 2, respectively, as func- 
tions of q’. 

It is seen that the integrands oscillate wildly as q? approaches 0.75. The in- 
tegrals are subdivided into four ranges of g?. For the smaller values of g? numerical 
integration is satisfactory, making use of the elliptic functions. However, for 
the larger values of g? it is convenient to use 26, as the integration variable and de- 
velop approximate expressions for d(8*) /d(26,) which will allow analytic evaluation 
of the required integrals. Thus C,‘?(3) is divided into the following four ranges: 
(1) 0 < @ < 0.52, (2) 0.52 < q < 0.68, (3) 0.68 < q? < 0.748, and (4) 0.748 < 
q? < 0.75. In the first range, the integrands were evaluated for q? = 0, 0.04, 
O08; 0.12, ....; 0.40, 0.42, 0.44, .. . , 0.52, using four-point Lagrangian interpola- 
tions together with the elliptic integral tables of Byrd and Friedman.'* The con- 
tributions to Cy (3) and C,)(3) from this range were found to be 0.2449 and 0.1136, 
respectively. In the second region, an expression for d8?/d(26,) is obtained by dif- 
ferentiating the empirical equation 


Berit? — B? = 1.062129 — 0.0640403(20,) — 0.0092255(26,)? + 
0.0006563(26,)*. (68) 


The constants in this cubic equation were obtained by requiring that the values of 
8? for q? = 0.52, 0.58, 0.64, and 0.68 be the values calculated directly from the tables 
of elliptic integrals using four-point Lagrangian interpolation. Spot checks for 
different values of q? indicated that equation (68) gives values of Berit? — 8? accurate 
to within 0.05 per cent through the second range. The contributions to C)“(3) 
and (;(3) from this region are 0.1020 and 0.0690, respectively. In the third 
range, an expression for d8?/d(26,) is obtained by differentiating the empirical equa- 
tion 

Berit? — B® = (—30.419 + 7.5241(20,) — 0.21952(20,)? exp (—4,). (69) 
The functional form of equation (69) was suggested by the asymptotic behavior of 


B?, and the constants were adjusted so as to agree with the direct calculations for 
q? = 0.68, 0.727, and 0.748. Spot checks for different values of g? indicate that 
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equation (69) is accurate to within 0.5 per cent within the third range. The con- 
tributions to C,(?(3) and C,(3) from the third range are, respectively, 0.1217 
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‘1G. 5.—Integrand of C, (3) from equation (67). 
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Fic. 6.—Integrand of C, (3) from equation (67). 





and 0.0444. In the fourth region, use is made of the true asymptotic expansion of 
the elliptic function, F(¢, k). It is convenient to express F(¢, k) in the form™® 


F(¢g, k) = 2 K(k) — F(a — ¢,k). (70) 
As q? > 0.75, k > land g¢—~ 27/3. It is easy to show that throughout the range 
0.748 < q? < 0.75, only a negligible error is made in replacing F(r — ¢, k) by its 
limiting value F(2/3, 1) = In (2 — (3)').. When k is close to unity, K(k) can be 
expanded in the form'® 


t 4 ? 
Ki) = (1+) (7) -* See (71) 


where ¢ = (1 — k?)'. Thus, in the fourth range, we find that to within the re- 
quired accuracy 
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dp? 

d(26;) 

Using equation (72), we carried out the integrations analytically to obtain the 

contributions 0.0026 and 0.0022 to C\"?(3) and C;,(3), respectively, due to the 

fourth range. Thus, all together we find that C,((3) = 0.4712 and (,(2(3) = 
0.2292. 

To evaluate the Co‘0(3), as defined in equation (60), it is convenient to change the 

variable of integration from 6? to p*®, where p is the first positive root of the denomi- 


= 36(48)'* (2 — 3°”) exp(—@) = 17.528 exp (—4,). (72) 


nator in equation (8). 
B* = p[3(p? — 1))-. 


Since y,,(8) = p, it follows that 


3 ‘/s p 3 3 ate f 
6, = ( - ) | | - PH ——P_ yy | dy, 
g- I oi” 2 = 1 


which is another elliptic integral of the first kind.'"7. Thus 


where 

g’ = 2[p(p? — 1)-' (4p? — 3)'7]-", 
k’? = 1/, + 1/o(2p? — 3) (4p? — 3)~ 
joing sfeePocetitt cans an 


a 2(3 — p?) 
Equation (60) for the C.°0(3) becomes 


(3 — p*) 


l 3 ‘ 
CL (3) = ; ; - (1 — CR arene 1°  d(p*) (79) 
= 1 


26 
Figures 7 and 8 show the integrands for 1 = 1 and / = 2, respectively, as functions 
of p?. 

Here again, the integrations were carried out in several parts, with the method 
chosen to fit the behavior of the integrand in the various ranges of p?. The four 
ranges are (1) 1 < p? < 1.025, (2) 1.025 < p? < 1.050, (3) 1.050 < p? < 2.90, 
and (4) 2.90 < p? < 3.00. In the first region, one can express cos x in a power- 
series expansion, 


2 


1 — cos x = (" ) (p — 1)? + 239.5 (p — 1)? +..., (80) 


where the constant 239.5 was obtained by requiring that, when p? = 1.025, equa- 
tion (80) gives the same angle of deflection as is obtained by the use of the tables of 
elliptic integrals. Making use of equation (80), the integrals over the first range 
were evaluated analytically to give 0.0015 and 0.0029 as the contributions to C2? (3) 
and C,°)(3), respectively. In order to evaluate the contributions from regions 2 
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and 3, the integrand of equation (79) was evaluated for p? equal to 1.025, 1.050, 
1.075 1.200, 1.250, 1. 300, 1.4, 1.5,..., 2.9. First, a four-point Lagrangian 
integration formula was used which corresponded to assuming that the integrand 
is the cubic polynomial in p? which fiis the calculated values at p? = 1.025, 1.050, 
1.075, and 1.100 and integrating this cubic from p? = 1.025 to p? = 1.050. This 
gave 0.0026 and 0.0052 as the contributions from the second range to C2‘?(3) and 
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The contributions to these two integrals from region 3 were 


C, (3), respectively. 
In region 4, cos (262) oscillates 


found by Simpson’s rule to be 0.1996 and 0.2267. 
rapidly, but the maximum value of the integrand is so small that it is easy to make 


a sufficiently accurate approximation to the contributions, Js? and J4, to the 
The largest possible value that J,‘ might have is obtained 
Thus 


two collision integrals. 
by setting cos (262) = 1 in the integrand of equation (79). 
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3 
IO Clow ae f (3 — p®) (p? — 1)7°? d(p?) = 0.00026. (81) 
2.9 


Since the cosine term of equation (79) has infinitely many oscillations in this 
region, to a sufficient degree of accuracy, the factor (1 + cos 26.) can be replaced 
by its mean value of unity in the integral for J,. Thus a good estimate of J,” 
is 1/9 Jmax‘?, and the contribution of region 4 to C2(3) is taken to be 0.00013. 

The best estimate of J, is obtained in a similar fashion. Here the maximum 
value of (1 — cos? 26) is unity, so that the maximum value of Jy = '/oJ max !. 
The average value (1 — cos? 24) is one-half, so that our best estimate of J,‘ is 
'/sJ max‘), or the contribution of region 4 to C23) is taken to be 0.00007. 

Thus, all together, we find that C."9(3) = 0.2088 and C,@ (3) = 0.2349. Com- 
bining these values with those of the previously obtained C,?(3) and C\(3), we 
finally obtain the results 

A(3) = 0.6750,  A®(3) = 04641. 


) 6 = 6; Attractive Potential.—The collision integrals for an attractive point- 
force power-law potential with 6 = 6 can be estimated from Mason’s‘ calculations 
of the collision integrals, Z“ ©, for the Buckingham (6-exp) potential in the limit 
that the temperature approaches zero. In this low-temperature limit, the expo- 
nential repulsive part of the Buckingham potential does not affect the values of 
the collision integrals. Thus we can consider the Buckingham potential to be 
equivalent to the inverse sixth-power attractive potential, @(r) = —a/r* witha = 
rm® €/(1 — (6/a@)). Interpreted in this fashion, the Z* for 7’ = 0 can be ex- 
pressed in terms of the A‘(6) by the relation 

_ 4(6)” ‘(+ DI(s + (5/3))A (6) | 


Zi s — ee Sa 9 
(s+ 1d +21 — (-1)) -) 


Thus, from the tabulated values of Z“ ® for 7’ = 0, we obtain 
A“(6) = 0.434, A® (6) = 0.328, A®(6) = 0.502, A‘ (6) = 0.420. 
Mason estimates that these integrals may be in error by 1 per cent. 


* This work is being done under Office of Naval Research Contract N7-onr-28511. 
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60%) as given in eq. (1) is obtained from Q by use of the equation 


kT \' Paar 
27, 0 


where y? = wg?/2kT (Hirschfelder et al., op. cit., p. 525). 

’ The evaluation of A“(2) and A‘(2) was carried out by John S. Dahler with the help of 
Donald W. Jepsen. 

§ Tables of Sine, Cosine, and Exponential Integrals, Vol. 2(New York: Federal Works Ageney 
Works Progress Administration for the City of New York, 1940). 

°K. T. Whittaker and G. N. Watson, A Course of Modern Analysis (4th ed.; Cambridge: At 
the University Press, 1952), p. 128. 

” TD. Bierens de Haan, Nouvelles tables d’intégrales définies (New York: G. E. Stechert & Co.., 
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* sin px bis > ; 7 : 
‘ dx = — ,Ci(pq) sin pg — Si(pq) cos pq}, 
oe Fr qd 


with 
Ci(pq )= 
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(Berlin: J. Springer, 1954), p. 88, eq. (241.00). 
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ON THE APPROACH TO THE FINAL PERIOD OF DECAY IN 
ISOTROPIC TURBULENCE ACCORDING TO 
HEISENBERGS TRANSFER THEORY 


By W. H. Retp 


BALLISTIC RESEARCH LABORATORIES, ABERDEEN PROVING GROUND, MARYLAND 
Communicated by S. Chandrasekhar, June 9, 1956 


Introduction.—The final period of decay is usually defined! as that period in the 
decay of the turbulence for which the nonlinear terms in the Navier-Stokes equa- 
tion are negligible. Although this formulation of the problem does give useful in- 
formation about the asymptotic form of the energy spectrum and the corresponding 
correlation functions, it cannot, by its very nature, deal with the approach to this 
asymptotic state during which the transfer of energy is not quite negligible. 
Furthermore, at the onset of the final period, the Reynolds number of the turbu- 
lence, though small, is not negligible, being of the order of 5 or so,? and this suggests 
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that a transfer of energy may exist which is compatible with the observed final- 
period decay law, 


u2 « (t — &)~”. (1) 


Thus, if one regards the final period of decay as being defined by the decay law (1), 
rather than by the absence of energy transfer, then one can obtain not only the 
known properties of the energy spectrum in the final period but also, for a given 
transfer theory, the asymptotic form of the transfer function. 

For the decay law (1), the dissipation length parameter is given by 


A? = 4r(t — to), (2) 
so that the Reynolds number of the turbulence, R, = (u?)'“A/v, is not constant 


in the final period as it is in the initial period but decreases according to the law 
Ry « (t — &)~”. (3) 
Final Period Similarity.—Consider now the usual (exact) dynamical equation for 
the rate of change with time of the energy spectrum function /(/:) in isotropic tur- 
bulence, 
Oh (hk) 
ol 


= T(k) — 2vk? E(k), (4) 


where 7'(/) is the transfer function usually denoted by this symbol. For purposes 
of the present discussion we shall use Heisenberg’s physical transfer theory*® to 
supplement equation (4); this transfer theory leads to an expression for 7'(k) in 
terms of £(k) which can be written in the form 


" d F E(k’) k 
7 k = —? , 4 Is! F > y . ” . , 
, , dk ra \ k’3 dh 3 h I (A ) dk ° 


where « is a constant of order unity. 
We now assume a similarity solution of the form 


E(k, t) = y *(t = 7 ae : R,? F(a, R,), 
and : : 
EE) = 9 (t— bh) Re UR), 
where 

x = (vy k*t)'”. (7) 


This type of similarity remains valid whatever may be the dependence of R, on the 
time and thus includes not only initial period similarity (2, = constant), which I 
have discussed elsewhere‘ along similar lines for small values of the Reynolds 
number, but also final period similarity with which we are concerned here. 

Thus, by substituting equations (6) in equations (4) and (5), and by using the 
fact that Ry « (t — to) ~’/*, we obtain 


S..2 ro) 
se ae ee 7 — 977 
(: ae =) F(x) + 4(a? — 1) F(x) = 2U(2) (8) 
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and 


| 
’ ) d | F(x’) , . wo ” ” 
U(r) = — 2«R, \ — dx 2"? F(x") dx". (9) 
dx J z’* 0 


The Asymptotic Solutions.—¥or small values of the Reynolds number, equation 
(9) shows that U(x) is of smaller order than F(x). Accordingly, we can assume a 
solution in the form of a power series in the Reynolds number of the form 


F(a, Ry) = > F,(x)(«R,)", 


n=0 
and 


U(z, R,) = > U,(z)(«R,)". 


n=1 


Since Ry, « (tf — ty) “‘, this power series expansion in the Reynolds number is 
equivalent to an asymptotic expansion in the time, and it is this fact which makes a 
solution of the form (10) meaningful in the final period of decay; for, with an in- 
creasing Reynolds number, we are simply moving backward in time. 

By substituting the series (10) into equations (8) and (9) and by equating like 
powers of R,, one is led to a system of equations from which the functions F(x) 


and U’,(x) can be obtained by quadrature. From the relation 


u? = 2/3 J,” E(k) dk, (11) 
the boundary conditions on the functions F’,(2) follow in the form 


na (3/8 (n=0), 
F(x) dx = - 12) 
Si ae oe ; 


which show that the general shape of /'(.x) is determined by F(x) and that the 
functions F’,(x), for n > 1, cause only a detailed modification of that shape as the 
Reynolds number is increased and leave the decay law (1) unaffected. 

The functions Fo(2) and U\(2) evaluated in this way are given by 


— 2x2 


F(x) — 4(2 7) ae e 


. 2 ‘d r p/s — , = we —2x"2 
U(x) = —16 alas Raa oe ae (14) 
T dx J x 0 


which are shown graphically in Figure 1. 
The Skewness Factor.—From the foregoing solution it is now possible to calculate 
the skewness factor of —O0u,/O2, 


and 


(Ou Or1)* 


Ve 15) 
[ (ou, dz,)? | i iis 


S=- 


What is perhaps most surprising about this final-period solution is that S does not 
tend to zero asymptotically in the final period as one might have expected but ap- 
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Fic. 1.-The distribution of energy transfer in the final period of decay according 

to Heisenberg’s transfer theory. Curve 1, the energy spectrum, Fo(.r); curve 2, the 

energy-dissipation spectrum, 0.805427 Fo(x); curve 3, the transfer function, 2U7,(x). 


proaches a nonzero constant value. This can easily be seen by writing equation 
(15) in the alternative form 


me . (30)? So k?T(k) dk 


[ fo” k? E(k) dk} 
In terms of the similarity functions, equation (16) becomes 


s = 2 (30) Sy aa U@) dx 
I4 R| So x? F(x) dx} 


which shows that, as R, > 0, 
S/« = 128/35 f5° 2? Ui(x) dx + O(«R,). (18) 


A numerical evaluation of the integral occurring in this equation then shows that, 
as Ry, > 0, 


—> 0.5764. (19) 
° 

Thus, provided that « is not Reynolds number dependent—a possibility that can 
hardly be admitted—the skewness factor approaches a constant value in the final 
period. 

It is of some interest to compare the value of S/x obtained here in the final period 
of decay with the initial period values‘ of 0.78 (R, ~ 0) and 1.52 (R, > ©). These 
results suggest, then, that the value of S/x varies by less than a factor of 3 through- 
out the entire life history of the turbulence, whatever may be the value of the 
Reynolds number. 

Finally, it should be remarked that in the final period the individual Fourier 
coefficients of the velocity field, though statistically dependent, do not interact 
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but simply decay through direct viscous action.2 The present calculation then 
shows explicitly that this statistical dependence between the different Fourier co- 
efficients does exist and persists but that, for Heisenberg’s theory, it has a universal 
character determined only by the transfer mechanism. 

'G. K. Batchelor, The Theory of Homogeneous Turbulence (Cambridge: At the University 
Press, 1954). 

2G. K. Batchelor and A. A. Townsend, ‘Decay of Turbulence in the Final Period of Deeay,”’ 
Proc. Roy. Soc. London, A, 194, 527-543, 1948. 

3 W. Heisenberg, “Zur statistischen Theorie der Turbulenz,’’ Z. Physik, 124, 628-657, 1948. 

‘'W.H. Reid, “Two Remarks on Heisenberg’s Theory of Isotropic Turbulence,’ Quart. Appl. 
Math. 14, 201-205, 1956. 

5 Cf. M. J. Lighthill, Nature, 173, 746, 1954. I am indebted to Professor Lighthill for some 
further illuminating remarks regarding this point. 


ON THE CLASSIFICATION OF REINFORCEMENT SCHEDULES* 
By W. N. ScHOENFELD, W. W. CumMING, AND E. Hearst 
COLUMBIA UNIVERSITY 
Communicated by C. H. Graham, June 11, 1956 


Schedules of reinforcement for operant behavior have ordinarily been classified 
under two major headings: those based on the passage of specified time intervals 
of nonreinforcement between reinforced responses and those based on the emission 


of specified numbers of responses by the organism. Skinner describes this dis- 


tinction as one between schedules which “are arranged by a system outside the 
organism [fixed and variable interval schedules] and those which are controlled 
by the behavior itself [fixed and variable ratio schedules].’’! Both types of schedule 
have in common a periodicity or intermittence of reinforcement, but some of the 
differences in their behavioral effects are striking.? In view of the differential 
effects observed on behavior and the different procedures followed by the experi- 
menter, the two categories are regarded as distinct by most, if not all, theorists. 
Nevertheless, it would seem worth while if such an integration could be made, 
and an attempt in this direction is offered in the present paper. 

We may consider that a schedule of reinforcement acts to maintain behavior 
through the differential reinforcement of a particular pattern of responses in time 
(or, equivalently, single responses having certain temporal characteristics with 
respect to other responses). From this viewpoint the ‘‘count” involved in a 
“ratio”’ schedule may be only incidental to the generation of a maximal rate of 
responding, with the crucial factor being the increased probability of reinforcement 
for responses following each other at short intervals. By implication, then, if 
external conditions were arranged to favor rapid responding by the organism, we 
might expect to observe the “bursts” (short periods of high response rate) and 
“breaks” (pauses in responding) so characteristic of “ratio” reinforcement. “In- 
terval” schedules, on the other hand, may differentially reinforce responses pre- 
ceded by relatively longer intervals of no responding, producing the response 
rates and temporal distributions characteristic of these schedules. Skinner has 
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noted that “schedules are simply rather inaccurate ways of reinforcing rates of 
responding. If these schedules are so regarded, there is no reason to believe 
that their integration into one conceptual framework is impossible. 

The projected classification of reinforcement schedules utilizes “time-sampling”’ 
contingencies, in the description of which we shall employ the terms ¢? and ¢°. 
The term ¢? represents a period of time during which a specified instrumental 
response of the organism may be followed by a reinforcing event; ¢* represents a 
period of time during which this response will not be followed by a reinforcing 
event. A fundamental case may be established with these three restrictions: 
(1) ¢? and ¢* are held constant, (2) ¢? and ¢° are alternated, and (3) only the first 
response in ¢? is reinforced. It may be noted that reinforcements can easily be 
“missed”? on such a schedule; if the organism does not respond at all during the 
?? period, it does not receive a reinforcement during that t? + ¢* cycle. 

The experimental domain suggested by these variables may be visualized as a 

co-ordinate system with the duration 

+0 of ¢? along one side and the duration 

40 A of t* along the other. Increasing either 

uti : t? or t* alone changes both the total 

cycle length, ¢? + ¢*, and the propor- 
tion of the cycle during which a re- 
sponse may be reinforced, ¢?/(t? + ¢*). 
The experimental domain may then be 
mapped as in Figure 1, wherein are 


993 





schematized certain sectors which may 
represent different behavioral effects 
correlated with classically defined ‘in- 





terval” and “ratio”? schedules. It 
should be remembered that the areas 
and boundaries of the sectors are 
merely suggestive; their existence and 
limits (probably not sharp) remain to 
be empirically determined in terms of 
the behavioral properties exhibited. 

Fixed Interval Schedules.—This case is marked by ¢?/(t? + ¢*) approximating 
or equaling unity (where {* > 0 and t? > 0) and comparatively long cycle lengths. 
If ¢?/(t? + t*) = 1, the organism can “miss” a reinforcement only in the sense 
that a whole cycle may pass without a response, with one reinforcement opportunity 
lost, though the next subsequent response is reinforced whenever it occurs. When 
0< t?/(t? + t*) < 1, “missing” of a reinforcement can take place in the same way, 
but it can also occur when the first response in a cycle is made in ¢(*. With ¢? + ¢ 
long enough, these schedules permit the development of a temporal discrimination, 
as shown by a varying “‘density”’ of response in which the organism’s interresponse 
intervals do not remain random but systematically decrease in any cycle. These 
schedules are schematically located in sector A of Figure 1. 














Fic. 1.—For explanation see text. 


Regular Reinforcement and Extinction.—‘‘Regular reinforcement” is obtained 
with ¢* = 0 and with cycle lengths (composed entirely of ¢?) shorter than the 
duration of a reinforcement. Thus, if a pigeon is allowed 3 seconds’ access to 
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grain as a reinforcement, a ¢? equal to 2 seconds and ¢* = 0 will insure that a new 
cycle has begun by the time the bird has finished eating, and every response will 
be reinforced (see sector B, Fig. 1). 

“Extinction” is the case where t? > 0 and t* > 0 (t?/(t? + t*) + 0). The 
subject is reinforced so infrequently that behavior cannot be maintained (see 
sector C, Fig. 1). 

Ratio Schedules. The incorporation of response contingency (‘‘ratio’’) schedules 
into the proposed framework apparently presents greater problems than the above, 
since the co-ordinates of Figure 1 are temporal in nature. Several considerations, 
however, help to resolve the apparent incompatibility between the two metrics 
of time and response number. By reducing ¢?/(t? + ¢*) to a sufficiently small 
value and keeping reasonably long cycle lengths, we favor the adoption by the 
organism of a high response rate, since on such a schedule reinforcement is more 
probable or frequent following short interresponse intervals (i.e., high rates). 
Sector D of Figure 1, by such reasoning, would include ¢? and ¢* values which might 
be expected to yield behavioral effects typical of fixed ratio and some variable 
ratio schedules. It- is an empirical question whether choices of such ¢? and ¢* 
values would result in “ratio”? behavior (‘“‘breaks” and ‘‘bursts,”’ ete.) and whether 
sudden large decreases in ¢?/(t? + t*) would reproduce the known effects? of 
switching an organism to a much higher fixed or mean variable ratio. Some pre- 
liminary data will be presented below indicating that such is, indeed, the outcome. 

Certain variable ratio schedules, which might be termed ‘“‘random ratio,” specify 
the probability of reinforcement for any particular response; a 0.25 random ratio, 
for example, indicates that there is a probability of 0.25 that any single response 
will be reinforced. The rates of responding generated by such schedules are par- 
ticularly stable and are free of the “breaks” after reinforcement.‘ In the present 
framework, as cycle lengths become shorter, such rapid alternations of t? and ¢* 
are eventually reached that ‘breaks’ after reinforcements should disappear, be- 
cause responses immediately after reinforcement are as likely to be reinforced as 
are any other responses. In this case, the probability of reinforcement for any 
response, as with “random ratio’? schedules, can be specified by the proportion 
t?/(t? + ¢*). For example, with short cycle lengths, perhaps of the order of less 
than a second, and with ¢?/(t? + t*) equal to !/2, there is on the average one 
chance in twenty of a given response’s being reinforced. Sector B, “regular re- 
inforcement,” fits into sector #, ‘‘random ratio reinforcement,’’ and accounts for 
cases where the probability of reinforcement is unity for each response. Sector 
C’, “extinction,” is also included in sector EF and covers the case where there is a 
probability of zero that any response be reinforced. Some preliminary results 
displaying ‘‘random ratio” response features, for schedules assignable to sector FE, 
will be presented below. 

It may be noted that “breaks” after reinforcement are observed in high “‘fixed 
ratio” schedules and ‘‘variable ratio”? schedules with a high mean ratio but become 
negligible during random and low mean variable ratio schedules of reinforcement. 
We would expect a similar difference to exist between sectors D and FE of Figure 1; 
in sector D (presumably fixed and high mean variable ratio schedules) the cycle 
length is long enough so that responses immediately after reinforcements are 


almost never reinforced, while in sector /, as remarked earlier, the short cycle 
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length makes “breaks” after reinforcement unlikely. Additionally, we might 
expect that once the cycle length goes below a critical small value, cycle length 
should cease to be an important variable; the organism’s responding is then 
affected only by the ¢? (t? + ¢*) value prevailing. 

The areas of Figure 1 as yet unmentioned cannot now be categorized in any 
systematic way with regard to their behavioral effects, and they may represent 
schedules of reinforcement not falling within conventional categories. Empirical 
determinations in these areas may indicate a blending of ‘interval’ and “ratio”’ 
behavior. 

Other Schedules.—The three restrictions underlying Figure 1 are not the only 
possibilities. Each represents a special case of one or more variables which might 
themselves be systematically manipulated. For example, a schedule may be con- 
structed in which ¢? and ¢* are not held constant but are varied according to some 
program consisting, perhaps, of a prearranged contingency table governing the 
probability at any time of passing from one part of the domain to another. If 
such a program were designed incorporating a high probability of arriving at and 
staying in the “interval” area of Figure | (sector A) and a low probability of staying 
within any single portion of that area, we should expect the resultant behavior to 


“ce 


approximate that appropriate to “variable interval” schedules. In like manner, 
a program which kept the schedule of reinforcement within sector D of Figure 1 
should yield behavior appropriate to a “variable ratio”? schedule though possibly 
not to “random ratio.”” Greater complexity in behavioral properties might result 
from a program not limited to a particular sector of Figure 1. 

We may wish, in addition, to alter the other two restrictions underlying Figure 1. 
For example, it would be possible to present randomly, rather than by regular 
alternations, a fixed ¢? and a fixed (*. Further, it would be possible to specify any 
desired probability of reinforcement of other than the first, or more than one, 


° >) . : . : . 
response in ¢?. In fact, the entire analysis up to now could be duplicated in any 


. . . . ° > « “e ’ és ° . 
?? period by including within ¢? itself any “interval” or “ratio” schedule so far 
discussed. Obviously, too, all the usual independent variables of behavioral 


research, such as drive, could figure as parameters of the present classification. 


SOME EXPERIMENTAL RESULTS 

With the end in mind of experimentally exploring the above schema, an ap- 
paratus was designed and constructed, yielding a range of ¢?/(t? + ¢*) values 
from less than 0.005 to 1.00, with eycle lengths variable from '/, second to ap- 
proximately 4 minutes. The subjects for these experiments were female white 
Carneaux pigeons about six years of age. In order to insure a constant “drive 
level,”’ the birds were maintained throughout the experiment at 80 per cent (+15 
gm.) of their weight under ad libitum feeding. The pigeons were trained to peck 
at a key, or trans-illuminated disk, at one end of the experimental cage; during 
preliminary training, each key-pecking response was reinforced by the presentation 
for 3 seconds of a small grain hopper several inches beneath the key. For a more 
detailed description of the experimental cage and conditioning technique see the 
account of these procedures given by Ferster.6 Response data were taken in the 
form of cumulative response curves, total number of responses and _ reinforce- 
ments per session, and polygraph records of response distributions in time. 
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Two types of studies from which the present data are taken are progressing in 
this laboratory. These data will be presented in greater detail in subsequent 
publications. In one, ¢?/(t? + ¢) is maintained at a constant value while cyele 
length is systematically varied. In the second, part of a doctoral dissertation by 
one of the authors, the length of ¢” is varied while total cycle length is maintained 
constant. On each schedule studied, the original three restrictions apply: t? and 
* are held constant, ¢? and ¢* are alternated, and only the first response in ¢? is 
reinforced. After initial training of the key-pecking response, all birds were 
placed on a schedule permitting the first response in any 30-second period to be 
reinforced (t2/(t? + t°) = 1.00, t? + # = 30 seconds), a schedule identical with a 
30-second ‘fixed interval.” 

After approximately fourteen days of preliminary training on the above inter- 
mittent schedule, birds in the first experiment were shifted to a new schedule, in 
which the total cycle length was 30 seconds and t? was 1.5 seconds (t?/(t? + t*) 
0.05). ¢?/(t? + t*) was held constant at 0.05 throughout this study. Data are 
being taken on cycle lengths of 30, 15, 7.5, 3.75, 1.88, and 0.94 seconds. Each 
bird is kept on a given schedule until its response rate has reached a steady state, 
as defined by the following stability criterion. The first seven days on any schedule 
are not considered in computing stability. For the next six days the mean of the 
first three days of the six is compared with that of the last three days; if the differ- 
ence between these means is less than 5 per cent of the six days’ mean, the bird is 
considered to have stabilized and is shifted to the next schedule. If the difference 
between submeans is greater than 5 per cent of the grand mean, another experi- 
mental day is added and similar calculations are made for that day and the five 
immediately preceding it. Such extensions of the experiment and calculations 
of stability are continued daily until the bird reaches the afore-mentioned 5 per cent 
criterion. Though each bird spends 20 minutes daily, weight permitting, in the 
experimental box, only the last 15 minutes are employed in calculating stability 
criteria. 

Figure 2 displays typical cumulative response curves of bird No. 8 for complete 
20-minute sessions late in training on each schedule with cycle lengths of 30, 15, 
7.5, and 3.75 seconds and respective ¢”’s of 1.5, 0.75, 0.38, and 0.19 seconds. Curves 
are not shown for cycle lengths of 1.88 and 0.94 seconds, since these, as it turned 
out, would be superimposed on the 3.75-second curve. Four birds are being used in 
this study, and all have exhibited consistent functions. With long cycle lengths 
the rate of responding is relatively low; as the cycle length decreases, rate of re- 


sponse increases, appearing to approach an asymptote at about 3.75 seconds. 
“Fixed interval’ responding is displayed when the total cycle length is 30 seconds. 
As the cycle is shortened to approximately 2 seconds, however, the bird’s response 
output develops features characteristic of ‘random ratio” behavior and approaches 


a mean value of 20 responses per reinforcement. 


The second current investigation involves an examination of the effects of de- 
creasing (?/(t? + ¢*) from 1.00 through 0.35, 0.05, 0.03, 0.02, and 0.013, while 
maintaining cycle length constant at 30 seconds. Each bird is used daily, weight 
permitting, for a total of 30 minutes. Figure 3 indicates representative cumulative 
response curves of bird No. 3 on the enumerated ¢?/(t? + ¢) values. When 
t?/(t? + t*) equals 1.00, the cumulative curve is marked by temporal discrimina- 
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Fic. 2.—Typical cumulative re- 
sponse curves for bird No. 8. The 
value of #2? + ¢4 is shown at the 
end of each curve. In all cases 
12/(t? + t4) = 0.05. Reinforce- 
ments are indicated by short diag- 
onal lines. 
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tions typical of responding under “fixed interval” schedules. As the ¢? period 
shortens, the response rates increase and the general shape of the curves becomes 
more similar to that under “fixed ratio” or high mean “variable ratio’? schedules. 
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Fic. 3.—Typical cumulative response curves for bird No. 3. In order to 
conserve space, the curves are displaced to the base line after approximately 
650 responses have accumulated. The value of t?/(t2 + ¢4) is shown for 
each curve. In all cases /? + (4 = 30 seconds. Reinforcements are indicated 
by short diagonal lines. 


F ’ . _ “ : ‘ . ) . . 
‘Breaks” after reinforcements appear to increase with a decrease in ¢?, a finding 


already known to be the outcome of successive increases in the magnitude of the 
fixed ratio. The last curve (#2/(#? + #*) = 0.013) resembles that obtained by 
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Skinner? after shifting subjects from one “fixed ratio” to a much higher “‘fixed 
ratio.” Three other birds employed in this second study show similar changes 
from “‘fixed interval’’-like response curves to “ratio’’-like curves and an increase 
in response rate as (2/(t? + f°) decreases. With continued shortening of e", 
perhaps beyond values used in this study, extinction of the response would even- 
tually occur, since the probability of coincidence between response and t? would 


fall so low (at any actually attainable rate of response) as to make the maintenance 


of behavior impossible. 


SUMMARY 


An attempt was made to show how operant reinforcement schedules may be 
viewed as special cases of one general conceptual framework employing certain 
temporal parameters as its major independent variables. Some preliminary data 
were presented to show how shifts from “interval” to ‘“‘ratio’’ response characteris- 
ties can be obtained by varying these temporal parameters. 
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